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For many years, the Luttinger liquid theory has served as a useful paradigm for the 
description of one-dimensional (ID) quantum fluids in the limit of low energies. This 
theory is based on a linearization of the dispersion relation of the particles constituting 
the fluid. We review the recent progress in understanding ID quantum fluids beyond 
the low-energy limit, where the nonlineajrity of the dispersion relation becomes essen- 
tial. The novel methods which have been developed to tackle such systems combine 
phenomenology built on the ideas of the Fermi edge singularity and the Fermi liquid 
theory, perturbation theory in the interaction strength, and new ways of treating flnite- 
size properties of integrable models. These methods can be applied to a wide variety of 
ID fluids, from ID spin liquids to electrons in quantum wires to cold atoms confined by 
ID traps. We review existing results for various dynamic correlation functions, in partic- 
ular the dynamic structure factor and the spectral function. Moreover, we show how a 
dispersion nonlincarity leads to finite particle lifetimes, and its impact on the transport 
properties of ID systems at finite tcinpcraturcs is discussed. The coiivciitioiial Luttinger 
liquid theory is a special limit of the new theory, and we explain the relation between 
the two. 
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I. INTRODUCTION 

The conventional description of the low-energy prop- 
erties of quantum condensed matter uses the notion of 
quasipartieles: elementary excitations behaving as free 
quantum particles with some energy spectrum which de- 
pends on the microscopic interactions. The low-energy 
properties of interacting electrons in normal metals, for 
example, are well represented by the theory of a Fermi 
liquid (Nozicres, 1997). Its elementary excitations are 
similar to free fermions. One may view the quasiparticle 
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states as those evolving from free fermions when adiabat- 
ically turning on the interactions. 

Quasiparticle states are labeled by their momenta, but 
their dispersion relation (i.e., the quasiparticle energy ^ 
measured from the Fermi level as a function of the mo- 
mentum k) differs from the one for free fermions. These 
states combine the best of the two worlds - on the one 
hand, their overlap with free fermions is significant. On 
the other hand, their lifetimes r become infinitely long as 
their energy ^ vanishes: in the absence of static disorder, 
T (X An electron easily tunnels into a metal, "dress- 
ing up" in the tunneling process to become a quasipar- 
ticle. Neglecting the small relaxation rate 1/r, an elec- 
tron entering a Fermi liquid with momentum k creates a 
single quasiparticle with a well-defined energy ^(k). In 
an inverse process, an electron may tunnel out, leaving 
behind a hole with well-defined energy. In either case, 
the tunneling probability per unit time of an electron 
with given momentum k and energy e or, more precisely, 
its spectral function A{k,e) is close to a delta-function, 
^(k, e) oc d{e — ^(k)). Residual interactions between the 
quasiparticles can be readily accounted for within con- 
ventional perturbation theory. A perturbative evaluation 
of the quasiparticle 's self-energy leads to a finite relax- 
ation rate 1 /r and to a slight broadening of the spectral 
function, transforming it into a Lorentzian with a width 
(5£~ 1/r. 

If the momentum is not conserved in the tunneling 
event (this happens, for example, if the electron tunnels 
through a point contact), then the electron extracted 
from the Fermi liquid leaves behind a superposition of 
holes, each of them having the same energy e. The 
area of the constant-energy surface in momentum space 
defined by the equation ^(k) = e determines the den- 
sity of states i/(e) oc / dkA{k,e) available for tunnel- 
ing (Nozieres, 1997) at energy e. Similar to the case of 
free fermions, the tunneling density of states at the Fermi 
level v{0) in the Fermi liquid is finite. 

One can find a more subtle example of a quasiparti- 
cle description in the Bogoliubov treatment of the ex- 
citations of a Bose gas with weak inter-particle repul- 
sion (Pitaevskii and Stringari, 2003). Bogoliubov quasi- 
particles are the bosonic excitations above the Bose con- 
densate; their spectrum differs qualitatively from the 
spectrum of "bare" bosons: due to the interaction be- 
tween the "bare" particles, the Bogoliubov quasiparticles 
are characterized by a sound-like spectrum, ^(k) (x k, at 
low energies. Still, a particle entering the system of inter- 
acting bosons easily "dresses up" to become a quasipar- 
ticle. The spectral function A(k, e) of a boson with mo- 
mentum k is close to a delta-function centered at energy 
e = C(k). At small energies, the lifetime of a Bogoliubov 
quasiparticle (Beliaev, 1958) with energy ^ diverges as 
T oc 1/^^, making the quasiparticle states well-defined. 
The rate 1/r determines the broadening of the spectral 
function. 



In the above examples, the affinity between the free 
particle and an elementary excitation of the many-body 
system is exemplified by the narrow energy width of the 
spectral fimction ^(k,e): the spectral weight is concen- 
trated around the quasiparticle energy £ = ^ (k) within a 
region 5s <C ^(k). While working well in higher dimen- 
sions, this picture fails in the case of a one-dimensional 
(ID) gas of quantum particles. This is clearly visi- 
ble, e.g., for ID fermions in the presence of even a 
weak interaction between them. The correction to the 
fermionic spectral function to the second (lowest non- 
vanishing) order in the interaction potential decays as 
l/[e — ^{k)] (Dzyaloshinskii and Larkin, 1974), transfer- 
ring the spectral weight away from the quasiparticle mass 
shell s = ^(fc). Along with the slow decay of the spectral 
function, the second-order corrections to the tunneling 
density of states v{e) at e = and to the momentum 
distribution function n(k) oc deA{k, e) at the Fermi 
points k = ikp, arc singular. 

The "magic bullet" effective in resolving many of the 
difficulties of the ID quantum many-body problem was 
suggested in a seminal paper by Tomonaga (Tomonaga, 
1950). It was noticed there that replacing the generic 
dispersion relation ^{k) of ID fermions with a linear 
one, ^(fc) = ±VF{k ^ kp), immensely simplifies finding 
the many-body dynamics of the system (here the up- 
per/lower signs correspond to the right-/left-moving par- 
ticles, and kp is the Fermi momentum). For free fermions 
with a linear spectrum, the energy E oi a right-moving 
excitation consisting of an arbitrary number of particle- 
hole pairs with a given total momentum q depends only 
on that total momentum, E = Vpq- This degeneracy al- 
lowed Tomonaga to encode the excitations of the Fermi 
gas (at fixed numbers of left- and right-movers) into the 
excitation spectrum of free bosons. These ID bosons are 
nothing but quantized waves of density of a ID Fermi 
gas, and their description is identical to that of acoustic 
phonons (Ziman, 1960). The beauty of the encoding is 
that it puts the free-fermion Hamiltonian and density- 
density interactions on an equal footing. Indeed, the full 
Hamiltonian of interacting ID fermions now becomes a 
bilinear form in bosonic creation-annihilation operators. 
Its diagonalization is standard and not different from 
the corresponding procedure for phonons (Ziman, 1960). 
This way, the Hamiltonian for interacting ID fermions 
with linear spectrum is diagonalized and cast in terms of 
free bosons with linear dispersion relation, u}{q) = v\q\ 
(we use units h= 1 throughout the text); the velocity v 
differs from vp because of the interactions. 

The bosonic representation introduced by Tomonaga 
makes the calculation of the density correlation func- 
tion of ID fermions straightforward. Indeed, the den- 
sity ffuctuation operator p{x, t) is linear in boson cre- 
ation/annihilation operators. According to the diagonal- 
ized form of the Hamiltonian, those represent excitations 
which propagate freely with a constant velocity v. This 
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is reflected in the dynamic structure factor (DSF), which 
is defined as the probabiUty per unit time to excite a den- 
sity fluctuation by an external source coupled to p. The 
DSF takes the form S{q,uj) cx |(j|(5(ci; — v\q\). 

The evaluation of the propagator of a fermion, its 
spectral function A{k,e), the tunneling density of states 
i^(e), or the distribution fimction n(fc) is somewhat more 
complicated. Luttinger (1963) attempted to evaluate 
the distribution function using the fermionic representa- 
tion; some aspects of that calculation were clarified later 
by Mattis and Lieb (1965). The spin-1/2 fermion prop- 
agator G{x,t) in space and time domain was evaluated 
for the Tomonaga-Luttinger model by Dzyaloshinskii and 
Larkin (1974). Their diagrammatic technique heavily re- 
lied on the linear dispersion relation of fermions. Luther 
and Peschel (1974) evaluated G{x, t) with the help of 
the bosonic representation of the Tomonaga-Luttinger 
model (Luther and Peschel, 1974; Mattis, 1974), and also 
showed that the tunneling density of states z/(e) displays 
a power-law behavior at energies e much smaller that the 
Fermi energy. Similar techniques have been used to cal- 
culate the long-distance behavior of correlation functions 
for ID bosonic systems (Efetov and Larkin, 1975). 

At low energies, the excitations of noninteracting 
fermions are particles and holes with momenta k in the 
vicinity of ifej?. The energy of, say, a right-mover is 
^(fc) = vpik - kp) + (k - kpf/i^m*). The second 
term here may be considered as the lowest-order expan- 
sion of a general nonlinearity of the dispersion relation 
(in case of particle-hole symmetry, such an expansion 
would start from the cubic (fc — kp)'^ term). For nonin- 
teracting fermions with Galilean-invariant spectrum, m* 
is equal to the bare mass. The quadratic term here 
scales as ^^/(to*w|,) at small ^: by a power-counting 
argument, this term is irrelevant. Indeed, it has been 
shown by Haldane (1981b) that the spectrum nonlinear- 
ity does not significantly affect the long-range behavior 
of the fermion propagator G(a;, 0) at fixed time (t = 0). 
That gives an incentive to dispense with the nonlinear- 
ity of the dispersion relation. After that, a wide vari- 
ety of ID systems with gapless excitation spectra can 
be mapped, at low energies, on the Tomonaga-Luttinger 
model. Thus, the Tomonaga-Luttinger model provided 
the foundation for the concept of the Luttinger liquid 
(LL), a phenomenological description of the low-energy 
excitations of interacting quantum particles confined to 
ID (Haldane, 1981a,b). In addition to liquids of fermions 
or bosons, these include also the low-energy excitations of 
half-integer spin chains, i.e., spin liquids (Haldane, 1980). 

After a real system is replaced by the corresponding 
LL, the "good" low-energy excitations appear to be waves 
of density of the corresponding liquid which have a lin- 
ear spectrum w(g) = v\q\. These excitations propagate 
along the a;-axis with fixed velocities ±.v and without any 
dispersion: a perturbation created at some point x propa- 
gates without changing its shape to the points x±vt. The 



presence of the formally irrelevant term {k — kp)'^ /{2m*) 
in the spectrum results in a dispersion of the propagating 
perturbation. In analogy with single-particle quantum 
mechanics, the width of the perturbation grows with time 
(X sjtjm* . The irrelevance of the quadratic term in ^(A;) 
means that the growth is slow compared to the rate of lin- 
ear displacement of the perturbation: yJTJm* / (vt) 
at f — >■ oo. The infinitely "sharp" DSF of the Tomonaga- 
Luttinger model, S{q,oj) oc |g'|5(a; — v\q\), reflects the 
propagation of the density perturbation with fixed veloc- 
ities ±v. The (X ^Jtjm* dispersion of the perturbation 
corresponds to some kind of broadening of the DSF: at 
given l^l « Lojv the characteristic width of the DSF is 

It may look like the designation of density waves in ID 

as the proper excitations parallels the Fermi-liquid idea 
in higher dimensions. Indeed, dispensing with irrelevant 
perturbations in these two systems, one finds a delta- 
fimction structure of 5'(g, uj) for ID bosons and similarly 
of >l(k, e) for Fermi quasiparticles in higher dimensions. 
Prom a dimensional analysis, one expects the irrelevant 
terms to broaden these delta- fimctions by cx o;^ and e^, 
respectively, at given q and k. However, the similarity 
stops there. In a Fermi liquid, irrelevant interactions lead 
to a self-energy with small imaginary part in the quasi- 
particle Green function, which transforms A(k, e) to a 
Lorentzian. The self-energy can then be evaluated within 
perturbation theory in the irrelevant terms (Abrikosov 
et al., 1963). In contrast, a naive attempt to use pertur- 
bation theory to evaluate the self-energy of bosons and 
thus "broaden up" S(q, lj) in the LL is doomed (Andrcev, 
1980; Pirooznia and Kopietz, 2007; Punk and Zwerger, 
2006; Samokhin, 1998). The symmetry (in fact, Lorentz- 
invariance) introduced by the linearization results in a 
degeneracy of the excitation spectrum. The terms de- 
scribing the curvature of the dispersion relation break 
that symmetry. The very same simplification which al- 
lowed Tomonaga (1950) to find the exact solution of the 
problem makes the perturbation theory in the curvature 
strongly degenerate. The nominally irrelevant terms re- 
move the degeneracy, lead to the emergence of new quali- 
tative behavior of the dynamic correlation functions, and 
to relaxation processes which do not exist in the linear 
LL. 

The perturbation theory in curvature is plagued by on- 

shell divergencies, and is in need of a proper resummation 
procedure, yet to be developed (Aristov, 2007; Samokhin, 
1998). A possible way of building a perturbation the- 
ory in boson representation may involve studying the re- 
sponses in an off-shell domain and comparing the results 
to the known limits, such as the free fermions (Pereira 
et al, 2007; Pereira et al, 2006; Teber, 2006, 2007). 

A progress, however, was achieved along a different 
route, via linking the properties of a nonlinear LL to the 
well-known Fermi edge singularity effect (Mahan, 1981; 
Nozieres and De Dominicis, 1969; Ohtaka and Tanabe, 
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1990). That connection was quite easy to notice for 
weakly-interacting fermions with a generic (nonlinear) 
spectrum (Pustilnik et al., 2006). Further development 
of that relation put the problem of edge singularities 
in a nonlinear LL into the class of so-called "quantum 
impurity" models (Affleck, 2009; Balents, 2000; Castella 
and Zotos, 1993; Castro Ncto and Caldcira, 1994; Cas- 
tro Neto and Fisher, 1996; Cheianov and Pustilnik, 2008; 
Imambekov and Glazman, 2008, 2009a,b; Khodas et al., 
2007a,b; Lamacraft, 2008, 2009; Pcreira and Sela, 2010; 
Pereira et al, 2008, 2009; Schmidt et al, 2010a,b; Sorella 
and Parola, 1996; Tsukamoto et al., 1998; Zvonarev et al, 
2009a,b). 

The use of the field-theoretical approach based on 
quantum impurity models has led eventually to a new 
phcnomcnological theory of nonlinear LLs. Remarkably, 
the threshold power-law singularities in the dynamic re- 
sponses [S{q,uj) and A{k,e) are the examples] occur not 
only in the vicinity of special points (k = ztkp for the 
spectral function), but at arbitrary momenta. The non- 
linear LL phenomenology relates the exponents charac- 
terizing the singularities to the properties of the threshold 
spectra, thus establishing a relation between two sets of 
independently measurable quantities. 

The phenomenology also provides effective tools for the 
evaluation of the exponents characterizing the dynamic 
responses: one may find first the energy spectra and then 
use the phcnomcnological relations to find the exponents. 
A class of systems for which such a program is especially 
attractive are the models exactly solvable by the Bethe 
ansatz. The thermodynamic Bethe ansatz is well suited 
for the evaluation of the spectra but not for the dynamic 
responses. However, the corresponding exponents can be 
found exactly by combining the thermodynamic Bethe 
ansatz with a field-theoretical description (Cheianov and 
Pustilnik, 2008; Essler, 2010; Imambekov and Glazman, 
2008; Pereira et al, 2008, 2009). Exactly solvable models 
provide stringent nontrivial tests of the field-theoretical 
approaches. 

The linear LL theory does not discriminate between 
integrable and generic ID systems: in either case, the 
original system is replaced by free particles devoid of 
any relaxation mechanisms. Methods emerging within 
the nonlinear theory pave a way of studying the kinetics 
of a ID quantum liquid and see the differences between 
generic and integrable systems. 

In summary, understanding ID quantum liquids out- 
side the sector of low-energy excitations requires break- 
ing the spell of linearization. The emerging theory, which 
accounts for the nonlinear energy spectrum of particles 
forming the liquid, answers that challenge. We review 
a number of methods of the nonlinear LL theory and 
expose relations between them. These methods have al- 
ready led to a progress in understanding the dynamic 
responses and relaxation of ID quantum liquids. The 
approaches we describe are controllable, yet versatile 



enough for application to a broad class of systems, from 
electrons in quantum wires (Deshpande et al., 2010), to 
spin liquids and cold atoms confined to ID. 

This review is organized as follows. In Sec. II, we in- 
troduce and develop in detail the general field-theoretical 
approach to describe the singularities of dynamic re- 
sponse functions in the momentum-energy plane. That 
approach is based on the phenomenology of effective mo- 
bile impurities moving in LLs. In Sec. Ill we combine this 
ficld-theorctical approach with the analysis of exactly 
solvable models, which allows one to obtain a plethora 
of new results for the latter, and provides stringent non- 
perturbativc checks of the phcnomcnological approaches. 
In Sec. IV, we illustrate the importance of the physics be- 
yond the linear LL theory for the kinetics and the trans- 
port in ID quantum liquids. 



II. SINGULARITIES OF THE DYNAMIC RESPONSE 
FUNCTIONS 

An adequate description of a quantum many-body sys- 
tem not only requires an understanding of the ground 
state, but also a characterization of its excitations. One 
of the most natural ways to probe the excitation spec- 
trum is to measure the dynamic responses of the sys- 
tem to external fields, such as electromagnetic radia- 
tion of a given momentum and energy. Within Fermi's 
golden rule, the scattering rate of such external fields 
is related to various dynamic response functions, such 
as the DSF or the spectral function (see precise defi- 
nitions below). This motivates the interest in study- 
ing qualitative features of the dynamic response func- 
tions of low-dimensional systems. In particular, we will 
mostly be interested in their behavior near the spectrum 
of elementary excitations. Many experimental techniques 
can be applied to probe the dynamic responses of ID 
systems, such as neutron scattering (Lake et al., 2005, 
2010; Masuda et al, 2006; Nagler et al, 1991; Riiegg 
et al, 2005; Stone et al, 2003; Tcnnant, 2009; Thiele- 
mann et al, 2009; Zheludev et al, 2008, 2002), angle 
resolved photoemission spectroscopy (ARPES) (Blumen- 
stcin et al, 2011; Claessen et al, 2002; Hoinkis et al, 
2005; Kim et al, 2006; Kondo et al, 2010; Sing et al, 
2003; Wang et al, 2009, 2006), and various forms of 
Bragg spectroscopy (Clement et al, 2009; Ernst et al, 
2010; Fabbri et al, 2009; Stamper-Kurn et al, 1999) and 
photoemission spectroscopy (Dao et al, 2007; Gaebler 
et al, 2010; Stewart et al, 2008). Additional interest 
in the response functions of ID systems is driven by 
a rapid progress in their numerical evaluation (Barthel 
et al, 2009; Feiguin and Huse, 2009; Kohno, 2010; Kokalj 
and Prelovsek, 2009; White and Affleck, 2008) based on 
time-dependent extensions of Density Matrix Renormal- 
ization Group (DMRG) techniques (De Chiara et al, 
2008; Schollwock, 2005; Vidal, 2003, 2004; White, 1992). 
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A nonlinear dispersion relation and interactions be- 
tween particles forming a ID quantum liquid modify in a 
nontrivial way all dynamical responses of the liquid, re- 
solved in energies and momenta. The DSF provides one 
of the examples. It is defined as 



S{q,uj) = 



dt 



dxe^^'^'-'^^\p{x,t)p{0,0)). (1) 



Here, p(x, t) is the density operator and the averag- 
ing (...) is performed over the Gibbs ensemble or the 
ground state in the case of finite or zero temperature, 
respectively. In the Tomonaga-Luttinger model, at small 
wavevectors q the DSF takes the form S'll oc |g|i5(a;— 
at any temperature. The dispersion results in a "broad- 
ening" of the delta-function. Accounting for a finite 
width of S{q, uj) even for small q is important (Pustil- 
nik et ai, 2003) for understanding Coulomb drag exper- 
iments (Debray et al, 2001, 2002; Laroche et ai, 2011; 
Yamamoto et ai, 2002, 2006). The broadening occurs 
even at zero temperature (T = 0), and we will concen- 
trate on that case. 

To illustrate the origin of the structure arising in 
S{q,uj) due to the dispersion, let us consider first the 
simplest case of free spinless fermions with a quadratic 
dispersion relation. 



2m 



(2) 



At zero temperature, the structure factor can be thought 
of as an absorption coefficient, i.e., the dissipative part 
of the linear susceptibility with respect to a perturba- 
tion 5% = U(x, t)p{x, t) by a potential U(x, t) varying in 
space and time with the wave vector q and frequency cj, 
respectively. In the case of free fermions, dissipation is 
caused by creation of particle-hole pairs by the perturb- 
ing potential, see Fig. 1. At q < 2kp , a. simple evaluation 
of Eq. (1) yields 

So{q, Lo) = (m/g)^?(gV(2m) - |^ - vpql) (3) 

with the Fermi velocity vp = kp/Tn. The two thresholds 
for absorption correspond to two special configurations in 
the momentum space of the particle-hole pairs, see Fig. 1. 
Specifically, the lower boundary, uj-{q) = vpq — q^ l{2m), 
corresponds to a particle just above the Fermi level, and a 
"deep" hole with momentum kp — q, moving with velocity 
Vh = vp — q/m, smaller than vp. Equation (3) and Fig. I 
allow us to make three interesting observations: 

First, at fixed q <^ 2kp, the width of the structure 
factor in the frequency domain Sco ~ uP' /(mvy) scales as 
(X up' . This is consistent with the power-counting argu- 
ment for the irrelevant curvature term in the spectrum 
S,{k). The limit of linear spectrum (which is also a trivial 
limit of an LL at zero inter-particle interaction) corre- 
sponds to taking m — ^ cx) at fixed value of vp. In this 
limit, (5w — 0. 





hp/ 




_^kp' + q 






hp 




/ki-- 
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FIG. 1 (Color online) Density structure factor S[q, u) for 
noninteracting fermions [see Eq. (3)] and weakly interacting 
fermions [see Eqs. (27), (30) and (32)]. In the noninteracting 
case, So{q,ijj) is constant in the dark shaded region and van- 
ishes otherwise. Left: Particle-hole configurations responsi- 
ble for the upper and lower thresholds. Right: S{q,ij) for 
fixed < q < 2kF. In the noninteracting case, S{q,uj) has a 
rectangular shape. Interactions turn the steps into power-law 
singularities, and S{q, lj) becomes nonzero in the light shaded 
region. 



Second, the dependence of So{q, cj) on its arguments is 
not analytic; besides, if one allows an arbitrary sign of the 
mass TO, it becomes immediately clear that Su oc 1/|to| 
at a fixed value of the Fermi velocity vp. Each of these 
two facts kills the hope to develop a simple perturbation 
theory in the irrelevant perturbation, i.e., the curvature 
of the particles dispersion relation. Indeed, it is clear 
from the above discussion that the perturbing part of 
the Hamiltonian associated with the curvature is propor- 
tional to 1/to. To obtain a broadened structure factor, 
one needs to evaluate the self-energy of the density prop- 
agator. The form of Eq. (3) tells us that the imaginary 
part of self-energy is cx 1/|to| and may result only from 
a summation of some infinite series in 1/m. 

Third, the specific structure of the particle-hole pair 
corresponding to the edge w_(g) gives us a hint at how 
weak interactions may modify Eq. (3). Indeed, suppose 
fermions weakly repel each other. Then, the created par- 
ticle would be attracted to the hole it left upon excitation. 
This interaction would lead to the Mahan (excitonic) 
singularity in the absorption coefficient (Mahan, 1981). 
That is, the step-like threshold at uj — uj-{q) would trans- 
form into a divergent power-law function with an expo- 
nent dependent on the inter-particle repulsion strength. 

The above simple picture establishes the relation, 
which is central for this section, of the nonlinear LL prob- 
lem to the well-studied problem of Fermi-edge singular- 
ities (Nozieres and De Dominicis, 1969). The latter is 
reviewed in great detail elsewhere (Gogolin et ai, 1998; 
Mahan, 1981; Ohtaka and Tanabe, 1990). We wih see 
that the power-law asymptote of the DSF at energies 
close to the threshold is a robust feature valid at arbitrary 
interaction strength and arbitrary q. We will also develop 
ways to evaluate the corresponding exponent, which does 
depend on these parameters. 

A nonlinear dispersion relation of interacting quantum 
particles confined to ID affects also their spectral func- 
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FIG. 2 (Color online) Spectral function A{k ~ kF,e) for 
a Luttinger liquid (LL) with linearized spectrum. Left: 
A{k, e) vanishes in the white regions. Right: A{k, e) for 
fixed k ^ kp- Interactions cause a power-law divergence 
at the mass shell e ~ v{k — kp) with exponent /i_ = 
1 — (JC + K^^ — 2) /4. A convergent power-law cusp with 
exponent /j,^ = — (A' + K^^ — 2) /4 emerges at the inverted 
mass shell e ~ —v{k — kp)- 




FIG. 3 (Color online) (a) An incoming particle with momen- 
tum k ~ kp and energy above the mass shell can tunnel into 
the system by creating a particle with momentum near kp 
on the mass shell and a low-energy particle-hole pair near the 
opposite Fermi point, (b) If the spectrum is curved, an in- 
coming particle with momentum k ^ kp and energy below the 
mass shell can tunnel into the system by creating a particle on 
mass shell and a low-energy particle-hole pair near the same 
Fermi point. 



tion A{k,e). The latter is defined as 



A{k,e) 



1 



IniG'(A;, e)signe 



(4) 



with the Green's function (Abrikosov et al.^ 1963) 

/oo poo 
-OO J — OO 

(5) 

Here ^{x,t) and '^\x,t) are the particle (fermion or bo- 
son) annihilation and creation operators, respectively, T 
denotes the time ordering, and the energy e is measured 
from the chemical potential. 

The spectral function may be thought of as a tunnel- 
ing density of states: the probability for a particle (hole) 
with given momentum k and energy e > (e < 0) to 
tunnel into a system is proportional to A{k, e). It can be 
measured using ARPES (Kim et al., 2006; Kondo et ai, 
2010; Wang et al, 2009, 2006) for sohd state systems or 
photoemission spectroscopy (Gaebler et al, 2010; Stew- 
art et ai, 2008) for low-dimensional ultracold atomic sys- 
tems (Gorlitz et at, 2001; Kinoshita et al, 2004, 2006; 
Paredes et al, 2004). It also determines electronic trans- 
port in systems with momentum and energy conserving 
tunneling (Auslaender et al, 2005, 2002; Barak et al, 
2010; Jompol et al, 2009). In the absence of interactions, 
a particle with a given momentum may tunnel only if its 
energy fits the dispersion relation of the particles consti- 
tuting the system, Ao{k,e) = 6[e — ^{k)]. The right-hand 
side here is the density of single-particle eigenstates with 
given energy and momentum. 

Before considering the effects of a nonlinear disper- 
sion, let us recall the behavior of A{k, e) in a fermionic 
LL. In the absence of interactions, the tunneling den- 
sity of states for, say, right-movers is A{k,e) = S[e — 
v{k — kp)]. Interactions between the particles forming 
the LL broaden the spectrum of energies at which tun- 
neling is possible. In the vicinity of the Fermi point 



+kp one has (Luther and Peschel, 1974; Meden and 
Schonhammer, 1992; Voit, 1993a,b, 1995) 



A{k, e) oc sign(e)- 



e — v{k — kp) 



X [e2_^;2(fc-fc^)2]K^+*)-^ . (6) 

The shape of A{k,e) for a linear LL is shown in Fig. 2. 
Here, the LL parameter K depends on the interaction 
strength [K < 1 for repulsion); free fermions correspond 
to the limit A' — >■ 1. The delta-function in the tunneling 
density of states of free particles got transformed into 
a power-law, divergent at the line e = v{k — kp) if the 
interaction is not too strong, see Eq. (6). Note that an 
important feature of the linear LL result is the particle- 
hole symmetry under the transformation 



(k^kp) -{k 



(7) 



which is a necessary consequence of the spectrum lin- 
earization. 

The nonzero values of A{k, e) outside the line e = v{k— 
kp) may be understood within the perturbation theory 
if one invokes the notion of tunneling probability. In the 
presence of interactions between particles, the tunneling 
free particle with energy e > v\k — kp\ may spend the 
extra energy on the creation of a particle-hole pair on the 
branch of left-movers and land on the mass shell £,{k) = 
v{k — kp) for the right-moving particles. This process is 
depicted in Fig. 3a. Similarly, the tunneling of a particle 
at energy — e > v\k ~ kp\ out of the system creates a 
right-moving hole and left-moving particle-hole pair. To 
summarize, in this perturbative picture tunneling of a 
particle into (out of) the system creates three excitations: 
a right-moving particle (hole) and a particle-hole pair on 
the opposite branch. 
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The above pcrturbativc consideration is easy to carry 
over to fermions with the nonlinear dispersion rela- 
tion (2). We see immediately the difference between the 
particle-likc {k > kp) and hole-like parts of the spec- 
trum. The hole-like part of the spectrum becomes the 
energy threshold ^(fc) < for tunneling at < kp- In- 
deed, because the hole's velocity \k/m\ < vp, it is not 
capable of emitting "Cherenkov radiation" of low-energy 
particle-hole pairs in either of the two allowed directions. 
At energies e < ^(A;) < 0, the tunneling of a hole is ac- 
companied by the creation of particle-hole pairs at both 
Fermi points. 

However, as shown in Fig. 3, the particle-like part 
of the spectrum (2) falls into a continuum of energies 
available for tunneling at given \k\ > kp. To see this, 
we may again concentrate on the timneling of a right- 
moving free particle with momentum k > kp. Unlike in 
the case of a linear spectrum, now a particle with en- 
ergy < £ < ^(fc) may tunnel by creating a comoving 
particle and a particle-hole pair. The total of three ex- 
citations should have the momentum k, but the sum of 
their energies is less than ^(fc) if all three momenta are 
within the region of width k — kp around kp. Within 
perturbation theory, we find the threshold in this region 
at — ^(fc — 2kp) for kp < k < 3kp, see Fig. 3b. Con- 
sequently, states at the free-particle mass shell ^(fc) at 
k > kp are not protected by kinematics: particles move 
fast enough (|fc/m| > vp) to allow Cherenkov radiation 
of particle-hole pairs. 

A comparison of the perturbative pictures for the linear 
and nonlinear dispersion relations reveals some substan- 
tial ramifications introduced by the nonlinearity. The 
nonlinearity destroys the particle- hole symmetry, Eq. (7), 
which existed in the LL. The hole-like part of the thresh- 
old morphs from a straight line into a nonlinear function; 
the nature of excitations created by a tunneling hole is 
not changed by the introduction of the curvature. How- 
ever, the particle-like part of the threshold changes dras- 
tically; the excitations defining it have no counterpart 
in the linear LL. In Sec. II. B we will review a univer- 
sal nonlinear LL theory valid in the vicinities of Fermi 
points (Imambekov and Glazman, 2009b). 

We should emphasize that at fc ztkp the range of 
energies in which A{k, e) is substantially modified com- 
pared to linear LL theory becomes narrow, since it scales 
as |£ — £,{k)\ oc {k — kp)^/m. At energies \s — ^{k)\ >• 
{k — kp)'^/m the linear LL theory does indeed describe 
the structure of spectral function. This is consistent with 
the curvature being an irrelevant perturbation (Haldane, 
1981b). However, the true threshold behavior of the spec- 
tral function is controlled by the nonlinear spectrum at 
any wavevector fc, close to or far away from ±kp. 

The above discussion based on the perturbation theory 
tells us that the presence of thresholds in the dynamic 
response functions is protected by kinematics: a slow- 
moving excitation (a hole in the above consideration) can 



not decay in ID. An inspection of Figs. 1 and 3b may raise 
a question, whether the response functions of a nonlinear 
LL have singularities within the spectral continuum. The 
answer is: "it depends" . In the case of integrable models, 
these singularities do survive, while in a generic liquid 
they are broadened and get progressively washed out if 
one moves away from the Fermi points. We discuss the 
singularities in the continuum in Sec. HI devoted to the 
integrable models, and the broadening of singularities by 
relaxation processes in Sec. IV devoted to the kinetics of 
nonlinear LLs. 

The remainder of this section is organized as follows. 
In Sec. II. A we make the above perturbative consid- 
erations quantitative and derive the effective mobile- 
impurity Hamiltonian for weakly interacting spinless 
fermions. The perturbation theory gives a clear hint on 
how to proceed with the calculation of the threshold sin- 
gularities at arbitrary interaction strength. In Sec. II. B 
we explain the theory of a nonlinear LL at arbitrary in- 
teractions in the vicinity of the points k = ±kp, consider- 
ing in detail the crossover between the generic threshold 
behavior and the linear LL asymptote. The adequate 
apparatus based on a mobile quantum impurity moving 
in a linear LL is extended further, and we develop the 
phenomenology of the threshold behavior of the dynamic 
responses for spinless fermions (Sec. II. C), spin liquids 
(Sec. II. D), bosonic systems (Sec. II. E), as well as for 
spinful fermionic systems (Sec. II. F). Next, Sec. II. G is 
devoted to effects which arise due to a finite system size 
and finite temperatures. Finally, in Sec. II. H we discuss 
the implications of the threshold singularities for corre- 
lation functions in the space-time domain and for the 
breakdown of conformal invariance. 

A. Perturbative treatment of interactions 

1. Dynamic structure factor (DSF) 

The DSF S{q, lu) characterizes the linear response of 
the density to an external field which couples to the den- 
sity. The absorption rate for quanta of momentum q and 
energy uj of such a field {e.g., photons) is proportional to 
S{q,Lj). The rate is given by Fermi's golden rule, 

%,^) = 5EI(/Kio)I''^(^-^/)- (8) 

l/> 

The equivalence of Eqs. (1) and (8) can be demonstrated 

using the Lehmann spectral representation. Here, the 
system length is denoted by L, and the Fourier compo- 
nents of the density are defined by = J2k ^k+q^k- 

In noninteracting systems, the groimd state |0) at 
zero temperature consists of a Fermi sea filled up to the 
Fermi momentum kp. A nonzero matrix element (/|pj|0) 
emerges only for final states |/) which contain exactly 
one particle-hole pair with momentum q. All possible 
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final states can thus be parameterized by the particle 
momentum kp and the hole momentum kh, 



l/) = *L*fe.|o), 



(9) 



where \kp\ > kp, \kh\ < kp, and kp — k^ = g. For the 
quadratic spectrum (2), the energy of such a state is e/ = 
{kp — k\)/{2m). The evaluation of the DSF according to 
Eq. (8) is then straightforward and yields for q < 2kF, 



m 



Suj{q) = w+{q) - = — . 

m 



(11) 



In view of a perturbative analysis, it is convenient to 
express the noninteracting DSF in terms of fermionic 
Green's functions. Using the fluctuation-dissipation the- 
orem, the DSF at zero temperature can be related to the 
susceptibility (Doniach and Sondheimer, 1998), S{q,co > 
0) = — 2Imx(5, w), where xil^'^) is the Fourier trans- 
form of 



x{x,t) = -ie{t) {[p{x,t),p{o,o)]) . 



(12) 



For Lo > 0, the imaginary part of the retarded density- 
density correlation function xil^'^) coincides with the 
imaginary part of the polarization diagram, 



^o(9,w) 



(13) 



where solid lines denote time-ordered fermion Green's 
functions, and the internal momentum qi and the internal 
energy cji arc integrated over. Therefore, for the nonin- 
teracting system So{q,LO > 0) = — 2ImPo(?)W), yielding 
again Eq. (10). 

Next, let us calculate the correction to the DSF for a 
weak density-density interaction. 

Hint = \ j dxdyp{x)V{x - y)p{y). (14) 



The first-order term 5S'^^^ {q, lo) can be cast into an RPA- 
like diagram and a vertex correction. 



5S^^\q,Lo) oclm 



+ lm 



(15) 



So(«,w) = j^^e[uj-uj.{q)]e[uj+{q)-w]. (10) 

It turns out that S'o((7,cj) is nonzero only within the in- 
terval LU- (q) < ui < {q) . The upper and lower edges of 
support physically correspond to final states |/) which 
have the highest and lowest possible energies, respec- 
tively, for a given momentum q. As shown in Fig. 1, 
for q > the final state with the maximum energy con- 
tains a hole with momentum kh = kp and a particle 
with momentum kp = kp + q. The energy of this state is 
co+{q) = Vpq + q^/{2m). Similarly, the density excitation 
with the minimum energy for q > contains a hole with 
momentum kh = kp — q and a particle near the Fermi 
point, kp = kp. This state has the lower threshold en- {li^) 
ergy ijJ-{q) = vpq — q^ /{2m). The width of support of 
SQ{q,ui) for fixed q is therefore 



The wiggly lines denote the interaction. Let us start with 

the first diagram and estimate its contribution towards 
the lower edge of support, uj « Lo^{q). Its imaginary 
part is proportional to VqlmVo{q,uj)ReVo{q,uj). How- 
ever, we showed that lTiiVQ{q,uj) is proportional to the 
noninteracting DSF, so lmVo{q,'-^) oc 9[lu — Lu^{q)] con- 
tains a threshold. The Kramers-Kronig relation (Landau 
and Lifshitz, 1980) therefore predicts a logarithm in the 
real part, Ke'Po{q,oj) oc ln{[cj — uj-{q)]/Suj{q)}. So, the 
first diagram diverges logarithmically towards the edge. 
The second diagram can also be calculated and leads to 
an identical asymptote but with a prefactor — Vq. As a 
result, the total first order correction near the threshold 
reads 



rn'(^g-Vb) 
7r|5|2 



.{q)]\n 



Suj{q) 



(16) 



Hence, a straightforward calculation of the first-order 
correction SS^^''{q,uj) leads to a logarithmic threshold 
divergence. The underlying physical mechanism is rem- 
iniscent of the Fermi edge singularity problem (Ander- 
son, 1967; Gogolin et al, 1998; Mahan, 1967, 1981; 
Nozieres and De Dominicis, 1969; Ohtaka and Tanabe, 
1990; Schotte and Schotte, 1969): the final state |/) for 
w = a;_ {q) and < g < 2kp contains a hole at momen- 
tum kp — q which generates a scattering potential. The 
infrared divergence is produced by scattering of particles 
near the Fermi points with small momentum exchange. 
In order to obtain a viable result, a partial resummation 
of the perturbation series is needed. 

The analogy between the physics of the Fermi edge 
singularity and the threshold behavior of the DSF can 
be exploited (Pustilnik et al., 2006) by using a method 
which is familiar from the solution of the Fermi edge 
singularity problem by Schotte and Schotte (1969). In 
the calculation of S{q,oj) near the threshold co w uj-{q), 
the hole with momentum kp — q assumes the role of a 
"deep" hole. The energy left for additional excitations 
\uj — oj-{q)\ is small, so all particle-hole pairs created due 
to the interaction are restricted to a small window of 
width /co <C g around the Fermi points. The Hamilto- 
nian can then be projected onto three subbands of width 
fco , one centered around kp — q containing the deep hole 
and two containing the Fermi points ±kp. We note that 
this projection into subbands is in fact very similar to a 
conventional procedure employed in the Fermi edge sin- 
gularity problem. In the latter, a full fermionic operator 
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is split into contributions from non-overlapping conduc- 
tion and core-hole bands. In our case, all fermions are in 
the same band, but since the momenta of the important 
states do not overlap due to kinematic constraints, split- 
ting the fermionic operator into subbands is a legitimate 
procedure. In contrast to the original Fermi edge prob- 
lem, the deep hole is mobile but this does not destroy the 
edge singularity (Balents, 2000; Ogawa et al., 1992). 

Let us illustrate the calculation for < q < 2kF- In 
this case, the fermion annihilation operator ^{x) is pro- 
jected onto the three subbands using 



*(a;) 



^fl(x) + e-*'=-^VL(2;) + e*'=^d(:r), (17) 



where k = kp — q lies within the Fermi sea and the pro- 
jected operators d{x) and ipR^Lix) have nonzero Fourier 
components only within the narrow bandwidth k^. Using 
this projection in the definition of the DSF and retaining 
only Fourier components close to q leads to 

S{q,uj) = j rfa;dte'"*-'«^(*t(a;,t)^r(a;,t)*t(o,o)*(0,0)) 

dxdte''^\S{x,t)iljR{xJ)^jj^p,{Q,())d{Q,Q)). (18) 



The next step is to project the interacting system 
Hamiltonian onto the narrow subbands. Because the re- 
duced bandwidth fco is small, the spectrum within each of 
the subbands can be linearized. This makes it convenient 
to bosonize using 



9(x)] 



(19) 



and to write the corresponding terms in the Hamilto- 
nian in the bosonic basis. In Eq. (19), the conventional 

bosonic fields 9 and (j) satisfy a canonical commutation re- 
lation (we use the notations of Giamarchi (2004) though- 
out the text). 



(a;),V6'(a;')] = iTT5{x-x'). 



(20) 



After a projection of the microscopic interactions onto 
subbands and bosonization, the Hamiltonian becomes 
H = Ho + Hd + Hint, where 



^o = g / dx[(yeY + {vm 



(21) 



Hd = j dxS{xMk)-ivdV]d{x), 

Hint = y '^^ [(^fc-feF - yo)PR + {Vk+kp - Vo)pl\ dd'^ . 

The term Hq describes the kinetic energy of the particles 
near the Fermi points. The energy of the impurity d 
is close to ^(fc), its motion is described by H4, and its 
velocity reads 



Vd = 



d^{k) _ k 
'dk 



= — = Vf 
m 



m' 



(22) 



Last but not least, the term Hi„t contains the density- 
density interactions between the impurity and the parti- 
cles near the Fermi points. The densities of right- and 
left-movers are given by 



Pa{x) 



1 
27T 



V{-(j) + ae), 



(23) 



where a = R, L = — . Interactions lead to the forma- 
tion of low-energy particle-hole pairs and are thus crucial 
for the shape of the DSF. Schotte and Schotte (1969) 
showed that Hint can be removed using a unitary trans- 
formation. Indeed, introducing the unitary operator 



U = exp 



dx 



27r' 



[e + mddH (24) 



one finds W{Ho +Hd + Hint)U = Hq + H^ by using the 
phase shifts 



Vd - Vf 



Vq — Vk+kf 
Vd + VF 



(25) 



Physically, the values of S±{k) correspond to the scat- 
tering phase shifts between the deep hole and the low- 
energy particles near the right and left Fermi points in 
the Born approximation. In order to calculate S{q, lu) us- 
ing Eq. (18), the same unitary transformation has to be 
applied to the operators tpR{x) and d{x). The impurity 
operator acquires a phase shift in the rotation. 



UU{x)U : 



d{x). (26) 



Note the similarity to the bosonization formula (19). In- 
deed, the "shake-up" of the particles at the right Fermi 
point caused by the interaction with the deep hole man- 
ifests itself as an additional phase in the bosonic repre- 
sentation of the operator tpii{x), and similarly for the left 
Fermi point. 

After bosonization and rotation, the expectation value 

in Eq. (18) thus factorizes into a term containing the 
bosonic fields <p and 9, and a term containing the impu- 
rity operator d. The dynamics of the fields (p{x) and 6{x) 
is governed by Hq and is therefore linear. Hence, the 
expectation values of exponentials of bosonic operators 
can be calculated straightforwardly (Giamarchi, 2004). 
The impurity dynamics is governed by H^ and leads 



to {d\x,t)d{0,0)) 



-iuj— (q)t 



6{x — Vdt). Upon Fourier 



transformation of the time-dependent correlation func- 
tion, the result for the DSF at its lower edge of support 

reads 



Sjq^uj) 
m/q 



uj-{q) 



for du){q) w — w_(g) > 0. 

(27) 



The threshold exponent of S{q,u!) depends on momen- 
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turn and interaction strength and reads 

2 



Mo(g) = 1-1 + 



TT 



m 



Trig] 



5_ 

(^0 - V,). 



(28) 



For generic repulsive interaction potentials Ho{q) > 
0, so the DSF has a power-law divergence at the 
lower threshold. The expansion of Eq. (27) for 
^o{q)^^{[oJ — oj-{q)]/Suj{q)} <C 1 coincides with the 
leading order logarithmic result (16). The present calcu- 
lation corresponds to a resummation of the leading log- 
arithmic divergences to each order in the perturbation 
series and thus yields the exponent lJ<o{q) to first order in 
the interaction potential Vq. 

For momenta \q\ ~ 2kF, the exponent (28) evaluated 
in the leading order of perturbation theory in the inter- 
action potential coincides with the corresponding limit 
of the linear LL exponents. Indeed, at g — > 2kF, the 
LL theory predicts a power-law divergence, S{q,ui) oc 
9[iO — \v{q — 2kF)W[oJ — \v{q — 2kF)\]^~^, with an expo- 
nent K — 1. For weak interaction, K can be calculated 
perturbatively, and one finds (Giamarchi, 2004) 



K^l-{Vo-V2k^)/{2nvF), 



(29) 



thus reproducing the prediction (28). At g ^ and V{x) 

decaying faster than oc l/x^, {Vq — Vo)/q 0. Then the 
exponent (28) vanishes, and for a fixed q a rectangular 
shape of S{q,oj) is recovered. The latter has a width 
5uj{q) oc q^, height m/q, and is located at the mass shell 
to = VFq- In the limit g — >■ 0, this peak indeed acquires a 
delta-shape as predicted by the LL theory. Later we will 
see that the relation /io (0) = for a short-range potential 
holds beyond the perturbation theory; the rectangular 
shape of S{q, u) at small fixed q is quite generic. 

The procedure that led to the DSF near its lower 
threshold exemplifies a rather versatile framework for 
the perturbative calculation of various dynamic response 
functions near singular thresholds. Generally, singular- 
ities appear whenever conservation laws allow that the 
entire energy of an incoming density- or single-particle 
excitation is transferred to a single particle or hole in 
the system. As illustrated above, the general procedure 
consists of the following three steps: (i) Identification of 
the "deep hole" configuration responsible for the singular 
behavior at the threshold of interest. This configuration 
always follows from momentum and energy conservation, 
(ii) Projection of the Hamiltonian onto a reduced band 
structure containing narrow bands around the deep hole 
and the Fermi points, (iii) Determination of the phase 
shifts due to the interactions between the deep hole and 
particles at the Fermi points by applying a unitary trans- 
formation. 

The shape of the threshold singularity of S{q,uj) for 
w « w+((j') can be obtained similarly. For g > 0, the 



configuration giving rise to this singularity contains a 
hole near the right Fermi point as well as a particle near 
kF + q (see Fig. 1). Projecting the Hamiltonian onto nar- 
row bands around the Fermi points and a narrow band 
around kF + q, and following essentially the same proce- 
dure as before, it was found that for du){q) \oj — u)+{q)\ 
(Pustilnik et al, 2006) 



S{q,ui) 
m/q 



where 



-A»o(9) 



q 

Amvp 



2kf 



V0-V2 
2tivf 



w < uj+{q), 
oj > U!+{q). 

(30) 

(31) 



Hence, power-law singularities with identical exponents 
— /xo(g) appear on both sides of the threshold a;+(g). In 
stark contrast to the noninteracting limit, S{q,io) no 
longer vanishes above the upper threshold. Also note 
that the prefactors are different on both flanks. 

The DSF S{q,uj) remains nonzero even above the up- 
per threshold because any excess energy u — w+(g) can 
be used for the creation of an additional particle-hole 
pair on the left branch. As the excess energy increases, 
the momenta of the two particles and two holes may be 
increasingly far away from the Fermi points. Hence, or- 
dinary second-order perturbation theory works well for 
CO — oj+{q) ^ 5oj. In this range, one finds 



S{q,uj) = 2u{q)- 



VFq 



v^q^' 



(32) 



The DSF for a weakly interacting system is depicted in 
Fig. 1. 

In conclusion, interactions lead to notable changes in 
S{q,u). Instead of the rectangular shape of S{q,oj) 
for any given q in the noninteracting case, interactions 
Iciad to the appearance of power-law singularities at the 
thresholds uj±{q). Moreover, the function no longer van- 
ishes above the upper edge w+(g). 



2. Spectral function 

The procedure employed for the perturbative calcula- 
tion of the DSF can also be used for the calculation of 
the edge singularities of the spectral function (Khodas 
et al., 2007b). In contrast to S{q,oj), the spectral func- 
tion A{k, s) characterizes the response of the system to 
the addition of a single particle or hole. It determines 
the probability for a particle with energy e and momen- 
tum k to enter (or, at e < 0, emerge from) the system 
in a momentum-conserving tunneling event; particle and 
hole sectors correspond to £ > and £ < 0, respectively. 
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FIG. 4 (Color online) Spectral function A{k, e) for interact- 
ing fermions, with notations for the edge exponents. The 
configurations determining the edge exponents fio,- and 
are indicated. In the weakly interacting case, the edge of 
support eth{k) at |fc| < kp coincides with the mass shell 
^(fc), see Eq. (2). In the noninteracting case, the spec- 
tral function A{k,e) is a delta-function at the mass shell, 
Ao{k,e)^S[e-ak)]. 

In the absence of interactions, A{k,£) = d[e — £,{k)], be- 
cause a particle or hole with momentum k can only be 
absorbed if its energy is on the mass shell ^(fc). In an 
interacting system, on the other hand, the spectral func- 
tion will generally become nonzero even away from the 
mass shell ^(fc), because incoming particles or holes may 
give up part of their energy and momentum to excite 
additional particle-hole pairs. 

For \k\ < kp and e < 0, the edge of support of A{k, e) 
coincides with the fermion mass shell. In this case, the 
calculation of the edge singularity is very analogous to 
the previous section. If a particle with momentum k and 
energy e ~ ^(fc) is extracted from the system, it leaves 
behind a hole with momentum near k on mass shell, as 
well as low-energy particle-hole pairs near either of the 
Fermi points. Therefore, it is again sufficient to retain 
narrow bands of widths fco ^ |fc| around k and ±fci?, 
and project the fermion operator as in Eq. (17). The 
remaining calculations closely follow Sec. II.A.I. 

The spectral function in the vicinity of ^(fc) < can 
be calculated as 

cx0[e(fc)-e][e(fc)-e]-^»-. (33) 

The exponent to lowest order in the interaction strength 
reads 

The phase shifts 6± are defined in Eq. (25). For k — > 
kp, the phase shift vanishes linearly for interaction 
potentials decaying faster than in real space. On the 
other hand, (5_ at fc — kp remains finite, = — (y2kf — 



Vq)/ {2vp). Therefore, in the limit k — kp, the exponent 
/io,- coincides with the Luttinger model prediction for 
weak interactions (Luther and Peschel, 1974). 

Note that to lowest order the exponent is quadratic in 
the interaction potential Vk, which is in contrast to the 
result (28) for the DSF. Moreover, the interactions of the 
impurity with left-movers and right-movers are equally 
important. 

Outside the regime e < and |fc| < kp, the edge of 
support of the spectral function no longer coincides with 
the mass shell, but it is still determined by kinematic 
considerations: the injection of a particle or hole and the 
ensuing creation of particle-hole pairs due to the interac- 
tions must respect momentum and energy conservation. 
For weak interactions, the edge of support can be de- 
termined quantitatively by using the Lehmann spectral 
representation (Abrikosov et at, 1963). Let us focus on 
the particle sector (e > 0) and the momentum range 
kp < k < 3kp, where 

A{k,e) = J2\{ml\0)\Hk^p^.,,5{e ~ Ef). (35) 
l/> 

The initial state |0) corresponds to the ground state of 
the system and the sum runs over a complete basis {|/)} 
of the Fock space. The energies and momenta of the final 
states are denoted hy Ej and Pf, respectively. In the 
absence of interactions, the ground state |0) is the filled 
Fermi sea and the only final state with nonzero overlap 
is I/) = *j.|0). This immediately leads to A{k,e) 

In an interacting system, on the other hand, the ground 
state |0), written in the basis generated by the operators 
^I'j. and may contain particle- hole pairs. Therefore, 
a nonzero overlap can also be achieved for final states 
which contain additional excitations. The simplest set of 
such final states is 

l/)=«*fc3lO), (36) 

and it can be parametrized by the momenta fci,fc2,fc3. 
The edge of support of A{k,e) can be determined by 
enforcing momentum and energy conservation and find- 
ing the configuration with lowest excitation energy Ef. 
Within perturbation theory, the energy Ef is determined 
using the noninteracting Hamiltonian. The result for the 
particle sector (e > 0) at /c > fc^? is fci = ^2 = fc^ and 
fca = 2kp — k, i.e., at the edge of support the entire energy 
is carried by a single hole with momentum 2kp — k. The 
energy of this configuration and thus the edge of support 
of A{k, e) is given by 

-a2kp -k)= vp{k - kp) - ~ . (37) 

The configuration yielding the threshold (37) remains the 
lowest energy state with total momentum kp < k < 3kp 
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even if compared to states with a higher number of 
particle-hole pairs than Eq. (36) . Hence, the edge of sup- 
port in this region coincides with the shifted and inverted 
mass shell —£^{2kp — k). 

The threshold configuration consists of a hole at mo- 
mentum 2kF — k and two particles with momenta in- 
finitesimally close to kp. Hence, the Hamiltonian must 
be projected onto narrow bands (of widths ko <^ k — kp) 
around these momenta. However, the initial particle with 
momentum k is outside this band structure, so the oper- 
ator \E'|. would vanish in a naive projection of the form 
(17). The solution is to use a Schrieffer- Wolff transfor- 
mation (Schrieffer and Wolff, 1966) to derive a projection 
which is of first order in the interactions. Then, it can be 
shown explicitly that the particle at momentum k creates 
the threshold configuration (Khodas et al., 2007b), 

*l 0^ $Z 'pRM'pRM^k.+k,, (38) 

fci ,k2 

where ki,k2 «C ko. The impurity operator d creates a 
hole near momentum 2kp — k, whereas the operators 
fei 2 create particles close to the right Fermi point. 
The fact that a single incoming particle now has to form 
three excitations in order to be able to tunnel into the 
system opens up the phase space available for the pro- 
cess, due to the free variables ki and k2 in the projection 
(38). The spectral function in the vicinity of the edge 
is proportional to the Fourier transform of (\I'fc(i)\E'^(0)). 
Neglecting the interactions between the subbands, a di- 
rect calculation of this correlator using Eq. (38) leads to 
A{k,e) cx [e + S,{2kp — fc)]'^, i.e., a convergent threshold 
behavior. An increased number of excitations in the pro- 
jection of the physical fermion operators, as required at 
momenta > kp, generally leads to more convergent 
power-laws. It is a recurring feature that all thresholds 
are characterized by configurations in which the entire 
energy is carried by a single particle or hole, while addi- 
tional particle-hole pairs reside close to the Fermi points. 

Interactions between the impurity d and the particles 
near the Fermi points again lead to a correction to this 
exponent. Using a mobile impurity Hamiltonian, the en- 
suing calculation is analogous to the previous discussion. 
The operator ^'^ can be bosonized as oc e^^^'^~^''d, 
and the spectral function near the edge behaves as 

A{k, e) oc e[e + ^{^kp - k)] [e + ^{2kp - k)]-''''+ , (39) 

where 



1^1,+ = 




The thresholds in the other sectors of the (fc, £)-plane 
can be derived similarly. The support of the spectral 
function in the weakly interacting limit and the lowest- 
energy configurations at the respective edges are dis- 
played in Fig. 4. 



The configurations which give rise to singularities at 
the edges of support are always stable: they represent 
the excitations of lowest energy for a given momentum 
and are thus protected from decay by conservation laws. 
On the other hand and in striking contrast to the non- 
interacting case, the mass shell ^(k) no longer forms the 
edge of support for momenta k > kp. Instead, it now lies 
within a continuum of excitations. Therefore, particles 
on the mass shell are generally subject to decay. In the 
spectral function, this will give rise to a broadening of 
the singularity at £ = ^(fc). This will be discussed more 
in detail in Sec. IV. A. 



B. The universal limit of nonlinear Luttinger liquids 

The analysis presented in the previous section predicts 
the dynamic response functions for weakly interacting 
Fermi systems at arbitrary momenta. In many realistic 
systems, however, the interaction energy can be of the 
same order as the kinetic energy, thus making perturba- 
tion theory inapplicable. To treat such systems, having a 
theory which accounts for the interactions exactly would 
be desirable. If the fermionic spectrum is strictly linear, 
all dynamic response functions at low energies can be cal- 
culated exactly (Dzyaloshinskii and Larkin, 1974; Luther 
and Peschel, 1974). Close to Fermi points, it is tempting 
to consider the band curvature as a small perturbation to 
the linear spectrum. The resulting corrections to single- 
variable correlation functions (for example in the fermion 
distribution function rife = (^E'l^'fe)) are indeed imiformly 
small. This is not the case, however, for the dynamic re- 
sponse functions. We will see here that the true values of 
the threshold exponents are different from the predictions 
of the linear LL theory even in the limit | A; | — > kp (Imam- 
bekov and Glazman, 2009b). The frequency domain near 
the threshold where these strong deviations take place, 
narrows down as (|A;| — kp)^. 

Phenomenological bosonization is the obvious ap- 
proach to tackle a strongly interacting ID system in 
the low-energy regime (Efetov and Larkin, 1975; Hal- 
dane, 1981a,b). This approach is based on rephrasing 
the fermionic problem in a bosonic language using the 
bosonization identities (19) and (23). Using this basis 
offers the advantage that a density- density interaction 
between the physical fermions produces a quadratic term 
in the bosonic variables. At low energies, where only de- 
grees of freedom close to the Fermi points are involved, 
bosonization allows an exact treatment of the interaction. 

However, for the quadratic spectrum (2) the kinetic 
energy becomes more complicated when expressed using 
the bosonic fields 0{x) and (j){x). The kinetic energy den- 
sity of an ideal gas in the ground state can be calculated 
by integrating the spectrum ^(fc) over k G [—kp,kp]. 
Local fluctuations of the left and right particle densities 
shift the Fermi points, kp'^{x) = ±[kp + wpii,L{x)], and 
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therefore change this energy density. By expressing the 
density fluctuations using Eq. (23), the kinetic Hamil- 
tonian can be derived. Near the two Fermi points, the 
spectrum for noninteracting fermions can be expanded 
as ^(fc) « vpi^k — kp) + {k^ kp)^ / {2m). Its linear com- 
ponent together with the interaction term produces the 
conventional LL Hamiltonian, 



LL 



V 

2^ 



dx 



(41) 



where v denotes the rcnormalized Fermi velocity, and K 
is the Luttinger parameter, which is in the interval < 
K <1 for repulsive interactions. For the noninteracting 
system, K ~ 1 and v = vp. The quadratic component of 
the spectrum, on the other hand, leads to cubic terms in 
the bosonic fields (Haldane, 1981b), 

H^i = -^j dx [{Wcj>f + 3(V(/.)(V0)2] . (42) 

As soon as the cubic band curvature terms are taken 
into account, an exact diagonalization of the Hamilto- 
nian ifLL + Hnl in terms of the fields (j) and 6 is no longer 
possible. The most obvious route is to treat H^u which 
is proportional to 1/m, as a perturbation. Diagrammat- 
ically, the terms in Eq. (42) correspond to three-boson 
interaction vertices. It turns out, however, that such 
an endeavor is far from trivial bccaiisc the bosonic self- 
energy diverges at the mass shell to = vk (Samokhin, 
1998). The physical reason is the linear spectrum of the 
bosonic modes: all bosonic excitations propagate with 
the same velocity v independent of their momentum and 
therefore - semiclassically speaking - have an infinite time 
to interact. 

This difficulty precludes a straightforward calculation 
of the DSF near the mass shell. Far away from the mass 
shell {lo ^ vq), on the other hand, the perturbation the- 
ory in the band curvature is convergent and a high-energy 
tail in S{q,Lo) emerges in the order (1/m)^ (Pereira et al, 
2007). For small nonzero interactions, this agrees with 
the perturbative result (32). 

An expansion of the free-fermion result (10) in orders 
of 1/m reveals that each individual term of the series 
diverges at the mass shell. Therefore, in order to ac- 
cess S{k, cj) close to the mass shell, an efficient resum- 
mation scheme is required. Standard procedures like the 
Born approximation fail because they still produce a di- 
vergence at w = vq. Various approximate schemes have 
been developed (Aristov, 2007; Pirooznia and Kopictz, 
2007; Pirooznia et al, 2008; Schonhammer, 2007; Teber, 
2007), but even the free-fermion result has been repro- 
duced in the bosonic basis only up to the order (1/m)^ 
(Pereira et al, 2007). 

Many of the complications which plague the bosonic 
perturbation theory can be avoided by using a ba- 
sis of fermionic quasiparticles (Mattis and Lieb, 1965; 



Rozhkov, 2005). Similar to the bosonic fields, the quasi- 
particles remain free in the case of a strictly linear spec- 
trum whereas a nonzero band curvature of the under- 
lying particles, together with the inter-particle interac- 
tions, leads to interactions between the fermionic quasi- 
particles. In contrast to the bosonic theory, however, the 
spectrum nonlinearity of the physical fermions also en- 
tails a band curvature of the quasiparticles. Hence, quasi- 
particles with different momenta propagate at different 
velocities. We will see that the scattering between them 
at momenta close to the Fermi points can be treated, for 
instance, within the Born approximation. Moreover, the 
use of a fermionic basis allows again the introduction of 
a mobile impurity Hamiltonian and thus connects to the 
method employed in the perturbative calculation. 

To illustrate the origin of the fermionic quasiparticle 
representation of the Luttinger model, let us diagonalize 
the Hamiltonian (41) by introducing the rescaled fields 



(43) 



e{x) = VK0{x), 4>{x) = 0(a;)/Vif 



Since this is a canonical Bogoliubov transformation, the 
fields 6{x) and 4'{x) are still canonically conjugate. In the 
new variables and (p, the Hamiltonian is indistinguish- 
able from the bosonized version of the Hamiltonian of 
free fermions with linear spectrum. These free left- and 
right-moving fermionic quasiparticles can be defined by 
using the bosonization identity (Haldane, 1981b; Luther 
and Peschel, 1974; Mattis, 1974) on the rescaled fields. 



^a{x) oc eyi\>{—i[a(f){x) — 9{x)\}, 



(44) 



for a = R.L 



Here * 



are creation 



and annihilation operators for quasiparticles on the right 
(left) branch, satisfying usual fermionic commutation re- 
lations [as Tisually, we didn't write out the Klein factors 
explicitly (Giamarchi, 2004)]. In terms of quasiparticles, 
ifLL becomes 



LL 



-IV 



dx 



: ¥j^{x)V^r{x) : - : ^1{x)V^l{x) 



(45) 



Colons indicate the normal ordering with respect to filled 
Fermi seas: for the right (Icift) branch all states with 
negative (positive) momenta are occnipied. The relations 
between pu(L) and Pr(^l) are linear, and follow from Eqs. 
(23) and (43). 

To proceed further, we need relations between ^r^l) 
and ^R(L)- Since both ^r and "^r carry momentum +kp 
and change the total number of particles by one, '^r 
should contain '^r and low-energy particle-hole pairs. 
Using the bosonization formula for expressing in 
terms of (j) and 6, and using Eq. (44) to "pull out" an 
operator ^r, one finds for right-movers 



^r{x) = Fr(x)'^r(x). 



(46) 
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An analogous expression holds for the left-movers. The 
string operator Ffi{x) is an exponential of the left-moving 
and right-moving quasiparticle densities (23), 

Fr{x) = exp |-i j dy[d+pR{y) + (5_pl(j/)]| . (47) 

The described refermionization procedure defines 
uniquely the two parameters 



2tt 



= 1- 



1 



K 



2VK 



<0, 



5_ 
27r 



1 



2VK 



K 
2 



(48) 



According to Eq. (46), the annihilation of a physical 

right-moving fcrmion in the liquid causes the annihilation 
of a right-moving quasiparticle. In addition it leads to a 
shake-up of the Fermi seas of the left-moving and right- 
moving quasiparticles which is described by Fr. The pa- 
rameters 6± can be interpreted as the phase shifts for the 
scattering of quasiparticles at ±kp off the right-moving 
hole in the quasiparticle distribution. 

The mapping between the interacting physical 
fermions and the noninteracting fermionic quasiparticles 
can also be performed directly via a unitary transfor- 
mation (Mattis and Lieb, 1965; Rozhkov, 2005), without 
invoking bosonization as an intermediate stop, and pro- 
duces identical results. The Hamiltonian (45) together 
with Eqs. (47) and (48) reproduces the usual results for 
the fermionic Green's function (Rozhkov, 2005). 

The quadratic spectrum (2) of the physical fermions 
leads to additional terms in the quasiparticle Hamilto- 
nian. Most importantly, a quadratic term emerges in the 
quasiparticle spectrum (Rozhkov, 2006, 2008, 2009). It 
is reflected in the Hamiltonian 



- 2^ / 



dx 



: (V*ijV*ii) : + : (V^LV^-l) : 



(49) 

The effective mass m depends on interactions, and us- 
ing the methods described in Sec. II. C, it is possible to 
express it via low energy properties as (Pereira et al, 
2006) 



1 V d r— 



m 



(50) 



where n is the chemical potential. From Eqs. (45) and 
(49), the quasiparticle mass shell for fc « fci? is given by 



i{k) = v{k - kp) + 



(fc - kpf 
2m 



(51) 



In addition, a spectrum curvature of and interactions 
between the physical fermions generally lead to inter- 
actions between the fermionic quasiparticles. One such 
term describes an interaction between quasiparticles on 
opposite branches (Rozhkov, 2006), 



H'int = ig 



E 

a=R,L • 



dxp-a 



: §UV*a) : - : (V§t^)§„ : 

(52) 



In order to understand the effect of this term, we consider 
the scattering phase shift between two quasiparticles with 
momenta p + kp {\p\ 4; kp) and —kp as in Sec. II. A. 
Fourier transforming Eq. (52) reveals that for small p, 
the interaction potential is proportional to p. This has 
to be compared with the difference in velocities, which 
according to Eq. (51) is finite and close to 2vp. Therefore, 
the term (52) produces only small additional phase shifts 
of order p/kp <C 1. 

Another interaction term can appear as a consequence 
of a momentum dependence of the physical interaction 
potential. It corresponds to a density-density interaction 
between quasiparticles on the same branch (Imambekov 
and Glazman, 2009b), 



TT" 

-"■int 



dpW{p)paip)pc.{-p). (53) 



The Luttinger parameter K already accounts for inter- 
action processes with momentum exchange close to zero, 
so the additional quasiparticle interaction potential must 
fulfill W{0) = 0. Moreover, it should be symmetric 
W{p) = W{—p). For small p, the scattering phase shift 
between two particles with momenta p + kp and kp is 
now of order rhW{p)/p. This phase shift therefore be- 
comes small for generic interaction potentials which fulfill 
W{p) (X p^ for p — ^ 0. This requirement is only violated 
for real-space potentials decaying 1/.T^ or slower. 

In particular, this excludes Haldane-Shastry (Haldane, 
1988; Shastry, 1988) or Calogero-Sutherland (Calogero, 
1969, 1971; Sutherland, 1971, 2004) type models. 

The above arguments establish that Eqs. (45) and 
(49), written in terms of fermionic quasiparticles, serve 
as the universal low-energy Hamiltonian which captures 
the leading role of the spectrum nonlinearity. The crucial 
advantage of the quasiparticle representation is that un- 
like in bosonic extensions of the LL theory, the universal 
Hamiltonian does not contain interactions. The nonlin- 
earity of the spectrum lifts the "accidental" degeneracies 
existing in the Luttinger model, and allows to treat small 
interaction terms (such as given in Eqs. (52) and (53)) 
in addition to the universal Hamiltonian within conver- 
gent perturbation theory. Surprisingly, the only addi- 
tional phenomenological parameter is the effective mass 
m, which sets the energy scale for nonlinear effects. Be- 
low we will work out universal predictions for the low- 
energy DSF and the spectral function beyond the linear 
spectrum approximation. We will show that the true val- 
ues of the threshold exponents differ from predictions of 
the linear LL theory, but nevertheless can be expressed 
as functions of the LL parameter K only. 

The relation between the total physical particle density 
and the total quasiparticle density is linear, Pl + Pr = 
VK[pl + Pr]- Therefore, the shape of the DSF S{q,Lo) 
for q <^ kp remains identical to the noninteracting result 
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(10), albeit with a renormalized mass m, 
Km . ^ 



S{q,uj) 



q 

2m 



vq\ 



(54) 



This agrees with the perturbative resuhs (27) and (30) 
in the hmit g — >■ 0. Power-law singularities at the upper 
and lower thresholds Lo±{q) — vqztq'^ /{2rh), respectively, 
as well as a high-energy tail would be produced by the 
quasiparticle interaction term (52) which becomes rele- 
vant only beyond the limit q/kp <C 1. 

Let us first investigate the spectral function A{k, e) 
in the presence of band curvature in the hole sector 
{0 < kp — k <^ kp and e < 0). In this region, the kine- 
matic edge of support of A{k,e) coincides with the mass 
shell (51), so we assume e « £,{k). According to Eq. (46), 
the extraction of a physical particle with momentum k 
from the system leads to the formation of a quasiparticle 
hole with momentum near k on mass shell, and the ex- 
cess energy is used for the formation of excitations near 
the Fermi points. In analogy to the previous section, the 
Hamiltonian (45)- (49) as well as the single-particle oper- 
ator (46) should be projected onto small bands around 
the momenta ±kp and k. The right-moving quasiparticle 
operator is projected as 



(55) 



and — > iPl{x) is used for left-movers; here, iPr^l{x) 
denote quasiparticles within a narrow momentum range 
around the Fermi points ^, while d{x) denotes an impu- 
rity with momentum close to k. Projecting the string 
operators Fr{x) requires a projection of the quasiparti- 
cle densities Pa{x) = ^^^{x)^ a{x) in Eq. (47). This pro- 
duces terms containing the quasiparticle densities in the 
narrow bands around the Fermi points, ipl^{x)4'a{x). In 
addition, the projection leads to impurity terms of the 
form d^{x)d{x) and mixed terms ^pl^{x)d{x). However, 
the former is a constant, and the latter can be neglected 
close to the edges because they require a higher energy. 
Therefore, projecting the string operators corresponds to 
replacing pa{x) in Eq. (47) with ipl^{x)'ipa{x) ■ 

A projection of the kinetic Hamiltonian i/LL + H'^^^ 
onto the three subbands leads to the mobile impurity 
Hamiltonian Hq + where 



Hf) = —iv I dx 



^I^I{x)V^r{x) : - : ^bl{x)VMx) 



Hd= I dxd){x) i{k)-ivd^ d(x), 



(56) 



^ Hereinafter, we reserve the upper-case symbols l,r{^) a-nd the 
lower-case symbols iPl,r{^) for operators before and after the 
projection, respectively 




FIG. 5 (Color online) Band structure for right-movers used 
for the calculation of A{k, e) for < kp — k <^ kp and e < 0, 
see Eq. (55). It contains the impurity band near momentum k 
and a low-energy band near the right Fermi point +kp. The 
injection of a hole with momentum k and energy e into the 
system leads to the formation of a hole (empty circle) with 
momentum near k on the mass shell, as well as a particle-hole 
pair near the Fermi point. For |e— ^(fc)| ^ (k—kp)^ /{2m), the 
subbands are well separated and the mobile-impurity Hamil- 
tonian can be applied. 



where the impurity velocity is Vd = d^{k)/dk. Note 
that the requirement that the impurity band be sepa- 
rated from the low-energy bands means that this mo- 
bile impurity Hamiltonian can be used to calculate the 
spectral function at energies e satisfying \e — £,{k)\ <C 
{k — kpY / {2'm). This is illustrated in Fig. 5. 

Then, the spectral function in the region < kp — k<^ 
kp and £ < becomes 

A{k,e)= I dtdxe-^^\d\x,t)d{Q,Q)) 



X {Fl{x,t)FR{Q,Q)) 



(57) 



The impurity correlation function can easily be calcu- 
lated using the noninteracting Hamiltonian Hd- The 
free string correlation function can most conveniently 
be derived by bosonizing and Fr. One finds that 
A{k, e) cx 6i(|(fc) -e)[e- £,{k)]-^'°^- , where 



1 - 



(58) 



As in the previous section, one of the consequences 
of a nonlinear spectrum with m > is that the edge 
of support of the spectral function in the particle sec- 
tor no longer coincides with the quasiparticle mass shell 
^(A:) but rather with the shifted and inverted quasipar- 
ticle spectrum. For instance, for < k — kp <^ kp, 
the threshold is given by —£,{2kp — k). The threshold 
configuration generated at this edge contains a quasi- 
hole with momentum 2kp — k < kp, as well as two 
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quasiparticles with total momentum 2k p- The power- 
law at this threshold can be derived using the same 
mobile-impurity Hamiltonian (56) but the projection of 
the quasiparticle operator has to be done analogous to 
Eq. (38). Due to the creation of additional particles, the 
spectral function at this edge is convergent. One finds 
A{k,£) (X e{s + i{2kF - k))[e + i{2kF - k)]-^'^'+, where 



Ml, 



1 



27r / 



2tt 



(59) 



Similar to the weakly interacting case, the singularity at 
the mass shell for k > kp and e « ^(fc) now lies within a 
continuum of excitations. Since interactions neglected in 
the universal Hamiltonian generally lead to a nonzero de- 
cay rate for quasiparticles on mass shell, one may expect 
that the singularity may be smeared. Within pertur- 
bation theory, the decay rate close to the Fermi points 
scales as r oc (k — kp)^ (Khodas et al., 2007b), and 
will be discussed in detail in Sec. IV. A. However, the 
power-law behavior is expected to be observable within 
an energy window of width (fc — kp)^/{2rn) around the 
mass shell (Imambekov and Glazman, 2009b). There- 
fore, the singularity is indeed resolved for /c — kp. Using 
the mobile-impurity Hamiltonian, one finds power laws 
with identical exponents on both sides of the mass shell, 
A{k, s) oc |e — £,{k)\~'^°-- , where the exponent is given by 
Eq. (58). Despite the fact that the exponents on both 
sides of the singularity are identical, the prefactors are 
not. It can be shown that (Imambekov and Glazman, 
2009b) 



A[k,i{k)+Se] 
se^o A[k,i{k)-Ss] 



sm 



lim 



sm 



1-K j 



(60) 



The exponents at the mass shell are identical in the par- 
ticle and hole sectors, and they differ from the predic- 
tions of the linear LL theory. As has been established in 
Sec. II. A. 2, the expansion of the exponent (58) to lead- 
ing order in X — 1 coincides with the Luttinger model 
prediction (62). However, a difference to the LL expo- 
nent emerges in the order [K — 1)^, and is of order one 
for strong interactions (|ii' — 1| ^ 1). 

The result of the linear LL theory is recovered for en- 
ergies further above the mass shell, e — ^(/c) 3> (fc — 
k pY / {2fii) . As illustrated in Fig. 5, in this regime the 
momentum bands encountered in the projection start to 
overlap, and the mobile-impurity Hamiltonian ceases to 
be applicable. From the point of view of an incoming 
particle, the band curvature becomes irrelevant at these 
energy scales. Therefore, the spectral function can be 
calculated using 



A{k,£) c<. j dtdxe 



let^—ikx 



X (^Jj(a;, t)Fl{x, t)FniO, 0)^0, «)), (61) 



and by assuming that the quasiparticle spectrum is lin- 
ear. In this case, the time evolution of the quasiparticle 
densities becomes simple, pa{x,t) = pa{x — avt) and al- 
lows for a calculation of the dynamics of the string opera- 
tor Fji{x,t). This makes it possible to calculate Eq. (61) 
in the basis of fermionic quasiparticles (Rozhkov, 2005). 
Alternatively, Eq. (61) with a linear spectrum may be 
calculated by bosonizing it. As a result, for e — ^{k) ^ 
{k - kpf/{2fh) one finds A{k,£) oc [e - |(fc)]-'*LL with 
the LL exponent 



Mll = 1 



2-.) 



(62) 



Compared to the exponent /ig,- near the mass shell, see 
Eq. (58), the 5+ term is missing. The phase 5+ is the 
scattering phase shift for interactions between the impu- 
rity at momentum k ^ kp and particles near the right 
Fermi point. Because the impurity band around k and 
the band around +kp start to merge for energies fur- 
ther away from the mass shell, this phase is absent in 
Eq. (62). The phase in contrast, is brought about by 
interactions with particles near the left Fermi point, and 
continues to be present. 

Understanding the crossover between regions with ex- 
ponents (59), (58) and (62) requires a calculation of 
A{k, s) at intermediate energies. Due to kinematic con- 
straints, only the nonlinearity of the spectrum of right- 
movers is important for k > kp and 6 > ^{k). There- 
fore, an analysis of the exact dynamics of the string 
operators Fji{x) is needed. Similar correlators have at- 
tracted attention recently in connection with the non- 
linear quantum shock wave dynamics of a free Fermi 
gas (Bettelheim et al, 2006a,b; Bettelheim et al, 2007; 
Bettelheim et al., 2008). Even though the quasiparti- 
cles are nonintcracting, this is a highly nontrivial prob- 
lem due to the nonlinear spectrum ^(fc). However, 
since it is essentially a single-particle problem, it can be 
mapped onto a calculation of certain infinite-size deter- 
minants (Abanin and Levitov, 2004, 2005) and then tack- 
led numerically (Imambekov and Glazman, 2009b) . The 
spectral function for momenta fc « fcj? at arbitrary ener- 
gies £ can be written as a function of a single variable, 



A{k,e) <x A{x), where a; = 



e — v(k — kp) 



(k- 



kpY/{2m) 



(63) 



where the function A{x) depends only on the Luttinger 
parameter K. In this sense, for momenta close to the 
Fermi points, the universality present in the linear LL 
theory persists even in the presence of band curvature. 
The results of a numerical evaluation of A{x) for a par- 
ticular value of K are shown in Fig. 6. 

Note that the width of the energy window where the 
exponent differs from the predictions of the linear LL the- 
ory is proportional to {k — kp)"^. The fact that this width 
vanishes quadratically for fc — >■ fcj? is consistent with the 
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FIG. 6 The numerically evaluated function A{x), see Eq. (63), 
for a fixed value of the Luttinger parameter K = 4.54 (corre- 
sponding to 5+/{2'k) = —0.3). This function determines the 

behavior of A{k, e) in the crossover region between the thresh- 
old (left end) and the linear Luttinger liquid (LL) asymptotes. 



irrelevance of the band curvature in the universal limit 
k kp- 



C. Phenomenology beyond the low-energy limit: the 
mobile quantum impurity model 

In Sec. II. A, the dynamic response functions were cal- 
culated perturbatively for weak interactions and at ar- 
bitrary momenta. A complementary result, valid for ar- 
bitrary interactions but only close to the Fermi points, 
was derived in Sec. 11. B. In this section, the basic ideas 
of the two approaches, refermionization and the use of a 
quantiim impurity model, will be combined in order to 
obtain phenomenological relations determining the expo- 
nents of the threshold singularities at arbitrary interac- 
tion strengths and momenta. 

To be specific, let us discuss the edge of support of 
the spectral function in the momentum region < kp 
in the hole sector, i.e., for energies e < 0. For a given 
k, the creation of a hole will generally require a nonzero 
minimum energy. At zero temperature, this entails the 
existence of a sharp edge of support which we will denote 
by £th{k), see Fig. 4 for an illustration. The extraction 
of a physical particle is impossible for |e| < lethC^)!- For 
a generic system with arbitrary interaction potential, the 
exact shape of eth(^) is not known, because it is related 
to the exact eigenenergies which crucially depend on the 
microscopic nature of the interactions. 

For noninteracting systems at zero temperature, the 
Fermi momentum kp is defined as the momentum of the 
highest occiipied single-particle state in the Fermi sea. 
The concept of a Fermi momentum can be extended to 
gapless interacting systems by defining kp as the small- 
est positive momentum k for which A{k,e = 0) 7^ 0, i.e., 
at momentum ±kp, the system is capable of absorbing 
even infinitesimal quanta of energy. Applied to this defi- 
nition, the Luttinger theorem (Blagoev and Bedell, 1997; 
Luttinger, 1960; Yamanaka et al., 1997) ensures that the 



value of kp is independent of the interaction strength. 
Therefore, the edge of support still satisfies ethi^kp) = 
even in an interacting system. 

In Sec. 11. A and Sec. II.B, it was established for weak 
interactions and for momenta close to Fermi points for 
arbitrary interactions, that the behavior of the spectral 
fimction near the edges of support can be understood 
by introducing effective models of impurities moving in 
LLs. It is important to emphasize that, similar to Fermi 
liquid quasiparticles, such impurities have definite overlap 
with the original fermionic operators, and possess the 
same quantum numbers. Extracting a physical fermion 
with momentum k and energy close to £th(fc) creates an 
impurity with momentum near k which carries almost the 
entire energy, as well as additional low-energy excitations 
near Fermi points. It is not the overlap between the 
impurity and the fermion, but the "shake-up" of low- 
energy modes near Fermi points which causes the power- 
law divergences in the spectral function as a function of 
the excess energy eth(fc) — e > 0. 

Combining these two limits, we may expect that even 
away from the low-energy limit for arbitrary interactions, 
the threshold properties of response functions near the 
edges of support can be described in terms of mobile- 
impurity Hamiltonians (Pereira et al., 2008). For this 
purpose, the interacting Hamiltonian is projected onto 
a band structure which contains two subbands around 
the Fermi points -iikp and one mobile-impurity subband 
around the momentum k. By continuity from the limit 
of the universal nonlinear LL, this impurity must carry 
the same quantum numbers as the fermion, and its over- 
lap with the fermion must be finite. The width of all 
subbands is small compared to the distance between k 
and the Fermi points, so the spectra in all subbands can 
be linearized. The low-energy excitations near -iikp can 
be modeled using a linear LL Hamiltonian with Fermi 
velocity ?; and Luttinger parameter K . In analogy to the 
previous section, the mobile impurity Hamiltonian thus 
becomes for < kp, 



H, 



dx 



K{V0f + ^{Vcl>f 



Hd = J dx d^x) [sth{k) - ivdV] d{x). (64) 

It is assumed throughout the paper that microscopic in- 
teractions decay faster than oc 1/x, so that compressibil- 
ity and sound velocity are finite. The impurity velocity 
is given by Vd = dsth{k)/dk. In general, there are inter- 
actions between the impurity and right- and left-moving 
particles near the Fermi points, 

Hint = I dx [VR{k)pR{x) + VL{,k)pL{x)] d{x)d\x) 



dx 



VR(/e)V- 



27r 



Vl(A;)V- 



27r 



rfrf^ (65) 
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Unlike in the perturbative regime, Eq. (21), where Vii{k) 
and VL{k) can be directly evaluated, here they have to be 
treated as momentum-dependent phenomenological pa- 
rameters. Similar to Sec. II. A, they determine the expo- 
nents of threshold singularities in all dynamic response 
functions. 

The Hamiltonian (64)- (65) was derived from simple 

phenomenological considerations, and provides a gener- 
alization of the LL theory beyond the low-energy limit. 
We note that effective quantum impurity Hamiltonians 
have been used before in the literature to describe the 
response properties when the effective impurity belongs 
to a band different from the one which forms the LL (Ba- 
lents, 2000; Carmelo et at, 1999; Castella and Zotos, 
1993; Friedrich et al., 2007; Furusaki and Zhang, 1999; 
Lamacraft, 2009; Mishchenko and Starykh, 2011; Sorella 
and Parola, 1996, 1998). Here we extend their applica- 
tion to the case when both bands coincide. Similarly to 
the case of a conventional LL, the Hamiltonian (64)-(65) 
predicts not only the edge exponents in terms of VR,{k) 
and V^ik), but also the structure of the finite-size theory 
(see Sec. II. G). For integrable models, some finite-size 
properties can be analytically extracted from the exact 
solutions (Cheianov and Pustilnik, 2008; Imambekov and 
Glazman, 2008; Karimi and Affleck, 2011; Pereira et ai, 
2008, 2009; Shashi et al, 2010; Zvonarev et ai, 2009b), 
and in combination with Eqs. (64)- (65) lead to nonper- 
turbative predictions for the edge exponents. Various 
examples of this procedure will be reviewed in Sec. III. 
In the rest of this section, we will describe an alternative 
procedure to phenomenologically relate the edge expo- 
nents to eth{k) only, which works for generic Galilean 
invariant systems (Imambekov and Glazman, 2009a). 

Basically, the parameters Vr(A:) and VL(fc) can be fixed 
by considering the shift of the edge eth(k) as a reaction 
to a uniform density variation and to a Galilean boost. 
These variations change and p^, and the resulting 
change in eth{k) is determined by Vn{k) and Viik). First, 
consider the effect of a uniform variation of the system 
density 5p. Because the edge of support is measured with 
respect to the chemical potential p,, its position shifts by 



in {Vcf)) . The total shift in the energy of the state reads 



SE 



dp dp 



Sp- 



dp 



K 



Sp. (66) 



The second equality follows from the general result for 

the compressibility of a LL (Giamarchi, 2004). Next, 
SE is calculated using the mobile impurity Hamiltonian 
(64)-(65), where a uniform density variation leads to a 
nonzero expectation value (V0) = —nSp. We use this 
Hamiltonian to calculate the energy of a state containing 
an impurity at momentum k, and compare the energies 
with and without the density variation Sp. A difference 
only emerges in (Hint) because it contains a term linear 



SE = -Ym±^sp. 



(67) 



Equating the energy shifts (66) and (67) leads to Vr(/c)-|- 
Vt(fc) = -2[9eth(fc)/9p + nv/K]. 

A second equation is needed to fix the difference 
Vnik) — VL{k). For Galilean invariant systems, it can 
be derived by considering the shift of the edge position 
as a response to a boost with a small velocity 6u (Baym 
and Ebner, 1967). Let us consider a hole at the edge 
of support in the moving frame which has momentum 
k' and energy e' = eth(fc'). Galilean invariance requires 
that in the rest frame it has momentum and energy 

k = k' + mSu, e = e' + k'Su + mSu'^/2. (68) 

Since the hole was at the edge in the moving frame and 
the system is Galilean invariant, this also corresponds to 
the edge in the rest frame. Thus the shift of the threshold 
in the rest frame to linear order in Su equals 

SE' = [eth(^ — mSu) + {k ~ mSu)Su + mSu^/2 — £th(fc)] 
deth{k) 



k — m- 



dk 



Su. 



(69) 



Again, we can now calculate the same energy using 
the mobile impurity Hamiltonian (64)- (65). The boost 
changes the Fermi momenta of the right-moving and left- 
moving particles, kp'^ = zLkp + mSu, and thus leads 
to a difference between the densities of right- and left- 
movers. This, in turn, leads to a nonzero expectation 
value (V0) mSu. Substituting this into Eq. (65), we 
calculate the energy of a state with impurity at momen- 
tum k with the help of Eqs. (64) and (65). Subtracting 
the corresponding energy for Su = 0, we find the energy 
shift due to the boost. 



SE' = -Ym^mSu. 
27r 



(70) 



Combining Eqs. (69)-(70) leads to Vnik) - ^^(fc) = 
2TT[detii{k) / dk — k/m\. Thus, both interaction potentials 
^R.Lik) can be expressed in terms of derivatives of the 
threshold ctii(fc) for arbitrary momenta —kp < k < kp. 

In order to calculate; the dynamic correlation functions, 
the interaction Hamiltonian Hint, see Eq. (65), is re- 
moved like in Sec. II. A using a unitary transformation 
U^{Ho + Hd + H,nt)U, where 



U = exp 



(71) 



and 6 and (p are defined in Eq. (43). Such a transfor- 
mation removes Hint if the momentum dependent phase 
shifts S+ (k) and S- (k) are chosen such that 



S_{k) 



Vl 



S+{k) 



Vr 

Vd-v' 



(72) 
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where ^^(fi) arc the couphngs which are obtained after 
rescahng the bosonic fields using Eq. (43). They are re- 
lated to parameters of Hint by 



(73) 
(74) 



{Vl-Vr)/VK = Vl-Vr, 
(Vl + Vb,)^ =Vl + Vr. 



Using the expressions for Vii{k) ± Vl(/c), they can now 
be expressed in terms of derivatives of eth{k), 



S±{k) 
27r 



2|--[- ^gth(fc) 



± VK 



dk ^ 
Idsthik) 

TT 



k_ _ dethjk) 
m dk 



dp 



V 

K 



(75) 



These phase shifts have the following symmetry property: 



6±{k) = -S^{-k) for |fc| < kp- 



(76) 



The knowledge of the phase shifts allows the calcula- 
tion of the edge singularities of A{k, e) at arbitrary mo- 
menta. In the interval k E [— /cp, fci?], the edge of support 
for negative (positive) energies is located at e = eth(^) 
(e — —Ethik)). A periodic continuation of the functions 
±£th(A;) yields the edge of support for arbitary k, as de- 
picted in Fig. 4. Such thresholds, which can be con- 
structed from £th(fc) in the main region k € [—k^.k^] 
by using shifts and inversions, are called shadow bands. 
Hence, for (2n - l)kF < k < {2n + l)kF (n € Z), 
the lower and upper edges of support are at the energy 
£ = ±eth{kn), respectively, where 



kn = k — 2nkF € [—kp, kp]- 



(77) 



Let us first discuss the hole sector, e « £th(fcra) < 0, at 
arbitrary k. The configuration responsible for the edge 
singularity can again be determined by requiring energy 
and momentum conservation. It contains an impurity 
hole on mass shell near momentum fc„, which carries 
almost the entire energy e. The remaining momentum 
2nkF is absorbed by creating n low-energy holes near 
the right Fermi point and n low-energy particles near 
the left Fermi point. Therefore, one projects the fermion 
operator as follows, 



^{x) e"''=[iplipRrd oc e"=^ exp[-2in(f>{x)]d{x) , (78) 



where we used the bosonization formula 4'R,l oc 
exp[—i{±(l) — 6)]. In this sector, the spectral function is 
the Fourier transform of (^'^(0, 0)^'(.x, t)) and its thresh- 
old singularities can be calculated using Eq. (78). 

An analogous calculation is used in the particle sec- 
tor, £ « — £th(fcn) > at arbitrary k. Here, the 
threshold configuration is a generalization of the one 
used in the perturbative calculation in Eq. (38). For 
kp < k < 3kp, it contains an impurity hole at momen- 
tum — fci = 2kp — k & [—kp,kp] on mass shell, as well 



as two particles close to the right Fermi point. Similarly, 
for general (2n — l)kp < k < (2n -|- l)kp, the impurity 
carries the momentum —kn, and there are n+1 particles 
at the right Fermi point and n — 1 holes at the left Fermi 
point. The corresponding projection now reads, 

vft(a;)^e^'=-[V>JjVL]"-V4V^l?rf 

oc e**^ exp[2m^(a;)] eyip[-2i0{x)]d{x). (79) 

The spectral function can now be found by calculating 
the correlation function (\E'(a;, t)4''''(0, 0)) with the help 
of the diagonalized mobile-impurity Hamiltonian (64)- 
(65). Combining the particle and hole sectors, the spec- 
tral function near the edges of support T£th(^n) at arbi- 
trary momentum k reads, respectively, 

A{k,e) cx 0[£th(A„) ±£]kth(fcn) ± £|-'^-±W. (80) 

Here the threshold exponents are given by (Imambekov 
and Glazman, 2009a) 



Mn,±(fc) = 1 - Mn.i.fl - tJ-n,±,L, 



(81) 



where jJLn.±,R and fi„_±.L denote the contributions due to 
left and right Fermi points. They are given by 



n2 



nVK -Sa,{kn)/i2TT) 



nVK + a/VK + (5„(-fc„)/(27r) 



(82) 

(83) 



where a = i?, L = +, — . In conclusion, the spectral func- 
tion at the edge of support still displays a sharp power- 
law even for momenta far away from the Fermi points. 
However, in contrast to the linear LL theory, the expo- 
nents become nonunivcrsal and momentum-dependent. 
In particular, they depend on the microscopic interac- 
tions not only through the Luttinger parameter K, but 
also through the phase shifts (5±(±fc„), which accord- 
ing to Eq. (75) are determined by the shape of the edge 
eth(fc). 

Let us now turn to the DSF S{q,uj). The key to the 
calculation of its edge singularities is finding the correct 
projection of the density operator p{x) = '^^ {x)'ii (x) . It 
was shown in Eq. (54) that for (7 — > 0, S{q,Lo) has a 
rectangular shape and the lower edge of support is at 
—^{kp — q) ^ vq — q^ /{2fa). The configuration which 
gives rise to this threshold consists of a hole on mass 
shell and a particle near the Fermi point. Therefore, by 
continuity the threshold for general momenta < g < 
2kp is at — £th(fc_F — ?), and the density operator must 
be projected to p(x) oc Vjj(a;)ci(a;). Here, '^Jj creates a 
particle near the right Fermi point and d creates a hole 
at momentum kp — q. 

Let us generalize this argument to arbitrary momenta, 
2nkp < q < 2{n+l)kp (n G Z). Because S{q,uj) is sym- 
metric, let us focus on g > 0. In analogy to the previous 
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paragraph, the configuration responsible for the edge sin- 
gularity at momentum q now contains an impurity hole 
d at momentum 

g„ = (2n + l)kF - g e [-kp, kp]- (84) 

Moreover, it contains n + 1 particles near the right Fermi 

point and n holes near the left Fermi point. Therefore, 
the density operator for momenta near q is projected as 

p{x) ^ e^''-V^;[V^UL]"d 

oc e"?^ cxp[i(2n + l)4>{x)] eyiY>[-i6{x)]d{x), (85) 

and S{q, u) can be found by calculating the correlation 
function (p(a;, t)p(0, 0)) using the projection (85). The 
total energy of this configuration is close to — £th(5rt)- 
Hence, clS cL result one finds that 

S{q./oj) cx e[uj + eaMn)][^ + eth(9«)]"^"^'^^ 
with a momentum-dependent exponent 
A«n(«) = 1 - - A^n.L, 

"(2n+l)v^ , a , (5c«((?n)'' 



(86) 
(87) 



2\/K 



+ 



2-K 



where a = R,L = and Hn,R and /i„,L denote con- 
tributions due to left and right Fermi points. 

For momenta close to ^kp and interaction potentials 
decaying faster than the results of this section re- 

produce the universal phase shifts of Sec. II. B at arbi- 
trary interaction strength. Moreover, for V{x) oc 
the corresponding phase shifts confirm those found from 
the Bcthc ansatz solution of the Calogero-Sutherland 
model, see Eq. (158). For generic weak interaction po- 
tentials, the threshold position Sth{k) can be calculated 
perturbatively and the phase shifts d±{k) can be derived 
from Eq. (75). In this limit it is possible to recover the 
results of Sec. II. A. 

The crucial step in the calculation of the exponents 
we outlined above is the identification of the fermionic 
operator at the edges in terms of the corresponding im- 
purity operator, see Eqs. (78) and (79). A comparison 
with the solvable cases above shows that such an identi- 
fication indeed holds in the vicinity of the Fermi points 
for any interaction strength, as well as at arbitrary mo- 
menta for weak interactions. However, the identification 
of the state which corresponds to the edge of support 
may have to be modified if the interactions affect the 
impurity spectrum too much, so that even in the region 
|fc| < kp the impurity band doesn't correspond to edge of 
support. The simplest scenario would be if at some k the 
true edge of support contains a single impurity with mo- 
mentum close to k, as well as a low energy particle-hole 
pair at one of the branches. We expect that it doesn't 
happen as long as 

dsth{k) 



dk 



< V for < kp. 



(88) 



Equation (88) guarantees that the phases in Eq. (75) are 
continuous functions of momentum, and the state which 
corresponds to the edge of the spectral function in the 
basic region will contain a single impurity. We note that 
the condition (88) can be violated, e.g., for a microscopic 
model which describes ultracold fermions with resonant 
interactions (Imambekov et ai, 2010). 

D. Phenomenology of spin liquids 

Besides fermionic systems, spin chains constitute an- 
other noteworthy branch in the family of ID quantum 
liquids. There is a plethora of experimentally accessi- 
ble compounds which form antiferromagnetic Heisenbcrg 
chains and whose spin structure factor can be probed 
using neutron scattering. In addition to neutron scat- 
tering, electron spin resonance (ESR) can be used to 
probe some of the results discussed here in the pres- 
ence of Dzyaloshinskii-Moriya interactions (Karimi and 
Affleck, 2011; Povarov et at, 2011). In the present con- 
text, antiferromagnetic spin- 1/2 systems are of particular 
importance because their low-energy degrees of freedom 
often fall into the universality class of a linear LL (Gia- 
marchi, 2004). The goal of this section is to extend the 
phenomenological considerations of Sec. II. C towards re- 
sponse functions of ID spin liquids. Beyond the low- 
energy limit, various microscopic models are used for 
the theoretical description of spin-1/2 chains. The phe- 
nomenological results of this section apply to a wide va- 
riety of model Hamiltonians of the form 



H 



' ^ ^ Jn—n' 

n,n' 



QZ QZ ' 



(89) 



where S„ = {S^, S^, S^) are spin-1/2 operators located 
on the lattice sites n. The coupling between spins on 
different lattice sites is denoted by J„. The exchange 
is anisotropic for A„ ^ 1. An applied magnetic field 
in z-direction is denoted by h. Importantly, we assume 
throughout this section that the parameters J„, A„ and 
h are chosen such that the system is a liquid, i.e., the 
spectrum is gapless. 

In order to introduce some general concepts of spin 
liquids on a simpler model, we start the discussion from 
the spin-1/2 XXZ model, 



n+li 



(90) 



The XXZ model can be solved exactly using the Bethe 

ansatz (Korepin et ai, 1993; Orbach, 1958) and this as- 
pect will be discussed in more detail in Sec. III.D. At 
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A = the Hamiltonian (90) is easily mapped onto a 
model of free spinless fermions, see Eq. (94) below. For 
A = 1 the nearest-neighbor couplings are isotropic, and 
this case is referred to as XXX model. If additionally 
ft = 0, it also becomes SU{2) invariant. For |A| > 1 
and h = 0, the system is gapped and the ground state 
becomes ferromagnetic or antifcrromagnctic, depending 
on the sign of A. On the other hand, the system re- 
mains gapless for |A| < 1 and sufficiently small magnetic 
fields (Haldane, 1980; Korepin et al., 1993; Takahashi, 
2005). The low-energy sector in the corresponding pa- 
rameter range can be modeled as an LL, characterized 
by a Luttinger parameter K > 1/2. The limit i^' = 1/2 
is reached for the isotropic Heisenberg model (A = 1, 
h = 0), while the noninteracting case (A = 0) leads 
to K = 1. Certain properties of XXZ model which are 
specific to its integrability can be calculated using ap- 
proaches which combine the phenomenology with the re- 
sults of the Bethe ansatz (Cheianov and Pustilnik, 2008; 
Pereira et al, 2006, 2008, 2009), see Sec. III.D. In this 
section, however, we will focus on its generic properties. 

The XXZ model and its extensions we consider here 
preserve rotation symmetry in the xy plane. The dy- 
namic responses of this type of spin chain are therefore 
encoded in the transversal and longitudinal spin struc- 
ture factors, respectively, 

S-+{q,uj) = 5^6-'"" [ dte'-' {S-{t)S+{Q)) , (91) 
S''iq,Lo) = ^e-^«" / dte^"* (SmS^iO)) , (92) 

n ■' 

where = iS^ and we assumed an infinite num- 
ber of lattice sites. Due to the presence of the lattice, the 
quasimomentum q is bounded, < tt. Let us also men- 
tion that although we focus on the spin structure factors, 
a similar phenomenology can also be used to find other 
correlation functions, e.g., the spin-exchange structure 
factor (Klauser et al, 2011) which determines the rate of 
resonant inelastic X-ray scattering (Ament et al, 2011). 

The spin operators can be mapped onto spinless 
fermions using a Jordan- Wigner transformation (Gia- 
marchi, 2004), 



fermions if the density is measured with respect to the 
half-filled band. 

The XXZ Hamiltonian has the following fermionic rep- 
resentation. 



5+^(-l)"4e'"'*", S-^(-l)"c„e- 
1 

2' 



(93) 



where cjj (c„) is a creation (annihilation) operator for a 
fermion on lattice site n and the exponential of (pn = 
^""l^cjcj denotes a Jordan- Wigner string. The lat- 
ter is needed to ensure the proper anticommutation re- 
lation for fermions on different lattice sites, {cl^,Cm} = 
Snm- Because the operator has the simple represen- 
tation (93) in terms of the fermionic density p„ = c|,c„, 
S^^ {q, oj) is identical to the DSF of a system of spinless 



Hxxz = ^(4c„+i + h.c.) -h^(pn-^ 

n n 

+ JA^(^p„-i) (^Pn+i-iy (94) 



In the fermionic language, a nonzero A corresponds to 
a short-range interaction potential. The magnetic field 
h plays the role of the chemical potential. The Jordan- 
Wigner transformation can also be applied to the more 
general Hamiltonian (89) and leads to higher-order, pos- 
sibly non-local interaction terms in addition to Eq. (94). 

Due to this mapping on spinless fermions, many of 
the methods presented in the previous sections can be 
generalized to spin- 1/2 chains with gapless spectra, but 
with certain peculiarities. First, one sees that unlike for 
Galilean invariant systems, the kinetic part in Eq. (94) 
can have different signs of the nonlinearity. This results 
in a much wider variety of threshold behaviors. Sec- 
ond, for a general filling fraction {h ^ 0) the concept 
of an edge of support is, strictly speaking, not defined. 
This is a consequence of the presence of the lattice which 
results in "foldings" of the shadow bands of Sec. H.C 
into the reduced zone scheme. However, the contribu- 
tions to response functions from higher shadow bands are 
suppressed exponentially in their order. Indeed, in the 
language of fermions, the foldings come from umklapp 
processes which create additional particle-hole pairs and 
change the momentum of fermion system by multiples of 
the Brillouin zone period. If the Fermi wave vector is in- 
commensurate with it, the creation of pairs is capable of 
reducing the energy of the "deep" hole to an arbitrarily 
low value (Pereira et al., 2009). The lower the resulting 
energy value, the more particle-hole pairs need to be cre- 
ated. The creation of these additional particle-hole pairs 
makes the Anderson orthogonality catastrophe stronger 
and thus suppresses the corresponding contributions to 
the response functions. 

In the most natural case of /i = 0, Eq. (89) is Z2- 

invariant, which corresponds to a particle-hole symmet- 
ric Hamiltonian in the fermionic representation (94). As 
a consequence of particle-hole symmetry, the leading cor- 
rection to spectrum is cubic in momenta. This new fea- 
ture appears due to the presence of a lattice and has ram- 
ifications regarding the universal description of Sec. II. B. 
At the same time, the edge of support is still well-defined, 
since shadow bands fold onto each other (Pereira et al, 
2009); we will focus on the h = case in the bulk of 
this section, deferring to its end the consideration of the 
/i 7^ case. 
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FIG. 7 (Color online) Longitudinal spin structure factor 
S^{q, ui) for the isotropic XY model, see Eq. (100). For fixed 
q < TT, the function has a step-function threshold at the lower 
edge of support ujl{<i) and diverges as an inverse square-root 
at the upper edge of support uju{q)- 



1. XY model at zero magnetic field 

First, we will illustrate some of the generic proper- 
ties of spin liquids and their field theoretical description 
by considering in detail the isotropic XY model which 
corresponds to A = in the Hamiltonian (90). Its coun- 
terpart (94) describes free fermions. Their single-particle 
spectrum reads 



^(fc) = — J[cos(fc) — cos(fcf-)], k e [— TT, tt] 



(95) 



and at h — 0, it is particlc-hole symmetric, kp — 7r/2. 
In order to evaluate the longitudinal structure factor 
S^Y{q,uj) with the help of quantum impurity model, we 
need to find the appropriate projections of the operators 
S^. The lower threshold for S'xy(9,w) is reached at the 
energy 



^L(<z) = -e(^/2-g)-e(V2 + g). 



(96) 



The two corresponding configurations consist, respec- 
tively, of a deep hole with an accompanying low-energy 
particle, and a finite-energy particle accompanied by a 
low-energy hole (the low-energy particle and hole, respec- 
tively, are near the right Fermi point). Introducing the 
operator di creating a hole with a momentum close to 
tt/2 — q, and the operator creating a particle with mo- 
mentum in the vicinity of 7r/2-f g, we may write (Pereira 
et al, 2008) 



S„ 



i/^idi 



(97) 



Having two rather than one threshold configurations in 
the projection Eq. (97) is a consequence of the particle- 
hole symmetry. Similar to the free-fermion density struc- 
ture factor (10), the threshold frequency uJL{q) reaches 
zero at g = and q = 2k p = tt. The threshold en- 
ergy as a function of q is determined by the free-fermion 
spectrum (95), and the corresponding exponent fx~ = 0. 

The nonzero values of {q, oj) are confined to a finite 
energy window (Miiller et ai, 1979, 1981). The upper 
threshold at /i = is generated by particle-hole pairs with 



particle and hole momentum in the vicinity of (tt -I- q) /2 
and (tt — g)/2, respectively. The projection of onto 
these states reads (Pereira et al, 2009) 



S„ 



and the threshold energy is given by 

UJuiq) ^ + q/2) ~ ~ q/2). 



(98) 



(99) 



The structure factor is proportional to the convolution 
of the particle and hole spectral functions, resulting in 
the exponent ^+ = 1/2, independently of q. A direct 
calculation of the structure factor yields (Miiller et al, 
1979, 1981) 



(g,w) = 



Hit) 



(100) 



a;2]i/2 

The threshold energies and the edges of support of 
Sxyiq,!^) are depicted in Fig. 7. 

The evaluation of the transversal structure factor 
•^XY even at A = 0, is more complicated due to the 

presence of the string operators in the fermionic repre- 
sentation of S^, see Eq. (93). The Jordan- Wigner trans- 
formation shifts the locations of zero modes of S^Y{q, uj) 
to g = TT ± 2nkp, as in the conventional LL theory (Gia- 
marchi, 2004) (note that zero mode ai q — n is present at 
any value of h). The thresholds of S^yiq^^) deter- 
mined by the same configurations that define the thresh- 
olds of the spectral function A{k, oj) for spinless fermions 
at momenta k = tt — q + kp. For simplicity, we will refer 
here to the vicinity of q = tt, i.e., < k — kp <C kp. Note 
that in the case of /i = 0, which we are still interested in, 
the threshold exponents obtained below are valid for the 
entire interval |g| < tt. 

The mapping onto spinless fermions allows us to use 
the results of Sec. II. B for the threshold configurations. 
The singularity in S^Y{q,uj) at energy ui = £^{3tt/2 — 
q) — —£,{q — 7r/2) is the counterpart of the "mass shell" 
singularity in the spectral function of free fermions, see 
Sec. II. A. 2. The particle-like configuration contributing 
to the projection of spin operators takes the form 



-ikn 



exp 



dypiv) 



(101) 



where d^ creates an impurity particle at momentum kd = 
37r/2 — g. As A = 0, there is no interaction between 
the impurity and excitations near the Fermi level, and a 
straightforward calculation yields the threshold exponent 
ji^ — 1/2. Remarkably, even at A = the nonlinear 
LL theory yields a result for S^ylq, to) which is different 
from the result /ill = 3/4 obtained from an analysis of 
a linearized fermionic spectrum (Giamarchi, 2004). The 
qualitative behavior of S'xy iUustrated in Fig. 8. 

Turning now to the peculiarity in S^Y{q, cj) at the up- 
per threshold ^(Tr - q/2) - £,{q/2) = 2^(7r - q/2), we no- 
tice that the configuration representing and consist- 
ing of a minimal number of excitations should contain a 
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FIG. 8 (Color online) Transversal structure factor S^y(1:'^) 
for the isotropic XY model for q ~ n at h — 0. A divergent 
power-law with exponent jj.^ — 1/2 is found at the lower 
edge ^{3n/2 — q) = ^£,{q — 'r/2). Note that the linearized 
Luttinger liquid (LL) theory predicts only one singularity with 
a different exponent /iLL = 3/4. At higher energy, ^(tt — g/2) — 
^{q/2) — 2^(7r — q/2), the transversal structure factor has a 
step-like feature, jj,^ = 0. 



particle- hole pair at momenta fc = 7r/2 ± (tt — g)/2 and a 
low-energy particle near the Fermi level, 



d\ d2V'i? exp 



dyp{y) 



(102) 



The time evolutions of the two impurities' creation and 
annihilation operators are independent of each other and 
of the evolution of the rest of the product, which in- 
cludes the creation operator of a particle near the Fermi 
point and the string operator. The time evolution of 
the impurities is controlled by the corresponding single- 
particle Hamiltonians. The third component of the prod- 
uct can be treated by the standard bosonization. This 
way, the spin correlation function in real space-time fac- 
torizes into a product of three functions. Taking the 
proper long-time asymptotes and performing a Fourier 
transform, one finds fi^ — 0. Note that because of 
the presence of the string operator in Eq. (102), the 
structure factor S^Y{q,uj) remains finite after the drop 
at CO = 2^{n — q/2), see the sketch in Fig. 8. The 
closed-form analytical expression for the correlation func- 
tion (S'~ (<)S'(j~(0)) has been obtained recently (Zvonarev 
et al, 2009c) as a Fredholm determinant, and in princi- 
ple contains the full information about the properties of 
the spectral function S'^y('?j^)- 



2. Generic spin chains at zero magnetic field 

Moving away from the XY model to the general Hamil- 
tonian (89) introduces interactions between the Jordan- 
Wigner fermions. The location of zero modes of the spin 
liquid (g = and q = tt at /i = 0) is preserved, but 
the interactions affect the detailed shape of the thresh- 



olds col{q) and uJu{q), and generically lead to smear- 
ing of uju{q) away from the low-energy regions. The 
particle-hole symmetry inherent to the case of zero field 
brings several new qualitative features which are ab- 
sent for models in the continuum. As a consequence of 
particle-hole symmetry, the quadratic term in the band 
curvature vanishes (m = oo) and the leading spectrum 
nonlinearity at momenta close to ±kp — ±7r/2 becomes 
cubic. This makes a direct application of the universal 
results of Sec. ILB to the investigation of the singular- 
ities at uj — LOL{q) impossible. One still may use per- 
turbation theory of Sec. II. A if the interaction constants 
A„ are small, and, e.g., for the XXZ model one obtains 
li~ ~ 2A/7r (Pereira et al., 2008). If the phase shifts 
are known non-perturbatively, then the exponents for 
S^^{q,oj) and 5" '"(g, w) can be evaluated using Eqs. (87) 
and (112), respectively. For the case of isotropic exchange 
(A„ — 1), the SU{2) invariance of the equal-time correla- 
tions dictates (Giamarchi, 2004) K ~ 1/2 irrespective of 
J„. The universal values of the threshold exponents can 
be established similarly with the help of SU (2) symmetry 
in conjunction with the phenomenology. 

The quasiparticle configurations with energies near the 
upper threshold 0Jij{q), which is well defined in the low- 
energy region {q <^ tt), involve two mobile impurities. 
Generically these two interact with each other, but as 
long as ft, = and particle-hole symmetry for fermions is 
preserved, there is no interaction of the two- impurity con- 
figuration with the low-energy excitations (Pereira et al., 
2008). That brings universality to the exponents /i+ and 
The response function S^^{q,ui) is associated with 
the creation of a particle (d|) and a hole ((^2) having re- 
spective momenta Tr/2 + p ± q/2, where p <^ q. Here, 
q and p are the total and relative momentum, respec- 
tively. The relative velocity of particle and hole motion 
vanishes at p = 0, so the particle-hole interaction effect 
is strong even at infinitesimal interaction. The DSF near 
the threshold is proportional to the probability density 
w{q,p) to create a particle and a hole at the same point, 
and to the joint density of states. 



S{q,uj)cc I dp 5[e{q,p) - uj\w{q,p). 



(103) 



Here, s{q,p) = ujuil) + E^ip) is the energy of the pair 
near the threshold. The energy of the relative motion 
Er{p) =p^/{2M) depends on q via 1 /M cx d^uu /dq^ oc 
q, while w{q,p) oc \il^p{x = 0)^ is found from the solution 
ipp{x) = cos (p I a; I +1}) of the Schrodinger equation for the 
relative motion in the pair, Hr^p^x) = Er{p)ipp{x). Here 



l~Lr 



'2Md^ 



Ul2S{x) 



(104) 



^ is the effective interaction between the 
< 0). The solution 
at small \Er{p)\ < \M\Uf2 « q^ yields ippix = 0) oc 



and U12 oc q 
particle and hole (note that M 
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^/Er{p)/{\M\U1^). Therefore, using Eq. (103) at < 
'^u{q) — w <C \uiu{(i) — vq\, one obtains 

S^^{q,uj) oc q-'^/'^e[iOu{q) - cjW^u{q) - w 

+ regular terms. (105) 

The momentum independent threshold exponent /i+ = 

— 1/2 differs from its value at A„ = 0: interactions cause 
a discontinuous change in the edge exponent (Pereira 
et al., 2008). 

To describe 5* '"(g, w) near the upper threshold, we 
may use projection of spin operators of Eq. (102). The 
real space-time correlation function factorizes into a 
product of two functions. The first one describes the 
evolution of the particle-hole pair, identical to that ap- 
pearing in the S^^{q,uj) correlation function. The sec- 
ond function in the product comes from the evolution 
of the string operator multiplied by the creation oper- 
ator of a fermion near the kp = 7r/2 point; that func- 
tion depends on K. The resulting convolution yields 
M+ = -3/2 + l/(2X). 

As we may see from the derivation of Eq. (105), 
the asymptotic behavior of S^^{q,u!) at small positive 
'-^u{q) — w is independent of the sign of the particle- 
hole interaction. However, at [/12 > (we assume that 
M < 0, as it is in the XXZ model), a bound particle-hole 
state is formed. It would show up, for instance, as an 
additional peak in S^^{q,oj) at some value ui > 0Ju{q). 
Generically the energy of such bound state is finite for 
all momenta (as is illustrated by the XXZ model (Pereira 
et al., 2009)), and can be of the order of bandwidth for 
strong interactions. The possible bound state peak as 
well as the upper threshold ijJjjW) (away from special 
points g = and q — ±7r) appear on the background of 
the spectral continuum, and will be generically smeared 
out. 

Considering the upper threshold, we find /z+ = /Lt+ at 
SU{2) invariance, K = 1/2. This is not a coincidence 
but rather a consequence of SU{2) symmetry which en- 
forces S'='={q,oj) = Syy{q,uj) = S''^{q,oj). Now we will 
use that constraint on the correlation functions in order 
to obtain their threshold behavior at = LJi^(q). To- 
gether with the constraint, phenomenology alone is suf- 
ficient to establish momentum-independent /Uj = 1/2 as 
the universal value independent of the microscopic de- 
tails. Using the mobile impurity Hamiltonian (64)-(65) 
and the projections similar to Eq. (97) and (101), the 
edge exponents of both S ^{q,Lj) and S'^^(g,cj) can be 
calculated independently as functions of the phase shifts 
^±((7). Moreover, since the system is half-filled, umklapp 
scattering is allowed and can produce terms of the form 
[V'lV'IV'fl^-R]™ ("^ G ^) ill the projections of and S^. 
The edge exponents of the longitudinal and transversal 
spin structure factors, //7rn ^'^^ l^xrm become functions 
of m. As a consequence, the structure factors at the edge 



are characterized by sums of the form, 

S^^{q,cj)<xJ2\^-^L{q)\-^'-, 

m 

S-+{q,oj)<xY,\^-^L{q)\-''^-. (106) 

n 

A calculation similar to Eqs. (81)-(83) combined with 
K = 1/2 leads to 

_ 1 /_4to- 1 (5+ -h<5_y 

_ _^ If 6+-S_ Y If 4n + l ^ S+ + 6_ Y 
^-'^ ~ 2\ 2n ) 2\ ^ ^ 2n ) ■ 

As a consequence of SU{2) invariance the entire series for 
and have to coincide at each q. This require- 
ment leads to the momentum-independent phase shifts 
S± = =p7r/\/2 and exponents (Imambekov and Glazman, 
2009a), 

S'%q,oj) = lS-+{q,w) oc \uj - UL{q)\-^/\ (108) 

It coincides with the result of linear LL theory for 
5-+(g,w) at g ^ TT, /iLL = 1 - 1/(4/^) = 1/2 (Gia- 
marchi, 2004). Note that similar to equal-time correla- 
tors (Affleck et al, 1989; Singh et al, 1989), Eq. (108) 
will also have logarithmic corrections due to existence of 
marginally relevant terms in the Luttinger Hamiltonian. 

3. Generic spin chains at finite magnetic field 

A finite magnetic field {h ^ 0) violates the Z2 invari- 
ance of the spin liquid and destroys the particle-hole sym- 
metry of the corresponding Hamiltonian in the fermionic 
variables. It shifts the Fermi points ±kp away from 
±7r/2; their new positions ikp = ±7r{{S^) + 1/2), are 
fully determined by the induced magnetization (S^) . The 
positions of two of the low-energy regions, q = for 
S^^{q,u)) and g = tt for 6'~+(g, w), do not shift. The 
other two regions are shifted from q = tt to q = ±2kp 
in S^{q,oj), and from q = to q = ±{n — 2fcj?) in 
S '"(g, w) (Giamarchi, 2004). As we already mentioned 
earlier in this section, folding of the shadow bands at 
h makes the thresholds in the response functions 
at arbitrary q ill-defined. However, sharp threshold lines 
in the (q, ix))-plane remain intact in the vicinities of the 
w = singularities of the response functions. Moreover, 
now the expansion ink—kp of the elementary excitations 
spectrum does contain the quadratic term, 1/rh 7^ 0, 
so one may use the universal nonlinear LL theory, see 
Sec. II. B. This theory is valid as long as the interaction 
potential J„A„ decays faster than l/n^. For k ikp, 
the phase shifts d±{k) reach the universal values (48). 
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The threshold behavior of S^^{q,U!) at q 
cal to Eq. (54), 

^2 



5-(g,a;) = 



K\m\ 



2|m 



z-j- — |w — vg'l 



is identi- 



(109) 



and the exponents = are independent of the interac- 
tions. Here v and K are the characteristics of the hnear 
LL, and 1/m = {v/K^/^)dv/dh + v^/{2K^^^)dK/dh is 
phenomenologically related to these quantities (Pereira 
et al, 2006). Similarly, for g' ^ tt (Imambekov and Glaz- 
man, 2009b) 



S ~^{q,uj) oc 



v\q — 7r| ± 



2m 



(110) 



Note that the exponents /j.^ = ±{1/VK) - 1/{2K) 
differ from the linear LL prediction /xll = 1 — 
l/(4i^)(Giamarchi, 2004). 

We also note that the region of validity of Eqs. (109)- 
(110) in q becomes narrow (oc h), as the magnetic field 
decreases. The interesting question of whether there is a 
universal description of the crossover from the results of 
Sec. II. D. 2 to the universal description at nonzero mag- 
netic fields remains to be addressed. The behavior of the 
response functions in the {q, to) plane away from the zero- 
frequency special points is generally model-dependent. 



E. Bosonic systems with or without spin 

The interest in interacting ID bosonic systems, sec the 
recent review by Cazalilla et al. (2011), has soared in past 
years mostly thanks to an increasing variety of experi- 
ments. The ID bosons can be realized using Joscphson 
junction arrays (Fazio and van der Zant, 2001), helium 
in nanopores (Del Maestro and Affleck, 2010; Del Mae- 
stro et at, 2011; Wada et at, 2001), as an effective de- 
scription of magnons in magnetic insulators (Giamarchi 
et al, 2008), and most prominently in systems of ultra- 
cold atomic gases. In the latter case, one quite naturally 
often ends up with exactly solvable models, and these 
will be reviewed in detail in Sec. III. Here we will fo- 
cus on a general phenomenological description of bosonic 
systems, both spinless and spinful. 

Let us start with the discussion of the spinless case. 
In the absence of interactions, ID bosonic systems are 
very different from fermionic ones. The former have 
a quadratic spectrum of excitations, while the excita- 
tion spectrum of the latter is linear. However, even for 
infinitesimally small repulsive interactions, bosons ac- 
quire a finite compressibility and their spectrum also 
becomes linear. The description of such bosons within 
the linear LL theory is not that different from the de- 
scription of fermions (Efetov and Larkin, 1975; Gia- 
marchi, 2004; Haldane, 1981a). The slowly varying com- 
ponent of the bosonic field is given by \l/(a;) ~ Q^^i^) _ 



The difference between this equation and, e.g., the 
bosonized expression for the right-moving fermionic field 
(19) is the factor exp [i{-!Tpx — cp)], which is nothing but 
a bosonized version of the Jordan- Wigner string opera- 
tor exp[z7r \E'^(y)\E'(2/)dy] discussed in Sec. II. D. As 
a consequence of that, the edge of support of the bosonic 
spectral function, which we denote by A{q,e), is shifted 
hy kp compared to fermionic case. For e > 0, it coincides 
with that of the DSF, and for historical reasons we shall 
denote it by £2(9) > in the region < q < 2kp. For 
weakly interacting bosons, £2(9) corresponds to excita- 
tion spectrum of dark solitons, see Sec. III.B. 

The Jordan- Wigner string docs not appear in the 
bosonized form of the density operator. Therefore, the 
low-energy projection of the density operator does not 
depend on statistics of the particles (bosons or fermions) . 
If one considers the DSF, it is not possible to distinguish 
spinless bosonic systems from spinless fermionic systems. 
Hence, the phenomenological relations of Sec. II. C de- 
termine the DSF exponents (87) as well as the phe- 
nomenological phase shifts (75), assuming that the ef- 
fective threshold position is identified as 

£th(fc = kF-q) = -£2(9) < for < g < 2kF. (HI) 

The same phase shifts and effective Hamiltonian also de- 
termine the edge exponents of the spectral function, when 
combined with proper operator identifications similar to 
Eq. (101) for spin chains. Similarly to the equal-time field 
correlator, only contributions arising from the Jordan- 
Wigner string from the vicinities of the Fermi points 
have to be treated differently compared to the fermionic 
case. Near the thresholds in the regions 2nkp < q < 
2{n + l)kp, the exponents for £ ^ are given by, (Imam- 
bekov and Glazman, 2009a) 



M^,±(9) = l-<±,ij-<±,i, (112) 
,-,a = [(2n + 1)VK/2 + 5„(g„)/(27r)] ' , 

M^+.a = [(2n + 1)Vk/2 + a/VK - (5„(-g„)/(27r) 

where a = R, L = +, —, and g„ = (2n -I- l)kp — q. The 
notation is illustrated in Fig. 12. Since the operators 
carry the same Jordan- Wigner string as the boson 
creation-annihilation operators, Eqs. (112) can also be 
used to express the exponents of S ''((/, w) in terms of 
the phase shifts. 

Next, let us consider a system of interacting bosons 
which possess (iso)spin, and we will focus on the discus- 
sion of the SU{2) symmetric case. Such systems in quasi 
ID configurations have already been realized with ultra- 
cold atomic gases (Higbie et al, 2005; McGuirk et al, 
2002; Sadler et al., 2006). In addition to the DSF and 
the spectral function, we will be also interested in the lon- 
gitudinal and transverse spin structure factors, 5^^(g, w) 
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and S^^ {q, to) respectively. These are defined by 

= J dxdte'^'^*-'^''^ {S^{x,t)S'^{0,0)) , 

S''{q,u) = y"da;die'('"*-«^)(5^(x,i)5"(0,0)), (113) 

where S'^{x) = S'^{x) ± iSy{x). In terms of the physical 
particles, the spin density S{x) = {S^{x), Sy{x), S^{x)) 
reads 



S{x) = Y,'^l{x)S,a,^,,{x), 



(114) 



where Sera' denotes the vector of spin matrices (half of 
the Pauli matrices for spin-1/2). 

Bose statistics requires that the total TV-particle wave 
function, which consists of a spatial part and a spin part, 
must be symmetric under exchange of the positions and 
spins of any two particles. Eisenberg and Lieb (2002) 
found that the unique spatial component of the ground 
state wave function is symmetric, so the spin wave func- 
tion has to be symmetric as well. As a consequence, the 
ground state of the spinful Bose system is ferromagnetic 
(we will assume that the magnetization is pointing in 
+z direction). This is in striking contrast to the case 
of spinful fcrmions, where the ground state is usually 
a spin singlet (Lieb and Mattis, 1962). Due to the spin- 
polarized ground state, the bosonic system responds very 
differently to external perturbations which conserve the 
total spin and perturbations which change it. For sim- 
plicity, in what follows we will discuss only the case of 
(iso)spin— 1/2. 

The longitudinal spin structure factor S^^{q,uj) does 
not change the total spin and is simply proportional to 
the DSF S(q,uj) for all {q,uj). The transverse spin struc- 
ture factor S ^iq,ui) vanishes because the system is al- 
ready in a state with largest possible S^, and the action 
of the operator S+{0) = *|(0)*4,(0) destroys it. On the 
other hand, 5"' {q,uj) involves the action of the oper- 
ators S~{0) = «'|(0)*^(0) and S+{x) = [S'-(.T)]t that 
flip the spin of a single boson. The excitation of lowest 
energy which is created when S~ (0) acts on the ferromag- 
netic ground state is a magnon. In an SU (2) invariant 
system, the magnon spectrum towards small momenta 
is quadratic, comil) = 9^ /(2m*), where m* denotes the 
effective magnon mass. Due to the quadratic spectrum, 
the threshold energy at small q of the transversal struc- 
ture factor is always lower than that of S^^{q,u)) and 
S{q,LLi). Note that the quadratic form of the magnon 
spectrum makes the calculation of the transversal re- 
sponse function based on the linear LL theory impos- 
sible (Akhanjee and Tserkovnyak, 2007; Kamenev and 
Glazman, 2009; Matveev and Furusaki, 2008; Zvonarev 
et al, 2007, 2009a). 

A magnon can be thought of a spin-down impurity 
moving in a liquid of spin-up particles. Unlike the mo- 



bile impurities introduced in previous sections as an effec- 
tive description, these magnons can be separately exper- 
imentally addressed, as has been demonstrated recently 
in experiments with ultracold gases (Catani et at, 2012; 
Palzer et al., 2009). An effective Hamiltonian similar to 
Eqs. (64)-(65) can be used to describe the singularities 
in the spin structure factors, if one introduces an oper- 
ator d){x) creating a bosonic spin-down impurity with 
momentum near q and energy near LOjn{q). Then the 
projection of the operator S~{x) onto the three-subband 
model reads 



5-(a;)cxe-^«"d^(a;)e*«("). 



(115) 



The phcnomcnological approach of Sec. II.C can also 
be directly generalized, and leads again to explicit 
predictions for the transverse spin structure exponent 
at the magnon spectrum fimiq) which were presented 
by Kamenev and Glazman (2009). The general expres- 
sions simplify somewhat in the limit \q\ <^ m*v due to 
the quadratic spectrum, and in order 0{q'^) are given by 



Kq^ 



l+(^)'(3 + 4a + <72) 

\m*vJ 



where a = —p/{2'm*)d'm* /dp. 



F. Fermionic systems with spin 

In this section, some of the concepts introduced in the 
previous sections will be applied for the exploration of 
the properties of spin-1/2 fcrmions in ID. Compared to 
spinless fermions, a complication arises due to the spin- 
charge separation encountered in interacting ID systems: 
a generic density-density interaction couples fermions of 
opposite spins and lifts the degeneracy between charge 
modes (excitations symmetric in spin-up and spin-down) 
and spin modes (excitations antisymmetric in spin-up 
and spin-down) present in the noninteracting system. 
Within the linear LL theory, the eigenstates of the inter- 
acting system arc then linear combinations of spin-up and 
spin-down excitations and can be interpreted as density 
waves which carry either only spin or only charge. These 
two types of excitations propagate at different velocities, 
Vs and Vc, respectively. 

The phenomenon of spin-charge separation is the hall- 
mark of the linear spinful LL theory (Dzyaloshinskii and 
Larkin, 1974): the fermionic Green's function is a prod- 
uct of two parts which describe excitations propagating 
with velocities it?;, and ±Uc, respectively. The fermionic 
spectral function which can be obtained from the Green's 
function by Fourier transformation thus has power-law 
singularities at the spin and charge modes (Meden and 
Schonhammer, 1992; Voit, 1993a, b). The spin and charge 
density structure factors, on the other hand, are delta- 
functions reflecting the linear spectra of the two types of 
bosonic modes. 
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The existence of two distinct excitation energies for a 
given momentum has already been observed in a num- 
ber of experiments with sohd state systems (Auslaender 
et al, 2005, 2002; Jompol et al, 2009; Kim et al, 2006). 
Moreover, experiments using ID ultracold fermions (Liao 
et al, 2010) could allow the observation of spin-charge 
separation in real space (KoUath and SchoUwock, 2006; 
Recati et al, 2003). 

In this section, we will review recent theoretical 
progress (Pcrcira and Sola, 2010; Schmidt et al, 2010a, b) 
in understanding the notion of spin-charge separation in 
ID models beyond the linear LL paradigm. This question 
has received some attention over the years. The existence 
of a certain form of spin-charge separation at all energy 
scales is implicitly built into the structure of the exact 
solutions of spinful integrable models which will be dis- 
cussed in more detail in Sec. III. For instance, Ogata and 
Shiba (1990) demonstrated that the eigenstates of the 
strongly interacting ID Hubbard model (see Sec. III.D) 
factorize into charge parts and spin parts. The former 
are identical to the eigenmodes of free spinless fermions, 
whereas the latter are solutions of the XXZ model. In an 
extension of this approach called pseudofermion dynam- 
ical theory (Carmelo et al, 1999; Carmelo et al, 2005, 
2004. 2006), some of the results of the present section 
have been envisioned before for the integrable ID Hub- 
bard model, but its field-theoretical basis and applicabil- 
ity for general nonintegrable models remain unclear. 

The approach of Ogata and Shiba (1990) can be ap- 
plied even without an underlying lattice or integrability, 
because strong repulsion makes the charge mode stiffer 
while softening the spin excitations. The enhanced rigid- 
ity of the charge spatial structure allows one (Matveev, 
2004a, b; Matveev et al, 2007a, b) to invoke the notion 
of a Wigner crystal (Wigner, 1934), in which fermions 
arrange on a ID lattice, as a starting point for the con- 
sideration of the dynamics, notwithstanding the absence 
of long-range crystalline order in ID. The dispersion of 
the spin excitations is suppressed in the amplitude of 
the repulsion potential. In the limit of strong repulsion 
it is tempting to dispense with the spin dynamics alto- 
gether, and to assume that arbitrarily oriented spins are 
attached to the sites of the Wigner crystal and that its 
excitations are sound waves. Such a model was dubbed 
spin-incoherent LL. The bandwidth of the spin excita- 
tions is zero in the spin-incoherent LL, and at any tem- 
perature the real-space electron Green's function decays 
exponentially in the spatial coordinate. This reflects the 
absence of any order in the spin system in the absence of 
exchange interaction between the spins (Cheianov and 
Zvonarev, 2004a,b; Fiete, 2006; Fiete, 2009; Fiete and 
Balents, 2004). We refer the reader to an excellent re- 
view of spin-incoherent LLs for details (Fiete, 2007) . 

In contrast to the notion of the spin-incoherent LL, 
here we concentrate on the generic case of comparable 
(but different) velocities of spin and charge excitations. 



As we shall sec, the concept of spin-charge separation sur- 
vives, albeit in a modified form, beyond the low-energy 
limit and the assumption of a linear spectrum. In the 
following, we shall focus on the case of repulsive interac- 
tions. In this case, the system remains gapless in both 
charge and spin sectors, and charge excitations propagate 
faster than spin excitations (vc > Vg). The parameter 
Kc characterizing the linear LL of charge excitations is 
in the range < < 1- For simplicity, we will only 
consider systems without Zceman splitting, although a 
small magnetic field might be beneficial to detect spin- 
charge separation in experiments (Grigera et al, 2004; 
Rabello and Si, 2002). The resulting SU{2) invariance 
generically leads to the spin LL parameter Kg = 1. 

The investigation of a nonlinear spectrum in the 
bosonic language runs into similar problems as in the 
spinless case (Sec. II. B). A nonlinearity leads to interac- 
tions between the bosonic modes. To lowest order, the 
corresponding self-energy diverges at energies uj = Vc^gk- 
It is still possible to investigate certain properties us- 
ing the bosonic approach (Brazovskii et al, 1993, 1994; 
Vckua et al, 2009), e.g., the dynamic response fimc- 
tions sufficiently far away from the thresholds (Pereira 
and Sela, 2010; Teber, 2007), but closer to the thresh- 
old perturbation theory fails. An efficient resummation 
scheme which eliminates these divergencies has not been 
developed yet. Therefore, it is easier to generalize the 
fermionic quasiparticle methods used in the previous sec- 
tions to spinful systems. The path to a fermionic descrip- 
tion is similar to the one followed in Sec. II. B for spinless 
fermions. For a linear fermion spectrum, the Hamiltonian 
becomes a sum of independent charge and spin parts, 
i?LL = He + Ha . The charge term He has the form of the 
linear LL Hamiltonian (41) and it is characterized by the 
charge velocity Vc and the Luttinger parameter K^- Its 
eigenmodes are bosonic charge density waves. The spin 
part Hg also contains a quadratic term of the form (41), 
characterized now by a different Luttinger parameter Kg 
and the spin velocity Vg. In addition, however, Hg gener- 
ally contains a sine-Gordon term, which describes spin- 
flip scattering between the physical fermions. For repul- 
sive interactions it is marginally irrelevant and vanishes 
in the limit of small bandwidth (Solyom, 1979). Then, 
the eigenmodes of Hg are bosonic spin density waves. 

Within the linear LL theory, the physical fermions 
can be expressed is terms of left- and right-movers using 
*<^(x) = e'^^''^R„{x)+e-^'''"=<^La{x) where a =t,4. de- 
notes the spin. The bosonization rules for the fermionic 
operators are given by (Giamarchi, 2004) 

^acr{x) (X e-i/^[°"t>cix)-ea(x)+aa^^{x)-a0s{x)] ^ (116) 

where a = R, L = +,—. Because H^ and Hg are formally 

identical to linear LL Hamiltonians, it is possible to ex- 
press them in terms of free quasiparticles by generalizing 
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Eq. (45), 



H. 



j dx : ¥^^{x)W^o.c{x) :, 

H, = -J2 ic^v, j dx : §L WV*a«(a;) :, (117) 

where a = R,L = +,—. The fermionic quasiparticles 
^ai^{x) [y = c, s) describe low-energy charge and spin ex- 
citations, which wc shall refer to as holons and spinons, 
respectively. As in the spinless case, see Eq. (46), the 
physical fermions {x) can be expressed in terms of the 
quasiparticles by using the bosonization identity (116), 
then rescaling the bosonic fields and finally refermioniz- 
ing them. For the right-movers, the result reads 

■^R^{x) CX ^ Rc{x)Fr^{x)^ r,{x)Frs{x), 

<^R^{x) (X ^R,{x)Fn,{x)¥^^{x)F],^{x). (118) 
The string operators are given by 

Fr,,{x) = exp |-i j dy [5+^pR,,{y) + 5_^pLi.(y)]| , 

(119) 

where Pau{x) = ^^^{x)'^ av{x) {a = R,L, v = c,s) de- 
note the right- and left-moving quasiparticle densities, 
respectively. Note that Eq. (118) contains charge and 
spin string operators, although we assume only Kc ^ 'i. 
The appearance of the string operators in Eq. (118) is an 
inevitable consequence of the attempted "splitting" of a 
spinful fermion into two spinless fermionic excitations. 
The phase shifts are determined by Kc, 




^ = 1- — — =0 
2-K v^' 2-K 



(120) 



Note that the Klein factors, which are required to ensure 
the correct fermionic anticommutation relations, were ne- 
glected in Eq. (118). This is justified as long as it is used 
to calculate expectation values of operators which con- 
serve charge as well as spin. This is indeed the case for 
all dynamic response functions. 

According to Eq. (118), the creation of a physical spin- 
up (spin-down) particle leads to the formation of a holon 
and the creation (annihilation) of a spinon. The string 
operators Frv{x) reflect the shake-up of the spinon and 
holon Fermi seas due to the addition of the physical 
fermion. The phase shifts (120) characterize the strength 
of the shake-up. Keeping the spectrum of fermions linear 
as in Eq. (117), and using Eqs. (118)-(120), one may read- 
ily reproduce the known results for the spectral function 



of fermions in a spinful linear LL (Schmidt et al, 2010b). 
The spectral function A{k,s) for fixed k > kp now dis- 
plays two peaks: divergent power-law singularities are 
located at the energies of the spinons and the holons, 
e = Vs{k — kp) and e = Vc{k — kp)- Within the linear 
LL theory, the exponents at both thresholds are deter- 
mined by the phase shifts 5±,j and thus by the Luttingcr 
parameter K^. (Meden, 1999; Meden and Schonhammer, 
1992; Voit, 1993a,b). 

A nonzero curvature of the physical fermion spectrum 
generally leads to interactions between the spinons and 
holons. Moreover, it also bends the spectra of both the 
spinons and the holons. For k —>■ kp, the leading correc- 
tion to the holon spectrum is quadratic. 



lc(fc) « Vc{k - kp) + 



(fc - kpf 
2^ 



(121) 



The effective mass m can be determined phenomenolog- 
ically and the result resembles the spinless case (50), 



d 



1 



rh y/2Kc dfJ- 



(122) 



The shape of the spinon spectrum ^,s{q)^ on the other 
hand, becomes rather different. From Eq. (118), it can 
be inferred that the spin- up/spin-down symmetry of the 
physical fermions entails particle-hole symmetry of the 
spinons. Therefore, there is no quadratic term in ^s(g). 
Instead, for g — )• the leading nonlinearity is cubic. 



is{q) ~ VsQ - for \q\ < kp, 



(123) 



with 7 > 0. The spinon spectrum near the Fermi points is 
therefore similar to the low-energy spectrum of the XXX 
model in zero magnetic field encountered in Sec. II. D. In 
both cases, the shape is a direct consequence of SU{2) 
symmetry. 

For repulsive interactions one has Vs < Vc, so for a 
given momentum k « -iikp exciting a spinon requires 
less energy than exciting a holon. Since the thresholds in 
all dynamic correlation functions should be continuous, 
they coincide with a shifted spinon mass shell S^sil) for 
arbitrary momenta. The precise shape of ^s(9) away from 
Fermi points depends on the microscopic details of the 
interaction and is generally unknown. 

The power-law singularities characterizing the dy- 
namic responses beyond the linearized theory can be de- 
termined by generalizing the mobile-impurity Hamilto- 
nian (56) to the spinful case. Within the linearized the- 
ory, all threshold singularities in the dynamic response 
functions are characterized by configurations where a 
spinon carries almost the entire available energy. Hence, 
the quantum numbers of the effective impurity should 
coincide with that of a spinon near the Fermi points. 
By continuity, this remains true also for momenta away 
from Fermi points. For a spinon impurity at momentum 
Qd = k — kp where k € [—kp, kp], we use the projection 
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^iis{x) — >■ e^'^'^^d{x) + ijjftsix). The spinons and holons 
near the Fermi points, as well as the impurity at momen- 
tum are then described by 



a 



dx 



■ '4>ic{x)^i>occ{x) : 



dx 



■ V'L(a;)VV'as(a;) 



Hd = j dxdlix) [is{q) - WdV] d,{x), (124) 

where a = R,L = +,— and Vd = d^s{<l)/dq. The inter- 
action between the impurity and the low-energy spinon 
and holon degrees of freedom leads to 

Hint = dx ^ ^ Va^{k)pai,{x)dl{x)ds{x). (125) 

a=R,L v=c,s 

The term Hint can be removed by using a unitary trans- 
formation similar to Eq. (71). 

The impurity causes phase shifts of the quasiparticles 
near the Fermi points, 



(126) 



where a = R,L = +,— and v = c,s. Let us first in- 
vestigate these phase shifts in the limit fc — kp. For 
interactions between the physical fermions which decay 
faster than oc 1/x^, the interaction potentials Vav{k) ful- 
fill VmA^) ^ (^ ^ kp)^. The denominators in Eq. (126) 
are determined by the cubic shape of the spinon spectrum 
(123). In the limit k — )• kp, one finds Vd±Vc — const, 
and Vd + Vg const., so the phase shifts AS±c and A(5_,, 
vanish. In contrast, Vd — Vg <x. {k — kp)^, so the phase 
shift AS+s remains nonzero even for A; — )• kp. 

This is a striking contrast to the spinless case, where 
all phase shifts vanish in the limit k ^ kp, see Eq. (56). 
The latter is a consequence of the quadratic spectrum of 
the spinless fermions near kp (Imambckov and Glazman, 
2009b). The spinon spectrum (123), on the other hand, 
is cubic due to SU (2)-symmetry. Thankfully, the phase 
shifts A5±s{k) can be determined for arbitrary momenta 
by exploiting the SU (2) symmetry. 

In order to determine these phase shifts, we use the 
same argument that led to the spin structure factor of 
the XXX model in Eq. (108): for an SU{2) symmetric 
system, the threshold exponents of the two components 
of the spin structure factors S ^{q,uj) and S^^{q,oj), de- 
fined in Eq. (113), have to coincide at arbitrary momenta 

q- 

As previously, the first step in calculating these ex- 
ponents consists in identifying the threshold configura- 
tion and projecting the operators S^{x) and S^{x) ac- 
cordingly. Then, the mobile impurity Hamiltonian (124)- 
(125) can be used to calculate the threshold exponents as 



fimctions of A5au{k). The requirement of identical expo- 
nents for S ^(q,ui) and S^''{q,ui) unambiguously entails 
that for arbitrary momentum k (Pereira and Sela, 2010; 
Schmidt et al, 2010a) 



5±s + AS±s{k) = 0. 



(127) 



Having fixed the phase shifts A^-i-j, by Eq. (127), 
only the values of A5±c{k) are left to be determined. 
For integrable models, A6±c can be extracted exactly 
from the finite-size corrections of the Bethe ansatz spec- 
trum (Carmelo et at, 2008; Essler, 2010) similar to the 
procedure discussed in Sec. II. G. For generic Galilean 
invariant models, a generalization of the phenomenolog- 
ical approach of Sec. II. C is possible. It relates AS±c in 
the interval |A:| < kp to the shape of the edge of support 
e^(fc) = fs(fc -kp) <0 (Schmidt et al., 2010a): 



k — kf 



AS±c{k) _ _^ ms/Kl 



± \fKr 



2 3ea(fc) , Oes(fc) 



dk 



27r 



2^2 ( 



dk + mif, 



(128) 



where m is the mass of the physical fermions. In the limit 
k — >■ ±fcj?, the low-energy expansion of the spinon spec- 
trum ^s(g) leads back to A(5±c = 0. This is a consequence 
of the conventional spin-charge separation of the linear 
LL theory. The equations (127) and (128) fix all param- 
eters in the mobile impurity Hamiltonian (124) and are 
thus the key for the calculation of edge exponents in the 
various dynamic response functions. 

Let us start with the spectral function. For the calcu- 
lation of the threshold behavior of A{k, e) for < A; < kp, 
one needs to consider the response of the system to the 
addition of a physical spinful fermionic hole with a given 
momentum k and the threshold energy es(fc) < 0. The 
incoming hole creates a spinon on the mass shell, which 
absorbs the entire energy, as well as a holon at a Fermi 
point. That spinon is protected from decay by energy 
and momentum conservation. If the energy of the in- 
coming hole is slightly below es(fc), the excess energy 
is used for the creation of low-energy particle-hole pairs 
in the holon sector. However, due to Eq. (127) it can- 
not be used for the creation of additional spinons! Since 
the velocity of the spinon impurity is smaller than the 
holon velocity, Vd = des{k)/dk < Vc, the spinon becomes 
spatially separated from the holons. While this effect 
is reminiscent of the conventional spin-charge separation 
of the linear LL theory, it should be pointed out that 
this new form of spin-charge separation survives for ar- 
bitrary momenta only for energies close to the threshold. 
Finding the threshold exponents of the spectral function 
A{k,s) at fc > requires a projection of the fermionic 
operator 'J'/j(j(x) in Eq. (118). Due to SU{2) symmetry, 
the result does not depend on the choice of a =t, i- For 
< k < 2kp, one can therefore project 
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FIG. 9 (Color online) Spectral function A{k, e) for spinful 
fermions {\k\ < kp), see Eq. (130). For repulsive interactions, 
the edge of support is determined by the spinon spectrum 
es{k). The power-law singularity at the holon mass shell ^c{k) 
is generally broadened away from Fermi points. 



For general (2n - l)kF < k < (2n + 1)/cf (n G Z), 
the excess momentum can be used to create additional 
particle-hole interbranch pairs similar to Sec. II. C. The 
mobile- impurity Hamiltonian now allows a calculation of 
the spectral function and the result is, 



A{k,e) oc 9[es{kr,) ± e]\es{kn) ± e\ 



(130) 



where k„ — k ~ 2nkp e [~kp,kp]. The location of the 
edges and the respective exponents are sketched in Fig. 9. 
The momentum-dependent exponents are given by 




(131) 



where A.S±c = A6±c{kn) and m± = (n + 1/2 ± 
1/2) mod 2. 

Farther away from the threshold, a second peak 
emerges at energies which correspond to the holon mass 
shell, £ sa Cc(fc)- Near this energy, the incoming physical 
particle triggers the formation of a holon on its mass shell 
as well as a low-energy spinon. Away from the Fermi mo- 
mentum, interactions between holons and spinous gener- 
ally lead to a nonzero decay rate for holons and thus a 
smearing of this peak. This will be discussed in more 
detail in Sec. IV.A. 

Let us conclude the discussion of the spectral function 
by a closer look at momenta k kp. In that limit, 
AS±c = and the exponent (131) coincides with the 
prediction of the linear LL theory, Hq _ — 1 — l/{4:Kc) — 
Kc/4: (Giamarchi, 2004). At k ^ kp, the peak at the 
holon mass shell becomes sharp. However, the exponent 
1 - [S-J{27r)f - [S+,s/i2TT) - 1]' - [5-,./(27r)]' found 



in the vicinity of the holon mass shell, \e — ^c{k)\ <C 
{k — kpY / (2rfi), is different from the LL prediction, just 
as in the case of spinless fermions. At larger detunings 
from the threshold, the power-law behavior with the LL 
value of the exponent is restored (Schmidt et al., 2010b). 
Unlike in the case of spinless fermions, here the crossover 
function between the two asymptotes of A{k,e) is not 
known. The spectral function A{k, e) and notations for 
exponents are illustrated schematically in Fig. 9. 

Next, we turn to the discussion of the spin structure 
factors. Due to SU{2) symmetry, they are related by 
S''^q,oj) = lS-+{q,uj). The calculation of S'^^(g,w) 
requires a projection of the spin density operator S^{x). 
The configuration of lowest energy for the momentum 
< q < 2kp contains a single particle-hole pair in the 
spinon sector. 



Rs- 



(132) 



For general momenta 2nkp < q < 2{n + l)kp, the projec- 
tion contains n additional low-energy particle-hole pairs 
which carry the momentum 2nkp. Using the mobile im- 
purity Hamiltonian, it can be shown that 



S''{q,u)^0[iu~\es{qM^~\^s{qd)\] 



(133) 



where qd ~ {2n + l)kp — q and the exponents are 



DSF 



nyTKc 
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27r 



V2 ^ ^ 

(134) 



The phase shifts are taken at momentum q^, i.e., 
AS±c = Ad±c{qd). Note that at g — the scattering 
phase shifts AS±c 0. Therefore, the exponent of the 
spin structure factor approaches the universal value 1/2, 
coinciding with the respective spin liquid exponent of the 
spin-1/2 XXX chain, see Sec. II. D. Similar to the spin 
chain model, the region of frequencies |w— |es((7rf)|| where 
Eq. (134) is applicable, shrinks as g^. The latter param- 
eter defines the width of the peak in the spin structure 
factor. Its detailed structure has not been investigated 
yet. In the linear LL, it is replaced by a i5-function at 
uj = Vsq. 

The quadratic dispersion of the holon spectrum leads 
to AS±c — 7- for A: — >■ zLkp and results in a rectangular- 
shaped peak in the charge DSF S{q,uj) similar to the 
case of spinless fermions, see Eq. (27). For uj w Vcq and 
|g| <^ kp, the result up to order q^ reads 



S{q,uj) 



2mKc 



2m 



Vcql 



(135) 



with an effective mass m. Because the width of the peak 
Suj^q) — q^/rh is proportional to q^ whereas its height 
scales as l/|g|, the limit q reproduces the linear 
LL resuh, S{q,uj) = 2Kc\q\S{uj - v^q). The DSF (135) 
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already satisfies the /-sum rule (Nozieres, 1997). Hence, 
additional features may exist with weights at most of 
order (q/kp)^- 

A second peak in S{q,oj) occurs at energies close to 
the spinon mass shell. Indeed, the coupling between 
spinous and holons has the remarkable consequence that 
the lower edge of support of the charge DSF now coin- 
cides with the shifted spinon spectrum. The weight of 
this additional peak can be estimated by using perturba- 
tion theory in the spin-charge coupling amplitudes k±, 
which can be also defined phenomenologically (Pereira 
and Sela, 2010), sec Eq. (199). 

For w ss VsQ, the total weight in the vicinity of w « VgQ 
equals Kc{a- + a+)^g^/12, where a± = k±/{vc ± Vs). 
At small q, the perturbation theory correctly predicts 
the peak with weight oc q"^ at the spinon mass shell. 
However, it is unable to predict the precise shape of the 
peak. An analysis using a mobile impurity Hamiltonian 
reveals again that at the lower threshold S(q, lu) has a 
power-law singularity (Pereira and Sela, 2010). The sin- 
gularity remains intact for arbitrary momenta (Schmidt 
et at, 2010a). The calculation of the exponent requires 
a projection of the charge density operator p{x). Using 
Eq. (118), one finds that the configuration with least 
energy for a momentum < q < 2kF reads p{x) — > 
e^'^^'i'RsPRsdsFiis and has the energy x \es{kF - q)\- In- 
terestingly, this configuration is identical to the threshold 
configuration for the spin structure factor (132). There- 
fore, near the edge of support S{q,oj) oc S^^{q,uj). In 
particular, the threshold exponent for S{q, oj) is given by 
Eq. (134). 



G. Finite-size and finite-temperature efTects 

One of the successes of the linear LL theory is its ability 
to easily predict finite-size and finite-temperature effects. 
This can be achieved because the Gaussian Hamiltonian 
of the LL is conformally invariant (Gogolin et al., 1998), 
which results in a universality of the finite-size correc- 
tions, such as in the oc 1/L correction to the ground 
state energy or the oc T correction to specific heat at 
low temperatures (Aflileck, 1986; Blote et al, 1986). It 
can be used as a powerful tool when combined with exact 
solutions, since the energies in the latter case can be eval- 
uated up to 1/L corrections: such an approach was used 
to evaluate the Luttingcr parameters of the Lieb-Liniger, 
the XXZ (Bogoliubov et al., 1987) and the ID Hubbard 
models (Prahm and Korepin, 1990). Conformal invari- 
ance also fixes the time dependence and the finite-size 
effects in correlation functions. 

For example, for spinless fermions at p|x±t!t| ^ 1, one 



has (Cazahlla, 2004) 

{^\0,0)^ix,t))^^. 



(136) 
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i{2n+l)kipx 



2z(-l)" {ip{vt + x) + Of"- {ip{vt -x) + Of ' 

where C„ are "nonunivcrsal" prefactors, and 

= (2^^ + l)^i^/4 ± (2n + l)/2 + 1/AK > 0. 

For a finite system with periodic boundary condi- 
tions on a circle of length L, conformal invariance dic- 
tates (Cazahlla, 2004; Shashi et al, 2011) that the nth 
term in Eq. (136) gets modified to 



n 



jT^(vt±x)/L \ ^^<-^' 



\*pisin^^±^+0^ 



• (137) 



The goal of this section is to promote the phenomeno- 
logical theory based on impurity Hamiltonians to the 
same status. It will not only provide new predictions, 
but will also serve as a calculational tool to extract in- 
formation from exactly solvable models (see Sec. Ill) and 
interpret the results of numerical simulations. 

Let us start with the discussion of finite-size effects. 
For concreteness, we will focus on spinless fermions and 
the vicinity of the edge of support. The finite-size spec- 
trum of oc 1/L corrections is that of a shifted Gaussian 
conformal theory (Tsukamoto et al., 1998), and can be 
determined using conventional techniques. The correc- 
tion to the position of the edge in standard notations (Ko- 
repin et al, 1993) is given by (Pereira et al, 2008, 2009) 



2'KV 
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4^ 



(AAT- 



Hmp 



y + K{D - dimp) 



(138) 

where AA'' and D are quantum numbers specifying the 
excitations, while riimp and dimp are related to phe- 
nomenological phases in Eq. (75) via 

6.{k)^S+ik) ^ 6+{k) 



2'k/VK 



5.{k) 



4v^7r 



(139) 



The quantum numbers AN and D are related to 
the numbers and Nj^ of fermions created at each 
Fermi point, which define the exponents as discussed 
in Sec. II. C. For fermionic models, they are given by 
AN = Nr + Nl and D = {Nr- Nl)/2. Thus, Eq. (138) 
allows one to calculate the phase-shifts S±{k) by, e.g., 
numerical tracking of the oc 1/L corrections to energies. 

We will now consider the finite-size oc 1/L structure 
of low- lying levels at fixed k, as well as the scalings 
of various matrix elements (Shashi et al, 2011). Us- 
ing a resolution of the identity in the expectation value 
(*^(0,0)*(a;,t)), we get 

(0, 0)*(a;, t)) = J2 e-*('=^'^-^»'*'|(s', N - l\^\N)f, 

(140) 
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where \s' , N —1) denote eigenstates with TV — 1 particles. 
In a finite-size system, the formfactors in this equation 
have to be matched with the field-theoretical predictions. 

To understand the procedure, let us first consider the 
scaling of the formfactors between low-energy states for 
ID fermions (Bogoliubov et al., 1987), which can be de- 
scribed by the conventional LL theory. We can now ex- 
pand terms in Eq. (137) using a Fourier series as 



r(^ + n±) 



r(/x)r(n± + i)- 



(141) 



In this equation, the summation only over n+ is implied 
for right branch contribution, and the summation only 
over n_ is implied for left branch contribution. Plug- 
ging Eq. (141) into Eq. (137), and comparing the result 
to Eq. (140), one can clearly identify contributions from 
excitations at the right (left) Fermi branches with ener- 
gies 2-Kvn±/L > 0, respectively. To accommodate the 
additional momentum — (2n -|- l)kF, one needs to put 
n+1 holes at the right Fermi point and n particles at the 
left Fermi point. The contributions from n_|_ = n_ = 
give the scalings of the "parent" formfactors (Bogoliubov 
et al, 1987) 



{n,N -1\^/\N)\^ 



2(-l)" \poL 



(2n + l)^ K'^+l 
2K 



(142) 



where \n,N — 1) is the lowest energy state of N — 1 
fermions with momentum — (2n + l)fci?. The nontrivial 
scaling of this formfactor with L is a consequence of the 
criticality of the LL. The studies of scalings of the form- 
factors serve as a tool to evaluate the nonunivcrsal prefac- 
tors Cn (Shashi et ai, 2011, 2010), which are usually not 
known except for a few cases (Astrakharchik et al, 2006; 
Gangardt, 2004; Gangardt and Kamenev, 2001; Jimbo 
et al, 1980; Lukyanov and Terras, 2003; Popov, 1980; 
Vaidya and Tracy, 1979). Within oc 1/L accuracy, for 
n± > 2 the excited states of TV — 1 particles are degen- 
erate, while the degeneracy within each "multiplct" is 
lifted by oc corrections due to the nonlinear spec- 

trum. The imivcrsal Hamiltonian of Sec. II. B predicts 
the distribution of the spectral weight within each "mul- 
tiplet", as has been shown by Shashi et al. (2011). 

Let us now apply a similar logic to the finite-size be- 
havior of the response functions near the edge of support, 
and for concreteness we will focus on the spectral func- 
tion A{k,e) for < kp and e < 0. In addition to the 
LL, we now also need to take into account the finite-size 
quantization of the momentum of the impurity moving 
with velocity \vd\ < v. For an infinite-size system, A{k,e) 
in the vicinity of the edge eth{k) < (see notations in 



Sec. II. C) can be written as 

A{k,e) = Ao_{k) j dxdte-'^^*D{x,t)L{x,t)R{x,t), 

(143) 

where Se = e — eth{k), D{x, t) = 6{x — v^t) is the impu- 
rity correlator, L{R){x,t) = {i{vt±x)+0)~^">'--^('i\ and 
^o.-(fc) is a "nonunivcrsal" prefactor, see Eqs. (80) and 
(81). After {x,t) integration, Eq. (143) results in 



A{k,s) 



2Tr0{-Se)Ao^-{k)\Ss\-t'°-- 



r(l - no-){v + Vd)''°---'^{v - Vd) 



(144) 



In finite-size systems, L{x,t) and R{x,t) get modified, 
see Eq. (141). Similarly, the change of D{x,t) to 
Snc 6^""°^'^""''*^/^ corresponds to the quantization of 
the impurity momentum. Since we are considering fixed 
fc, the total momentum of the excitations on the left and 
right branches should be equal to the inverse of the shift 
of the momentum of the impurity 2'jTnD/L, which implies 
no = n_ — n+. Combining these terms, we get for A{k, s) 



sUe- AE 



27rn+ , , , , , 

+ —i—{v-Vd) H ^[v + Vd) X 



n±>0 



27rn_ 



('27i-)^o.~.«+'^o.-.-£- + ^ 



ljfiO,-,R+I^O,-,L 



(145) 



At cx 1/L accuracy, the finite-size structure of the re- 
sponse function is given by the sum of two generically 
incommensurate frequency "ladders", with the relative 
spectral weights in each "multiplet" controlled by the 
phase shifts 6±{k). Equation (145) allows for an analyti- 
cal or numerical evaluation of Ao,_(fc) based on the scal- 
ing of the single formfactor with n+ = n_ = as 



\{k;N -1\^\N)\' 



^0, 



■ (fc) /27r 



,R+IJ'0,-,L 



where |fc;A^ — 1) in the eigenstate of iV — 1 particles 
corresponding to the edge of support at < kp- The 
structure described by Eqs. (142) and (145) can be ex- 
plicitly confirmed for certain integrable models (Kitanine 
et al., 2009a,b, 2011; Shashi et al., 2010) using known ex- 
pressions for the finite-size formfactors. This provides a 
stringent microscopic check of the phenomenological im- 
purity Hamiltonians and allows to analytically calculate 
varioiis "nonunivcrsal" prefactors for these models. In 
addition, it allows to calculate various prefactors pertur- 
batively (Shashi et al., 2011). 

Let us now comment on the effects of finite tempera- 
tures, which are quite different for nonlinear LLs com- 
pared to linear ones. In the latter, conformal invariance 
allows one to calculate (Giamarchi, 2004) most of the 
finite-temperature effects by simple substitutions such as 
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w/L — > iT. This cannot be done for nonlinear LLs. The 
full analysis of the finite-temperature effects, especially 
in kinetic problems (see Sec. IV), remains an open prob- 
lem. Let us make some general remarks here focusing on 
the response functions of spinless fermions. 

For the interacting case, a finite temperature smears 
the sharp edges of support of the spectral function. How- 
ever, far away from the Fermi points and at |eth(fc)| ^ T 
(we use fcs = 1 throughout the text), there is a large 
interval of energies where the effect of temperature can 
be captured by substituting the Fermi point correlators 
by their finite-temperature versions (Karimi and Affleck, 
2011), i.e., in Eq. (143) 



LT{RT){x,t) 



sin [2mT{t±x/v)+Q]^ 



At the same time, the impurity correlator D(x, t) can 
be kept as a delta-function, because the correction to 
the impurity "occupation number" is exponentially sup- 
pressed. These substitutions result in imiversal functions 
characterizing the smearing of the edge singularities 



AT{k,e) = A, 



o,-(fc) J ' 



dte-'^'*LT{vdt,t)RT{vdt,t). 



Note that these functions can be evaluated numerically 
and generically have a strongly non-Lorentzian shape. 
The temperature mostly affects the response functions at 
energies of the order of ~ T from the edges of support. 

In the vicinities of the Fermi points, one can use the 
universal Hamiltonian of Sec. II. B to take a finite temper- 
ature into account. A naive extension of the above argu- 
ment would imply a smearing of the nonlinear LL physics 
for the DSF at temperatures on the order of ^ q^/{fh). 
We note, however, that this is not the case, as can be il- 
lustrated by the DSF of noninteracting fermions. The lat- 
ter can be straightforwardly evaluated as in Sec. II.A.l, 
and in the limit q 'ti kp reads 



So{q,co) 



n 



vf- 



±m{co — Vpq) — (Z^/S 



where the Fcrmi-Dirac distribution function np{e) = 
l/[exp (e/T) -|- 1] replaces the step functions in Eq. (10). 
One sees that the zero-temperature result (10) survives 
up to temperatures on the order of ^ vpq. Similarly, the 
universal Hamiltonian of Sec. II. B implies that at finite 
temperatures Eq. (54) is replaced by 



9 4- 



±m{ui - vq) - /2 



(147) 



so that the DSF is barely different from its T = form 
as long as the temperature is small compared to vq. The 
mechanism of the DSF smearing by finite temperature 
expressed in Eq. (146) is the same as in Eq. (147). At 



q kp, the DSF involves only contributions from the 
right-moving quasiparticle and quasihole with momenta 
of the order ~ q around the Fermi point. A finite tem- 
perature smears their velocities by 5v ^ T/{rhv). The 
corresponding energy variation q5v ~ qT/[fhv) is small 
compared to q^ /rh as long as T <C vq. 

There is a substantial difference between the domains 
of apphcability of Eqs. (146) and Eq. (147) though. In 
the former, the limits g — )■ and T — )■ may be taken in 
any order. The latter is valid only in the limit when, in 
addition to keeping ifi{Lo~vq)/q^ constant at g — >■ as in 
Eqs. (54) and (63), T/{qv) is kept constant. As a result, 
Eq. (147) does not hold at fixed T and q ^ 0, which will 
be important in Sec. IV. B. 

Unlike the DSF, the spectral function in the univer- 
sal limit is a convolution of contributions from both left 
and right Fermi points. The kinematic considerations of 
Sec. II. A and Sec. II. B imply that the quasiparticles at 
the left branch have energies of the order ~ {k — kp)^/m,, 
so the finite-temperature smearing of the nonlinear ef- 
fects in the spectral function is a two-step process (Ma 
and Imambekov, 2012). 

First, at temperatures ~ {k—kp)^/m, the contribution 
from the left branch gets smeared out. The contribution 
from the right branch gets affected significantly only at 
temperatures on the order of ~ v{k — kp). It should be 
noted that the spectral function of chiral fermions at the 
edges of quantum Hall states (Altimiras et ai, 2009; Al- 
timiras et al, 2010; Chang, 2003; Granger et a/., 2009 
Heyl et al., 2010; Jolad et ai, 2010; Limdc et al, 2010 
Neuenhahn and Marquardt, 2009; Paradiso et al., 2011 
le Sueur et al., 2010) should be more robust to finite tem- 
peratures due to the absence of the counter-propagating 
branch. 



, (146) Real-space correlation functions 



In the previous sections, we showed that the power- 
law singularities of the spectral function A{k, e) near the 
edge of support can be described with the help of mobile- 
impurity Hamiltonians for ID fermionic, bosonic and spin 
systems. We will show in this section that the existence 
of these threshold singularities has important implica- 
tions for the space-time correlation functions in the limit 
of large x and t. In particular, we will elucidate the 
connection between the threshold singularities and the 
breakdown of conformal invariance, focusing on the case 
of spinless fermions. 

Using a Lehmann spectral representation (Abrikosov 

et ai, 1963), various space-time Green's functions can be 
obtained by Fourier transforming the spectral function. 
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n = 1 J- 







(b) 



FIG. 10 (Color online) (a) Calculation of the real-space cor- 
relation functions. The integration range for the calculation 
of (*''^(0, 0)*(a;,t)) from A{k,e) is shaded. For large t and 
fixed xjt — Vd, the contributions to the integral come from 
a region around e = 0, and from points at which the lines 
e — kvd ~ const, touch the edge of support, (b) Breakdown of 
conformal invariance in (vp^ (0, 0)^(a; = vat, i)) . In the shaded 
areas {\vd\ > v) only the linear Luttinger liquid (LL) power 
laws survive, see Eq. (136). In the white areas, for \vd\ < v, 
new power laws appear in addition, see Eq. (151). 



For example, 

(¥{0,0)^{x,t)) = ^ I dk r d£e''="e-*^*A(fc,e). 

(148) 

In a similar way, (^'(a;, t)^^(0, 0)) can be expressed via 
the Fourier transform of A{k, e > 0) in the particle sector. 

As is well known in the theory of Fourier transfor- 
mations (Bleistein and Handelsman, 1986), the non- 
analyticities of A{k, e) control the long space-time be- 
havior of its Fourier transforms. The spectral function 
is non-analytic in the vicinities of the Fermi points and 
their 2nkp images. We showed in the previous sections 
that the regions where A{k, e) deviates significantly from 
the predictions of the linear LL theory become narrow 
in the limit e — )■ 0. As a consequence, the effects of the 
spectrum nonlinearity are suppressed near e ~ when in- 
tegrating over k. Neglecting the nonlinear effects in the 
vicinities of the Fermi points produces the well-known 
space-time power-law behavior of the correlation func- 
tions in Eq.(136) at p\x ± vt\ 3> 1. This result is man- 
ifestly conformal invariant due to the conformal invari- 
ance of the LL Hamiltonian (41). If one considers the 
limit 



t ^ oo and Vd = x/t fixed. 



(149) 



Eq. (136) results in the same set of power-law tails in t 
irrespective of the ratio Vd/v. 

We will show now that in the limit (149), the threshold 
singularities in A{k, e) may generate a new set of power- 
law tails in t. For generic interacting spinless fermions. 



such power laws appear only for \vd\ < v and, moreover, 
the exponents depend on the ratio Vd/v. The origin of 
these power laws is very similar to that of van Hove sin- 
gularities in the density of states of noninteracting sys- 
tems (Landau and Lifshitz, 1980) and is illustrated in 
Fig. 10a. To study the limit (149), it is convenient first 
to make a rotation of coordinates in the [x, t) and {k, e) 
planes to 



X - Vdt 



k' = £/vd + fc. 



t' 



X + Vdt 



2vd 

e' = £ — kvd- 



(150) 



In these coordinates, one has kx — et = k'x' — e't\ and 
we are interested in the limit x' — Q, t' = t oo. There- 
fore, to study this limit one needs first to keep e' constant 
(tangential lines in Fig. 10a) and integrate over k' . The 
non-analyticities near the Fermi points and their 2nkp 
images produce the conventional LL power laws (136). 
In addition, after the integral over k' is performed, for 
\vd\ < V, the existence of sharp edges of support pro- 
duces new non-analyticities in e' from the vicinities of 
the touching points, see Fig. 10a. The origin of these 
non-analyticities is very similar to that of van Hove sin- 
gularities, but unlike the latter they also depend on the 
threshold exponents of the spectral function. The touch- 
ing condition implies that Vd is nothing but the impurity 
velocity. Each "shadow band" produces a separate power 
law and an evaluation of the Fourier transform (148) re- 
sults in 

{^\0,0)^{x - Vdt,t)) = (*t(o,0)*(i;di,<))LL (151) 



tl/2(t;t + x)^">->-f' {Vt - XY^.-.R 



where the correlation function of the linear LL 
(*'l'(0, 0)*('Udi, i))LL defined in Eq. (136) and the mo- 
mentum kd is defined by the touching condition Vd — 

^^^qI:^^ ■ The effective mass is given by l/rud — 

k=kd 

dvd/dk\k=kd- The exponents fJ,n,-,R{L) are defined by 
Eq. (82) and An.-{kd) are the non-universal prefactors 
from Sec. H.G. 

We note that for sufficiently weak interactions the new 
"nonlinear" tails in t decay slower than the linear LL tails 
for all \vd\ < V. Indeed, for noninteracting fermions and 
Vd — 0, the contributions from the Fermi points decay 
as cx 1/t, whereas the contribution from the bottom of 
the band decays only as cx due to a conventional 

van Hove singularity there (Gutman, 2008). For the in- 
tegrable Lieb-Liniger model, expansions similar to our 
Eq. (151) have been obtained from purely microscopic 
considerations (Kozlowski, 2011; Kozlowski and Terras, 
2011) and match the exact results for the exponents cal- 
culated in Sec. III.B. Similar results have been obtained 
for a gas of ID Lieb-Liniger anyons (Patu et ai, 2009). 
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In Sec. II. B, we have cstabUshcd that in the vicinity 
of the Fermi points, the response functions are univer- 
sal within the nonUnear LL theory. The results of this 
section then imply that for |i;d±t;| ^ v, the time depen- 
dence of the field correlator, and its crossover between 
linear and nonlinear regimes is universal as well. Us- 
ing the universal expressions (48) for the phase shifts 
in this limit and Eqs. (136) and (151), one can estab- 
lish that /Ul + M-R > I + A*o,-,i? + Mo,-,l- Therefore, the 
nonlinear results always decays slower in time than the 
linear LL result. Using a more careful analysis of the 
prefactors (Imambekov and Glazman, 2009b), one can 
also establish the scaling of the crossover time tc be- 

4-G'/K + 4K 

twecn the two regimes as pvtc ~ ^ "^^"^'^•' ^ 2~^vk+2k ^ 
For weakly interacting fermions, this condition reduces 

to pVptc ^ Vp/ {vp - Vd)"^- 

Let us now comment on the space and time be- 
havior of the transverse spin correlation function of 
SU{2) invariant spinful bosons. In Sec. II. E, it has 
been established that for small momenta (9,w) a 

(a;-g2/(2m*))i^-^«'/(27r"p2)_ Qf particular interest is the 
regime when repulsive interactions between spin-up and 
spin-down particles are strong. Then the spin-down 
impurity cannot exchange positions with other parti- 
cles and effectively becomes trapped (Zvonarev et al., 
2007). In the thermodynamic limit, the magnon mass 
TO* diverges and the bandwidth of spin excitations be- 
comes very narrow. This is very reminiscent of the 
narrow-band spinon excitations of strongly-interacting 
s = 1/2 fermions discussed by Matveev (2004a) ; Matveev 
et al. (2007a, b), and mentioned in Sec. II. F. In both 
cases, the large difference between the bandwidths of 
spin and charge excitations results in the existence of 
a new regime where a certain new universal behav- 
ior of correlations takes place. Performing the inverse 
of the Fourier transformation (113) within the saddle 
point approximation, one obtains (Zvonarev et al., 2007) 



(S'+(.T,f)5-(0,0)) cx 



K 



In {Ept) 



it 
2m* 



-1/2 



exp ■ 



is introduced to provide a short-time cutoff. For generic 

interactions, the logarithmic term can be ignored and one 
obtains the scaling oc t/m* characteristic of a single- 
particle wave packet spreading. However, in the limit of 
infinitely strong interactions to* diverges and one obtains 
a logarithmic scaling (npx)'^ ~ Kln^Ept). For large but 
finite to* a logarithmic scaling law is applicable in an in- 
termediate time interval, the length of which grows with 
the increase of to*. 



LL description. However, in ID there exists a class of 
exactly solvable (or integrable) models, for which energy 
spectra and thermodynamical properties can be calcu- 
lated exactly using the Bcthc ansatz. The calculation of 
their correlation functions, on the other hand, is a much 
more complicated task (Korepin et al., 1993), because 
it requires not only the knowledge of the exact energies, 
but also of the matrix elements (formfactors) . In addi- 
tion, one needs to be able to sum over all excited stated 
in the Lehmann representation in order to obtain an- 
swers in the thermodynamic limit. The exact expressions 
for the formfactors are usually known only in finite-size 
systems (Ha, 1996; Kitanine et al., 1999; Kojima et al., 
1997; Slavnov, 1989, 1990), and so far only few-spinon 
approximations (Bougourzi, 1996; Bougourzi et al., 1998; 
Karbach et al., 1997) or fully numerical summations over 
formfactors (Caux and Calabrese, 2006; Caux et al, 2007, 
2005; Caux and Maillet, 2005; Gritsev et al, 2010; Kohno 
et al, 2010) have been implemented for certain gapless 
models. 

In this section, we review some recent results for the 
correlation functions of integrable models which have 
been obtained by combining exact results with field the- 
ories beyond the linear LL theory. We will show that 
a plethora of new results can be derived using this ap- 
proach. More importantly, exactly solvable models also 
provide nontrivial checks of the phenomenology and pro- 
vide additional verification of the effective impurity mod- 
els. It should be noted that, historically, the analysis of 
exactly solvable models (Haldane, 1980, 1981a) played an 
important role in the justification of the linear LL theory. 
Now these models are proving their worth to extensions 
of the LL theory as well. 

The rest of this section is organized as follows. 
In Sec. III. A we review some recent progress for 
the Calogero-Sutherland model (Calogero, 1969, 1971; 
Sutherland, 1971, 2004). In Sec. III.B we discuss 
the Lieb-Liniger model and demonstrate several ap- 
proaches (Cheianov and Pustilnik, 2008; Imambekov and 
Glazman, 2008; Pereira et al, 2008) to extract parame- 
ters of the effective impurity Hamiltonians. In Sec. III.C 
we consider Yang-Gaudin models (Gaudin, 1967, 1983; 
Yang, 1967) which describe spinful multi-component sys- 
tems interacting via a contact potential. In Sec. III.D, 
we discuss ID lattice models, such as the spin- 1/2 XXZ 
model (which is equivalent to spinless fermions on a lat- 
tice), and ID Hubbard model (Essler et al, 2005; Lieb 
and Wu, 1968). 



A. Calogero-Sutherland model 



III. EXACTLY SOLVABLE MODELS 

In the previous section, we concentrated on the prop- 
erties of generic ID quantum systems beyond the linear 



In this subsection, we will discuss the Calogero- 
Sutherland (CS) model (Calogero, 1969, 1971; Suther- 
land, 1971, 2004), the most well studied model with in- 
verse square interactions (Haldane, 1988; Kuramoto and 
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FIG. 11 (Color online) (a) The DSF S{q,uj) of the Calogero- 
Sutherland (CS) model differs from zero only in a finite in- 
terval of frequencies tj_ < a; < At the boundaries of this 
interval, S{q,uj) exhibits power-law singularities, see Eq. (156) 
(b) Dependence of S{q,uj) on a; at a fixed q < 2kF and for 
repulsive interactions, A > 1. 



Kato, 2009; Kuramoto and Yokoyama, 1991; Shastry, 
1988). Such models are special among exactly solvable 
models, because their ground state wavefunctions can of- 
ten be expressed as Jastrow-type products, and their dy- 
namical correlation functions can be derived in a closed 
form as multiple integrals because of special properties 
of their excitations (Arikawa and Saiga, 2006; Arikawa 
et al, 2001, 1999, 2004; Ha, 1994, 1995; Ha, 1996; Ha 
and Haldane, 1994; Haldane, 1991; Haldane et al, 1992; 
Haldane and Zirnbauer, 1993; Kato, 1998; Lesage et al., 
1995; Pustilnik, 2006; Talstra and Haldane, 1994; Ya- 
mamoto and Arikawa, 1999; Yamamoto et al., 2000). The 
special quantum mechanical properties of a system with 
inverse-square (cx 1/r^) interaction potential can already 
be seen at the classical level, because a liquid of such 
particles admits a description using integrable hydrody- 
namics (Abanov et al., 2009; Abanov and Wiegmann, 
2005; Kulkarni et al, 2009; Stone et al, 2008). For con- 
creteness, here we will focus on the DSF of the CS model 
following Khodas et al (2007b) and Pustilnik (2006), and 
illustrate connections to the phenomenology discussed in 
Sec. H.C. Most of the related results on other inverse- 
square models are reviewed in a recent book by Kuramoto 
and Kato (2009). 

The CS Hamiltonian is given by 



1 



N 



92 



CS 



2m ^ dx-f 

i—l ^ 



^y(j:, - Xj). 



(152) 



Here m is the mass, and the interaction potential is 
A(A - l)/m 



Vix) 



A(A-l) ^ . 



{L/nY siv^ly-Kx/ L) mx^ 



where A > 1/2 is a dimensionless interaction strength 
defining the LL parameter as K = 1/A. 

Excitations can be simply classified using the language 
of quasiparticles and quasiholes with respect to a filled 
"Fermi sea" . The long range of the interactions results 
in a rather peculiar excitation spectrum. For quasiparti- 
cles with velocities u+ larger in absolute value than the 



sound velocity v = nXp/m, the spectrum is given by 
to(u^ — v'^)/2. In contrast, for quasiholes with veloci- 
ties < V, the spectrum is given by mX{v'^ — wi)/2. 
The discontinuity of the effective mass near the Fermi 
points can be expected from perturbation theory, since 
the Fourier transform of the l/r^ potential is nonana- 
lytic, Vk oc Because of that, the universal Hamilto- 
nian of Sec. II. B is not applicable to the CS model. 

A special feature of this model is that for rational 
X — r/s (where r and s are coprime), the operator Pg^Q 
acting on the ground state creates only right-moving ex- 
citations: s quasiparticles and r quasiholes (Ha, 1994, 
1995; Ha, 1996). This is not a generic property of inte- 
grable systems and has a profound effect on the dynamic 
response functions. In particular, since no left-moving 
excitations are created, the DSF S{q,uj) is nonzero only 
in a finite interval of energies, w_(q) < lu < Lo+{q), see 
Fig. 11. The existence of an upper threshold is not ex- 
pected for generic ID systems, see Sec. II.A.l. The upper 
threshold corresponds to a configuration where the entire 
energy is given to a single quasiparticle, while for q < 2kp 
the lower threshold corresponds to all energy being given 
to a single quasihole. One can show that for q > 0, 



uj+{q) = vq + q^/{2m), 

(q) ^ vq — Xq^ / I2m) for q < 2kp. 



(153) 
(154) 



The DSF can be written as 

S{q,uj) = q' j \{dv+,^dv_^,F,^r5{q-P)5{io-E), 

(155) 

where P and E are the total momentum and energy of 
the excitations, respectively, and the expression for the 
formfactor F^^r is given by (Ha, 1994, 1995; Ha, 1996; 
Haldane, 1995) 



1-1/A 



The analysis of the multidimensional integral in Eq. (155) 
performed by Pustilnik (2006) then yields a power-law 
behavior in the allowed regions as 



S{q,uj) 
m/q 



for \uj — uj± I ^ ujj^ — a;_ . 

(156) 



The prefactors of the DSF have also been evaluated fol- 
lowing Sec. II.G (Shashi et al, 2010). 

The analysis of the spectral function can be performed 
similarly (Khodas et al, 2007b), and for |fc| < /cj;-, e < 
results in the power-law behavior 

A{k, e) cx 0(eth(fc) - e){e,^{k) - e)i-(^-i)V(2A)^ (^57) 
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These exponents at generic edges of support ±eth{k) = 
=Fw_(A;f — k) can be simply recovered from the phe- 
nomenological considerations of Sec. II. C, with Eq. (75) 
resulting in momentum-independent phase shifts 

gcs ^cs i/A-l\ lfK-l\ , , 

We see that near the right Fermi point, the phase shift 
follows the prediction of Eq. (48) in Sec. II. B, while 
8'^^ differs from the universal prediction, as expected due 
to the slow decay of the inverse-square potential. 



B. Lieb-Liniger model 

Arguably the simplest exactly solvable model, the 
seminal Lieb-Liniger (Korepin et al, 1993; Lieb, 1963; 
Licb and Linigcr, 1963) model of ID bosons interact- 
ing via a contact potential, played an important role in 
the development of both Bethe ansatz ideas (Korepin 
et al, 1993) as well as the LL description (Efetov and 
Larkin, 1975; Haldane, 1981a). Although the studies of 
this model were mostly an academic exercise for more 
than 40 years, its experimental realization in ultracold 
atomic gases (van Amerongen et al., 2008; Haller et al., 
2009; Kinoshita et al., 2004; Kinoshita et al., 2005; Ki- 
noshita et at, 2006; Tolra et at, 2004) now allows for 
a parameter-free comparison of theoretical predictions 
with measurements. Generally, ultracold ID Bose gases 
are realized by loading a Bose-Einstein condensate into 
a deep 2D optical lattice in the y — z plane formed by 
perpendicular laser beams, or using atom chips. The 
tight transversal confinement inhibits the occupation of 
higher transverse modes and provides a clean realization 
of the Lieb-Liniger model (Olshanii, 1998). These experi- 
ments stimulated significant interest in the long-standing 
problem of calculating its correlation functions (Kore- 
pin et al., 1993), which can be measured using inter- 
ference (Donncr et al, 2007; Hofferberth et al, 2007, 
2008; Imambekov et al., 2007, 2008; Polkovnikov et al., 
2006), photoassociation (Kinoshita et al., 2005), anal- 
ysis of particle losses (Haller et al, 2011; Tolra et al., 
2004), density fluctuation statistics (Armijo et al., 2010; 
Jacqmin et al., 2011), time-of-flight correlation statis- 
tics (Hodgman et at, 2011; Imambekov et al., 2009; Manz 
et al., 2010), scanning electron microscopy (Guarrera 
et al., 2012, 2011), or Bragg and photoemission spec- 
troscopy (Clement et al., 2009; Dao et al., 2007; Ernst 
et al, 2010; Fabbri et al, 2009; Gaebler et al, 2010; Papp 
et al, 2008; Stamper-Kurn et al, 1999; Stewart et al, 
2008; Vceravalli et al, 2008). Recently many new theo- 
retical results were obtained in this direction (Calabrese 
and Caux, 2007; Caux and Calabrese, 2006; Caux et al, 
2007; Cazalilla et al, 2011; Cheianov et al, 2006a,b; 
Cherny and Brand, 2009; Deuar et al, 2009; Gangardt 



and Shlyapnikov, 2003a, b; Golovach et al, 2009; Imam- 
bekov and Glazman, 2008; Kheruntsyan et al, 2003; Kho- 
das et al, 2008, 2007a; Kitanine et al, 2009a; Kormos 
et al, 2011, 2009, 2010; Kozlowski, 2011; Kozlowski and 
Terras, 2011; Pozsgay, 2011; Shashi et al, 2010; Sykes 
et al, 2008), but a fully analytical calculation of the cor- 
relation functions is still lacking. In this subsection, we 
will review recent progress for the Lieb-Liniger model 
based on combining the phenomenology beyond the LL 
theory with the Bethe ansatz. 

The exactly solvable Lieb-Liniger model is given by 

I ^ q2 

3=1 i l<j<k<N 

where c > is the interaction constant and m is the par- 
ticle mass. In the thermodynamic limit, the ground state 
is fully characterized by the dimensionless parameter 

7 = 2mc/p, (160) 

where p = N/L is the density. The regime of weak in- 
teractions corresponds to 7 ^ 1, while strong repulsion, 
i.e., the Tonks-Girardeau limit (Girardeau, 1960) corre- 
sponds to 7 ^ 1. The LL parameter equals K = vp/v, 
where v is the sound velocity and vp = Trp/mis the Fermi 
velocity of a noninteracting Fermi gas of density p. The 
parameter K is uniquely defined by 7, with K « tt^~^^^ 
for 7 < 1 and « 1 -I- 4/7 for 7 > 1 (Cazalilla, 2004). 

Let us briefly review the solution of the Lieb-Liniger 
model to introduce the notation. We will mostly follow 
the conventions of Korepin et al. (1993). The ground 
state quasimomenta i/j (1 < j < N) are given by the 
solutions of nonlinear Bethe equations 

N 

Lvj + 9{i'j — Vk) — iwrij, (161) 
fe=i 

where 9{i') = 2 arctan 5^ is the two-particle phase shift 

and the ground state quantum numbers are rij = j — 1 — 
[N — l)/2. In the thermodynamic limit, this system gives 
rise to the linear integral equation 

Pi^) - ^ /"^ Kiu, n)pir,)dv = ^. (162) 

Here, p{i^) = limL^oo ^/{L{i^k+i — i^k)) is the density of 
roots, K{i/, rj) ~ Amc/[{2mc)'^ + {v — riY], and Q {-Q) is 
the highest (lowest) filled quasimomentum; Q is defined 
by the normalization condition p — J^g p(y)dv. 

Particle-like excitations (Lieb-I mode) can be con- 
structed by adding an extra quasimomentum |A| > Q, 
while hole-like excitations (Lieb-II mode) are obtained 
by removing a quasimomentum |A| < Q. Such excita- 
tions change the total number of particles, and it is c\is- 
tomary in the literature (Korepin et al, 1993) to change 
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FIG. 12 (Color online) (a) Dynamic structure factor (DSF) 
S{q,uj) and (b) spectral function A{q,e) for the Lieb-Liniger 
model. Shaded areas indicate the regions where the functions 
are nonvanishing. Lieb's particle mode £i{q) and hole excita- 
tion mode £2(9) are indicated. 



the boundary conditions for wavefunction from periodic 
to antiperiodic when the number of particles changes by 
±1. Since all quasimomenta are coupled to each other 
by Eq. (161), the addition of an extra particle or hole 
will shift all quasimomenta. A convenient way to take 
this change into account is to introduce a shift function 



1+1 



(163) 



where Vj are the new solutions with antiperiodic bound- 
ary conditions, and the upper (lower) sign corresponds 
to an extra particle (hole). In the thermodynamic limit, 
F(i>\X) satisfies an integral equation 



F(i.|A) 



1 

2n 



K{v, r])F{rj\X)dri 



-Q 



2tt ' 



(164) 



Due to the antiperiodic boundary conditions of i'j, the 
shift function F{i'\X) in fact corresponds to the fermionic 
Cheon-Shigehara model (Cheon and Shigehara, 1998, 
1999), which is dual to the Lieb-Liniger model. However, 
all results for the Lieb-Liniger model can be formulated 
using Jordan- Wigner strings and F{i'\X), so we will use 
the fermionic language for consistency with Sec. II. C. 

As will be shown below, the shift functions F{zLQ\X) 
play a crucial role in the calculation of the edge singu- 
larities, so we will investigate them in more detail. One 
can analytically derive the limiting behavior 



FiQlQ) = 1 
F{-Q\Q) = 



1 VK 
27k 
1 ^/K 

27W ' 



(165) 
(166) 



and F{zLQ\ — Q) = —F{'^Q\Q). Moreover, one can show, 
VK 2mcVK 



F{±Q\X) 



ttA 



for Q, mc < A. (167) 



Equations (165)-(166) have been derived by Korepin and 
Slavnov (1998), and Eq. (167) follows from an expansion 



of the right hand side of Eq. (164) combined with 



p{±Q) = VX/27r, 



(168) 



see, e.g., Eqs. (I.9.20)-(I.9.22) in Korepin et al. (1993). 

The shift function can be used to calculate the exact 
energies of Lieb's particle (ei > 0) and hole (£2 > 0) 
excitations as a function of the momentum q{X). They 
are given by ei,2(<z) — ±e(A), with e(A) defined by 



e(A)- 



1 

2tt 



K{X, ri)e{v)d'n = AV(2to) - /i, (169) 



where /i is the chemical potential and e(±Q) = 0. 
The momentum corresponding to a quasimomentum A 
is given by 



q{X) =± \ X-T:p+ e{X~ v)p{v)dv 



Q 



(170) 



Here the upper (lower) sign corresponds to a particle 
(hole) excitation with X > Q (|A| < Q), and q{Q) — 0, 
q{-Q + 0) = 27rp = 2kF. Equations (169) and (170), to- 
gether with the normalization condition mentioned ear- 
lier for p^v) provide the full set of equations to determine 
the form of £1,2(9): see Fig. 12. Lieb's particle and hole 
modes can be simply understood in the limits of weak 
and strong interactions. 

In the limit of weak interactions (7 ^ 1) ID bosons 
form a quasicondensate characterized by slow algebraic 
decay of real-space correlation functions at long distances 
(Mora and Castin, 2003; Petrov et ai, 2000; Popov, 
1980). On a semiclassical level, its state can be described 
by a macroscopic wave function 'i'{x,t) which is a solu- 
tion of the ID Gross-Pitaevskii equation (Pitaevskii and 
Stringari, 2003), 

1 



idt'^{x,t) + £-9^*(a;,t) + 2c(p - |*(x, t)n*(x, i) = 0. 

(171) 



2m 



One class of wave-like solutions of Eq. (171) describes Bo- 
goliubov quasiparticles with a dispersion relation £i{q) — 
w(7-\/T^P7W2"^J^5 where v — y/^p/m is the sound veloc- 
ity for the weakly interacting gas. A second class of solu- 
tions of Eq. (171) are dark solitons (Khodas et ai, 2008; 
Kulish et ai, 1976) which have energies £2(9) < £i{q)- 
They correspond to localized perturbations of the qua- 
sicondensate density and travel at velocities Vs(q) = 
d£2{q)/dq < V. The spectrum is defined implicitly by 



£2 



ipv , 3 

— sm^ 



q 



sm 



.)] 



(172) 



where the parameter dg G [0, 2tt] is related to the ratio 
of the soliton velocity Vs{q) and the sound velocity v by 
cosies/2) = Vs/v. 
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In the Tonks-Girardeau Hmit (Girardeau, 1960) of 
strong interactions (7 ^ 1), a large contact repulsion 
enforces the equivalent of a Pauli exclusion principle, 
so the system becomes analogous to weakly interacting 
fermions. For < g < 27rp, the Lieb-II mode approaches 
£2(9) = vq — /{2m). The Lieb-I mode approaches 
£i{q) = vq + q^ /(2m)) for any q. The exact solution de- 
scribed above smoothly interpolates between the limits 
of strong and weak interactions. The small-g behavior 
of the spectra preserves the form £1,2(9) = vq± q^ /{2m) 
at any 7, with v and m being functions of 7. We find 
m = 47ri/2m/(37i/*) at 7 < 1 from Eq. (50). Matching 
the above asymptotes with the spectra of the correspond- 
ing wavelike solutions of Eq. (171) determines the region 
of applicability of the asymptotes, q < /97^/* (Khodas 
et al, 2008). 

Let us now describe the response functions, and for 
concreteness we will focus on the vicinity of £2(9)- As 
has been explained in Sec. II. G, the phase shifts 5±{k) 
define the exponents of the dynamic response functions 
and can be extracted from the Bethe ansatz by various 
techniques. The first approach, which docs not rely on 
integrability, is based on Eq. (75), where £th(fc = kp — 
q) = —£2(9) < 0. Another technique is based on the 
calculation of the finite-size (cx 1/L) energy shifts defined 
by Eqs. (138) and (139). This was pioneered by Pereira 
et al. (2008, 2009) and results in (see Eqs. (138)-(139) for 
definitions) 

rQ 

nimp = / P%mp{v)dv, (173) 

J-Q 

/-Q roo 
Pimp{j^)diy- / Pimp{y)dv, (174) 
-00 Jq 

where Pimp{'^) is defined by 

1 /•« 

Pimp{v)-^ / K{v,i^)pimp{'n)dr] = -K{v,X). (175) 

271" J-Q 

Finally, Ghcianov and Pustilnik (2008) and Imambekov 
and Glazman (2008) showed that 

^±(fc) =27rF(±Q|A). (176) 

This result is based on an identification of the operator 

U in Eq. (71) with a boundary condition changing op- 
erator (AfHeck and Ludwig, 1994; Schotte and Schotte, 
1969) of the LL particles whenever they pass the mobile 
impurity d. The LL describes the low energy particles 
near Fermi points ±Q, and thus the phase shifts 5±{k) are 
proportional to the shifts of the quasimomenta of parti- 
cles near these points. Since the shift functions F{±Q\X) 
are proportional to the shifts of the quasimomenta due to 
the presence of a hole, they are proportional to the phase 
shifts 5±{k). The proportionality coefficient in Eq. (176) 
is fixed by requiring that the excitation of a particle near 



±Q to the next allowed energy state corresponds to a 
phase shift ±27r. 

It is not at all obvious that three approaches described 
above should lead to the same phase shifts S±{k). It was 
shown analytically by Pereira et al. (2009) that the pre- 
dictions of the latter two coincide, and it can be checked 
numerically that Eqs. (75) result in the same phase shifts. 
Thus the coincidence of the three difi'erent predictions of 
the phenomenological Hamiltonian given by Eqs. (64)- 
(65) provides an unambigiious microscopic confirmation 
for its validity for the Lieb-Liniger model. Since the 
eff'ective field theory near the edge of support did not 
rely on integrability, it is natural to assume that such 
an approach holds universally for a large class of micro- 
scopic non-integrable models as well. Equations (165)- 
(166) together with (176) arc nothing but the universal 
phase shifts given by Eq. (48). Here the universal phase 
shifts were derived in a purely microscopic fashion, which 
provides an independent check of renormalization group 
(RG) arguments of Sec. II. B. 

The phase shifts (176) evaluated at |A| < Q provide 
nonperturbativc expressions for the exponents of the DSF 
S{q,uj) and the spectral function A{k,e) at their edges 
of support, see Eq. (112). It should be noted however, 
that the response functions of integrable systems might 
also have protected singularities (or non-analyticities) 
within a continuum. For the Lieb-Liniger model, the 
most prominent of these occur at £ = -iiei{q). In ad- 
dition, various weaker "shadow" singularities occur at 
£ = ±e\{±q — 2nkF). The exponents can be calculated 
as in Sec. II. G by introducing the effective mobile impu- 
rity with Vd = dei{q)/dq > v, and they are given by (see 
Fig. 12 for notations) 

M = 1 — PO,R — P'0,L, 

where Po,r{l)j ^'^id Mo q: ij(L) defined by Eqs. (87) 
and (112), and one needs to use the phase shifts (176) 
evaluated at A > Q. 

For singularities which occur within a continuum, one 
can also calculate the "shoulder ratios" of the weights 
right above and below the singular line. For instance, for 
the DSF one obtains (Pereira et al., 2009) 

S[q,£i{q) + 6e] ^ sin(7r/iox) 
Se^o S[q,ei (9) - Se] sin(7rMo,fl) ' 

If the exponent JI is negative (i.e., there is a non- 
analyticity and not a singularity), then Eq. (178) de- 
scribes only the non-analytic parts. Similar to Eq. (60), 
the shoiilder ratios are determined only by the phase 
shifts. The perturbative result (31) can be also inter- 
preted in terms of perturbative phase shifts (25). In 
addition to the phenomenologically determined shoul- 
der ratios, for the Lieb-Liniger model one can com- 
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bine the known expressions for the finite size form- 
factors (Slavnov, 1989, 1990) with the field theoretical 
predictions of Sec. II. G to obtain predictions for the 
non-universal prefactors of the edge singularities (Shashi 
et al, 2010). Such a finite-size analysis also provides 
a microscopic justification for the existence of the singu- 
larities within a continuous spectrum for the Lieb-Liniger 
model. 



C. Yang-Gaudin models 

It has been established by Yang (1967) and Gaudin 
(1967) that the Hamiltonian (159) of the previous sec- 
tion remains exactly solvable using the so-called nested 
Bethe ansatz, even if one does not require the wave func- 
tion to be symmetric with respect to permutations of 
Xi and Xj. One can impose either symmetry or anti- 
symmetry with respect to permutations of certain sub- 
sets of Xi. This leads to a family of exactly solvable 
models for ID multi-component systems, such as spin- 
1/2 (Gaudin, 1967; Yang, 1967) and SU{N) (Sutherland, 
1968) fermions, as well as Bose-Bose (Fuchs et al., 2005; 
Guan et a/., 2007; Li et al, 2003) and Bose- Fermi (Batch- 
elor et al., 2005; Frahm and Palacios, 2005; Imambekov 
and Demler, 2006a,b; Lai and Yang, 1971) mixtures. In 
this section, we will review recent exact results for re- 
pulsive (iso)spin-l/2 bosonic and fermionic Yang-Gaudin 
models, illustrating connections with the universal phe- 
nomenological description of Sec. II. E and Sec. II. F, re- 
spectively. 

The general approach is based on the calculation of 
finite-size corrections to the edge state energies, and their 
interpretation in terms of phase shifts. This procedure is 
aided by the known finite-size structure of the effective 
impurity theories described in Sec. II.G. 

Let us start from the discussion of bosons (Zvonarev 
et al., 2009b), which have a ferromagnetic ground 
state (Eisenberg and Lieb, 2002). As in Sec. II. E, we 
choose the magnetization to be pointing in -|-z direction. 
The microscopic wave functions in the ferromagnetic sec- 
tor coincide with those of the Lieb-Liniger model, so 
all dynamic response functions which do not involve the 
spin-down state, such as S{q,uj), S^^{q,u)), and the spec- 
tral function for spin-up particles A^{q,e), coincide with 
those of the Lieb-Liniger model. The response functions 
which involve one spin-down state, e.g., (9,w), will 
have singularities at the magnon spectrum u;,„(g). The 
states which contain one magnon are characterized by a 
set of quasimomcnta {vi , . . . , h'N, A} which satisfy the set 
of equations (Gaudin, 1983) 

N 

Lvj + ^ e{vj - Vk) = 2Tmj + e{2vj - 2A) + tt. (179) 
fc=i 
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FIG. 13 (Color online) The function a{q) defining the trans- 
verse spin structure exponent for isospin-1/2 bosonic Yang- 
Gaudin model, see Eq. (181), is plotted for different values 
of the dimensionless coupling constant 7. The values of the 
Luttinger parameter K are indicated for each curve and cor- 
respond in increasing order to 7 = 00, 1.65, 0.56, 0.238 and 
0.109 respectively. Adapted from Zvonarev et al. (2009b). 

The total momentum P and the energy E are 

N 1 ^ 

For q <C TTp, the magnon spectrum can be expanded as 
^m{q) ~ 9^/(2™*), and the expression for m/m* as a 
function of 7 can be obtained analytically from the exact 
solution (Fuchs et al., 2005). It has the asymptotic be- 
havior 1 - 2^/(37r) for 7 < 1, and 277^/(37) for 7 > 1. 
An analysis of the finite-size corrections to the energy 
of the magnon allows one to derive equations similar to 
Eqs. (164) and (176) which define the phase shifts for ar- 
bitrary interactions and momenta, see Zvonarev et al. 
(2009b) for more details. Similar to the Lieb-Liniger 
model, the phase shifts evaluated from the finite-size 
corrections coincide numerically with the phenomeno- 
logical predictions (Kamenev and Glazman, 2009). In 
Fig. 13, we present the exact results of Zvonarev et al. 
(2009b) for the transverse spin structure exponent /im 
after reparametrization 

which is chosen such that a{q) vanishes at g = 0, 2kp. 

Considerably more complicated is the case of spin- 1/2 
fermions (Essler, 2010), since for repulsive interactions 
the ground state is a singlet (Lieb and Mattis, 1962) and 
both spin and charge Fermi surfaces are present, as was 
discussed in Sec. II. F. The existence of two Fermi points 
is implicitly built into the structure of the Bethe ansatz 
solution, because instead of a single set of quasimomenta 
Vi as in Eq. (161), one needs to introduce the spin quasi- 
momenta Kj which live in an auxiliary spin space. In 
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a finite-size system, periodic boundary conditions lead 
to (Lee et al, 2011; Yang, 1967) 

Luj = 2nlj - ^ e{2uj - 2A„), (182) 



N 



J2 ^(2Aa - 2uj) = 27r J„ - ^ ^(A„ - A^). (183) 

3 = 1 /3=1 

where in the first equation j = 1, ...,N, and in the second 
equation a = 1,...,N^^, while /j, Jq are integer or half- 
intcgcr depending on the parities of TV, 7V-|-. The energies 
and momenta of the eigenstates are given by Eq. (180). 

The ground state is characterized by two filled "Fermi 
seas", one for the quasimomcnta i>i, and another for 
the spin quasimomenta ■ Similarly to the Lieb-Liniger 
model, excitations can be constructed by creating holes 
in these distributions. At zero magnetic field, the edge 
of support for the spectral function at \k\ < kp corre- 
sponds to a spinon excitation, where a hole is created in 
the spinon Fermi sea while a holon is created at the Fermi 
surface, in complete accordance with the field-theoretical 
description of Sec. II.F. Equations (180) and (182)-(183) 
contain the full information about the excitation spec- 
trum of both holons and spinous. In addition, finite- 
size corrections to their energies can be analyzed simi- 
lar to Sec. III.B. When combined with the extension of 
Sec. II. G, they lead to explicit predictions (Essler, 2010) 
for the phase shifts in terms of microscopic parameters, 
and the obtained results coincide numerically with the 
phenomenological predictions (127) and (128). This co- 
incidence provides a nontrivial non-perturbative check 
of the RG arguments of Sec. II.F, justifying the effec- 
tive Hamiltonians of impurities with fractional quantum 
numbers. 

Finally, let us mention that recent experimental 
progress with alkaline earth ultracold atoms (DeSalvo 
et al., 2010; de Escobar et al., 2009; Fukuhara et al, 
2007; Kraft et al, 2009; Stellmer et al, 2009; Taie et al, 
2010; Takasu et al, 2003) which naturally posses a higher 
symmetry of interactions (Cazalilla et al, 2009; Gorshkov 
et al, 2010) calls for extensions of the present approach 
to SU [N) invariant Sutherland-type models (Sutherland, 
1968). 



D. Lattice models: XXZ, spinless fermions and ID 
Hubbard model 

As has been discussed in Sec. II. D, the presence of 
a lattice leads to much wider possibilities for threshold 
behaviors. In this section, we will review some recent re- 
sults obtained by combining field-theoretical approaches 
with the exact solutions of XXZ (Cheianov and Pustilnik, 
2008; Karimi and Affleck, 2011; Pereira et al, 2008), spin- 
less fermion (Pereira et al, 2009) and ID Hubbard mod- 
els (Essler, 2010). The main modification for a generic 



non-integrable model at arbitrary filling is that, strictly 
speaking, the edges of support disappear due to the pres- 
ence of the lattice, as has been discussed in Sec. II. D. 
For intcgrablc systems, this might not necessarily lead 
to a smearing of the singularities. Nevertheless, to avoid 
this possible complication, we will discuss here models at 
commensurate fillings, such as half-filling. 

The XXZ model is given by Eq. (90), and its basic 
properties were discussed in Sec. II. D. Using a Jordan- 
Wigner transformation it maps onto a Hamiltonian of 
spinless fermions with nearest-neighbor interactions, see 
Eq. (94). The structures of the exact solutions of both 
models are the same and, e.g., S^^{q,uj) of the XXZ 
model coincides with the DSF of the fermionic model, 
while S~^~{q,uj) and the spectral function differ due to 
the Jordan- Wigner string. For concrctencss, here we 
will focus on the XXZ model and refer the reader to 
Ref. (Pereira et al, 2009), where spinless fermions have 
been discussed in detail. 

Similarly to the Lieb-Liniger model, the XXZ wave- 
function is written as a combination of plane waves (Ko- 
repin et al, 1993; Orbach, 1958). It is convenient to char- 
acterize them in terms of rapidities A which are related 
to the bare two-particle phase shift 6{X = Ai — A2) and 



the bare momentum f)o(A) as ^ = iln 



sinh(2ir^-t-A) 
sinh(2i77— A) 



PO = 



i In 



cosh(A— 177) 
cosh(A-|-iT/) 



, where r] conveniently parameterizes the 

interaction via A = — cos 2r]. The solutions of the Bethe 
equations in terms of the rapidities can be imaginary, 
which generally leads to a number of complications, such 
as the existence of bound states discussed in Sec. II. D. 
Inside the gapless regime, however, the ground state is 
constructed similar to the Lieb-Liniger model out of real 
solutions which occupy a "Fermi sea" (—A, A). Spin wave- 
like excitations also can be constructed by creating holes 
and adding particles with real rapidities on the top of 
the filled "Fermi sea" . The density of ground state roots 
p{X) and the shift function F^vlX) satisfy the equations 



1 



27r.y_A 



\i,.,^,)p{^)drJ=^^^. (184) 



277 dv 

6{u - A) 
27r ' 



(185) 



where K{v,ii) = d6{v — ii)/di', and the normalization 
condition for N spin-down particles on a lattice of size M 
reads J^^p{X)dX = N/M. Away from half-filling, the A 
following from this equation is finite and the leading non- 
linearity of the spin wave spectrum is quadratic; the re- 
sults of Sec. II. D. 3 are applicable. The peculiarity of the 
half-filled case (M = 27V), expected from the considera- 
tions of Sec. II. D, manifests itself in the exact solution as 
A — >• 00. In this case, all integral equations can be solved 
analytically by Fourier transformation, which leads, e.g., 
to an analytical expressions for the Luttinger parame- 
ter K = (2 — 2arccos A/tt)"^ and the edge of support 
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ULiq) = vsm{q) = ^arlcofl sin(g). As expected, the 
leading nonlinearity of spectrum is cubic. 

The central objects which determine the exponents of 
the response functions arc the phase shifts S±{k), which 
similarly to Eq. (176) are given by 6±{k) = 27ri^(±A|A). 
However, one needs to take the limit A — )• oo, and there is 
an ambiguity in the way this limit should be approached. 
This has lead to conflicting predictions by Pereira et al. 
(2008) and Cheianov and Pustilnik (2008). The ambi- 
guity was resolved by Imambekov and Glazman (2009a) 
in favor of the former, based on a comparison with the 
universal results of Sec. II. B and the SU{2) symmetry 
arguments of Sec. II. D for A = 1. The resulting phase 
shifts are momentum-independent, and they are given 
by 6-/{2'k) = -5+/{2-k) = 1/{2VK) - y/K/2. The 
exponents of S^'^{q,uj) and S'"' {q,'-^) can now be ex- 
plicitly evaluated using Eqs. (87) and (112), and are 
given by (Karimi and Affleck, 2011; Pereira et al, 2009) 
/ij = 1 — K, fi~ = 2 — — K. These exponents 
are momentum-independent and interpolate between the 
results of the XY and the XXX models of Sec. II. D. 

The results for the phase shifts and exponents away 
from half-filling can be obtained by numerically solving 
Eq. (185). The behavior of the response function near 
the energy of the bound state can be also analyzed: such 
a bound state merges with the spinon excitation at fi- 
nite momentum, and the field-theoretical description of 
the singularity changes at this point (Karimi and Affleck, 
2011; Pereira et al, 2009). 

Finally, let us also briefly comment on the application 
to the ID Hubbard model (Essler et al., 2005; Licb and 
Wu, 1968), which describes spinful fermions on a lat- 
tice. The Lieb-Wu system of equations which determines 
the energies and momenta of the eigenstates is similar 
to Eqs. (182)-(183), but with sint-j substituting Vi in- 
side the phase shifts. The finite-size corrections to the 
energies of holon and spinon excitations have been calcu- 
lated (Essler, 2010), and were used in conjunction with 
the field-theoretical results of Sec. II. F to obtain predic- 
tions for some of the threshold exponents. For certain 
values of the parameters they coincide numerically with 
the results of Carmelo et al. (2008, 2004, 2006) obtained 
using completely different methods. We note, however, 
that the detailed analysis of the ID Hubbard model is 
rather complicated due to the rich kinematics, and some 
of the first steps in this direction have been made re- 
cently (Pereira et al, 2012). 



IV. KINETICS OF AND TRANSPORT IN A NONLINEAR 
LUTTINGER LIQUID 

In this Section, we will review some elementary pro- 
cesses of relaxation, as well as kinetic and transport phe- 
nomena emerging in a nonlinear LL. In statistical me- 
chanics, one assumes that a generic macroscopic system. 



even if isolated from the rest of the world, will even- 
tually reach a local thermal equilibrium. The density 
matrix of a finite-size part of such a system will reach 
the Gibbs distribution as long as that part comprises 
many particles. The parameters of the equilibrium dis- 
tribution are fixed by additive conserved quantities (par- 
ticle number, energy, momentum). Normally, we expect 
the approach to equilibrium to be controlled by a spec- 
trum of relaxation rates found from an appropriate ki- 
netic equation (Huang, 1987). However, there are promi- 
nent counterexamples to that common wisdom. The ap- 
proach to thermal equilibrium of a system of interacting 
particles, in any dimension, may be hindered by disorder, 
resulting in a "many-body localization" . This possibility 
was raised by Anderson (1958), analyzed in the contexts 
of disordered solid-state conductors (Basko et al., 2006) 
and atomic cold gases (Aleiner et al, 2010) recently, and 
currently receives a considerable attention, see Pal and 
Huse (2010) and references therein. Closer to the sub- 
ject of this review, the abundance of integrals of motion 
in a disorder-free system is also deemed to prevent equi- 
libration (Polkovnikov et al, 2011). Such a possibility is 
foreseen in quantum integrable ID systems (Sutherland, 
2004). Quite remarkably, an experimental investigation 
of the latter subtle roadblock to equilibration became 
possible in the context of cold gases (Kinoshita et al, 
2006). A restricted phase space for scattering events 
suppresses relaxation processes even in a generic (non- 
integrable) ID system. Recently, some peculiar features 
of the electron equilibration were found in experiments 
with quantum wires formed within a GaAs heterostruc- 
ture (Barak et al, 2010) and carbon nanotubes (Chen 
et al, 2009). 

Related but not identical to the eqiiilibration problem 
is the question about singularities in the dependence of 
the response functions on momentum and frequency. As 
we saw in Sec. III.B, integrability allows the functions 
A{k,e) and S{q,u!) to be singular within the spectral 
continuum of excitations, in addition to the "manda- 
tory" non-analytical behavior at the thresholds of the 
continuum. If integrability is violated, the singularities 
within the continuum vanish even at zero temperature. 
At small e or oj, the singularities are smeared but not 
washed out completely, being rather replaced by some 
finite-width peaks. Like in the theory of Fermi liquids, 
these widths may be associated with the inverse lifetimes 
of the quasiparticle states which approximately diagonal- 
ize the many-body Hamiltonian of the nonlinear LL. Peak 
broadening of A{k, e) may be measured, in principle, in a 
tunneling experiment. We will identify some important 
elementary relaxation processes specific for various ID 
systems in Sec. IV. A. 

Equilibration processes and the dynamic density re- 
sponses of the liquid determine some of its transport 
properties. The most studied of those are the linear con- 
ductivity and the linear conductance of electron liquids 
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subject to an external electric field. In a ID system, the 
relation between the conductivity and the conductance is 
not trivial. The conductivity is well-defined in a contact- 
less setting, for a homogeneous liquid filling the entire ID 
space. Contrary to that, the conductance is determined 
as the current flowing through a system attached to leads 
biased with some small voltage. The conductance does 
depend on the properties of the leads. In fact, the linear 
LL theory predicts that the dc conductance is determined 
by the properties of the leads and is independent of the 
parameters of the LL (Maslov and Stone, 1995; Pono- 
marenko, 1995; Safi and Schulz, 1995). 

The conductivity a{q,Lo) of a homogeneous liquid is 
related by the Kubo formula to the current- current cor- 
relation function (Mahan, 1981). The real part cr'(g, co) = 
Iic(7(q,uj) of the conductivity can be expressed, with 
the help of the continuity relation and the fluctuation- 
dissipation theorem, in terms of the DSF. The expres- 
sion for the DSF obtained in the linear LL theory at any 
temperature results then in (j'(0,w) oc d{uj), commonly 
referred to as the Drude peak. For Galilean-invariant sys- 
tems it is not destroyed by spectrum curvature or finite 
temperatures, regardless of the interactions between par- 
ticles (Sirker et al, 2011). However, predicting its fate at 
finite temperatures and in the presence of a lattice is be- 
yond the linear LL description. The umklapp processes 
which are caused by the lattice and are formally irrele- 
vant at T = 0, may smear the (5- function singularity in 
a' at finite temperatures. We briefly review this question 
in Sec. IV.B. 

Equilibration processes do affect the conductance G of 
a ID electron liquid. These processes, absent in the lin- 
ear LL, make the conductance temperature-dependent. 
We review various elementary processes leading to equi- 
libration and their effect on the conductance and other 
transport characteristics in Sec. IV. C. 

Concluding the introduction to this section, we wish 
to emphasize that all of the questions raised here are 
beyond of the realm of the linear LL theory. The latter 
is trivially intcgrablc and easily mapped onto a system 
of free bosons or free fermions, so one does not expect to 
flnd any finite relaxation. 



A. Relaxation processes of excitations in a nonlinear 
Luttinger liquid 

Like in higher dimensions, it is instructive to start the 

consideration of relaxation processes in ID by discussing 
the case of almost-free spinless fermions. At zero inter- 
action, the single-fermion excitations are the true eigen- 
statcs with no degeneracies in the single-particle sector, 
and the ground state is not degenerate or pathological 
(unlike in the case of free bosons). This is helpful in 
building a theory of relaxation processes using perturba- 
tion theory in the interaction strength. The main part 



of Sec. IV. A. 1 is devoted to the identification and evalu- 
ation of the elementary relaxation processes for spinless 
fermions. We will see that the lack of particle- hole sym- 
metry leads to drastically different relaxation rates for 
particles and holes at low temperatures. We will also in- 
vestigate the peculiarities of the energy and particle num- 
ber transfer between the left- and right- moving species. 

Similar to the relaxation rates in higher dimensions, 
the perturbatively evaluated relaxation rate in ID van- 
ishes when the particle's excess energy tends to zero. Due 
to phase space constraints, the rate is proportional to a 
higher power of the excess energy, see e.g., Eq. (188). 
This should help in building a full analogue of the 
Fermi liquid and the kinetic theories of the quasiparticles 
emerging in the universal description of the nonlinear LL, 
sec Sec. II. B. Such a program for spinless fermions has 
not been performed yet. 

A step in that direction for an actually more com- 
plicated case of spin- 1/2 fermions is described in 
Sec. IV. A. 2. The complication arises from the spin de- 
generacy of the free-fermion single-particle states. A 
harbinger of the difficulties is already seen within the per- 
turbation theory: the scattering cross-section evaluated 
in the basis of free fermions is divergent at low energies, 
leading to a relatively slow dependence of the relaxation 
rate on the particle's energy, see Eq. (191) in Sec. IV. A. 1. 
We will see in Sec. IV. A. 2 that upon proper removal of 
the degeneracy and introduction of spinous and holons, 
the decay of the latter branch is efficiently suppressed. 

Methods built in Sec. IV. A. 1 to consider the relaxation 
of fermions help in the investigation of the relaxation in 
a ID Bose liquid. We move to ID bosons in Sec. IV. A. 3. 
For a weakly interacting gas, the relaxation of particle- 
like excitations can be understood with the help of per- 
turbation theory, with some improvements required for 
taking care of the strong modification of the low-energy 
excitation spectrum. The relaxation of the other impor- 
tant branch of excitations - dark solitons - turns out to 
be similar to the relax;ation of holes in a Fermi gas near 
the bottom of the band, but requires a non-perturbative 
treatment. 



1. Weakly interacting fermions 

In order to identify processes important in relaxation, 
we first turn to the case of spinless weakly interact- 
ing fermions. The curvature of the dispersion relation 
^(k) introduces particle- hole asymmetry into the prob- 
lem. For relaxation processes, the importance of particle- 
hole asymmetry is already seen within perturbation the- 
ory. Indeed, it follows from Eq. (2) that the hole veloc- 
ity is smaller than the velocity of low-energy excitations, 
i.e., particle- hole pairs near the Fermi points. There- 
fore, according to the Cherenkov radiation criterion, a 
hole introduced into the system cannot emit these exci- 
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FIG. 14 (Color online) (a) Relaxation of a high-energy par- 
ticle due to three-particle scattering (b) Relaxation process 
for high-energy holes at nonzero temperatures. Filled states 
are depicted in blue (darker), empty states in white (lighter) 
color. 

tations and consequently cannot relax at zero tempera- 
ture. On the other hand, a particle moves faster than 
the low-energy excitations. The emission of particle-hole 
pairs by a moving particle is therefore allowed by en- 
ergy and momentum conservation laws. The emission of 
a single particle-hole pair is identical to a two-particle 
collision. In this case, energy and momentum conserva- 
tion can only be satisfied if the two incoming particles 
with momenta ki and k2 either keep their initial mo- 
menta, {ki, ^2) — ^ {ki, fc2), or switch their momenta with 
the other particle, (^1,^2) — > (fc2,fci). Neither of these 
options can cause relaxation. 

Scattering processes that result in a redistribution of 
momenta and thus potentially lead to a finite relaxation 
rate must involve at least three particles. In such three- 
body collisions three particles with momenta k = kp +P, 
kR = kp +pr, and fc^ = —kp +Pl in the initial state \i) 
end up in a final state |/) with different momenta fc' = 
kp +p' , k'^ — kp +p'ji, and k'^ = —kp +p'i^, see Fig. 14a. 
For a generic interaction the transition \i) — >■ |/) has a 
nonvanishing momentum-dependent amplitude A. 

In order to evaluate the T = relaxation rate of an 
extra right-moving particle with momentum k = kp -\- p 
(0 < p <C kp) due to three-body collisions, we note that 
the single-particle states Pr,Pl in the initial state of the 
transition \i) are below the Fermi level, while all three 
single-particle states in the final state |/) are above it. 
Applying now Fermi's golden rule, we find 

—-ex/ dp'dp'j^dpL dp'LdpR\A\'' (186) 

5[{p + pr+pl)- {p' + p'r + p'l)] 
X <5| [£_{kp +p)+£_{kp +pr) + £.{-kp+pL)] 

~ mr + P') + i{kp + p'r) + i{~kp +p'l)]], 

where A is the three-body collision amplitude introduced 
above, and the ^-functions express the energy and mo- 
mentum conservation. 

In writing Eq. (186) we took into account that for 
p = k — kp <^ kp the conservation laws cannot be sat- 



isfied unless the collision involves both right- and left- 
moving particles.^ Further analysis shows that the con- 
servation laws allow a small (< p^/m) energy transfer 
to the left-movers. Such a solution can be found by 
iterations. To zero order in p^ — p'l, the momentum 
conservation gives p — p' = p'r — Pr- The energy re- 
leased in the collision of two right-moving particles then 
is ({kp +p)+ ({kp + pr) - ({kp + p') - ({kp + p'j^) < 
p^jm. This energy is transferred to the left-movers, 
({—kp +p^) — ({—kp -\- fci) < p^/m, which corresponds 
to the momentum transfer pl—p'l S P^I {"mvp) ^ P- Ac- 
cordingly, energy and momentum conservation restrict 
the range of the momenta contributing to the integral in 
Eq. (186) to 

p',p'r,\pr\<P, PL,\p'L\<pVimvp). (187) 

The (5-functions in Eq. (186) remove the integrations over 
p'j^ and p'j^ . The remaining phase space constraints yield 
integration domains ~ p^ /{mvp) for and ~ p for p' 
and \pr\^ see Eq. (187). These three factors (one oc p^ 
and two cx p) yield the estimate (Khodas et ai, 2007b) 

-^^\Amk)/vp]\ (188) 

For a weak generic interaction, the nonvanishing three- 
particle collision amplitude A appears already in the sec- 
ond order in the interaction strength. In the case of a 
long-range potential, allowing one to neglect terms pro- 
portional to V{2kp), the relaxation rate is (Khodas et ai, 
2007b) 

where C = i^Tr/{5 ■ 2^) « 0.06 and v is the density 
of states. For a potential falling off faster than 1/a;^ 
in real space, Vq — Vk-kp oc {k — kp)^, which leads to 
1/Tp{k) (X {k — kp)^. Perturbation theory in the irrele- 
vant interactions Eqs. (52) and (53) to the universal-limit 
Hamiltonian hints that such behavior persists beyond the 
perturbation theory. 

In the special case of an integrable model, one may ex- 
pect A to be identically zero (Sutherland, 2004). Within 
the lowest-order perturbation theory, it was checked 
by Lunde et al. (2007) for the Cheon-Shigehara model 
(Cheon and Shigehara, 1998, 1999) and by Khodas et al. 
(2007b) for the Calogero-Sutherland model that indeed 
1/Tp{k) = for these models. 

A vanishing relaxation rate would entail the presence 
of a power-law singularity in the spectral function A{k, e) 



^ The relaxation rates presented below in Eqs. (188)-(191) assume 
\£,{kp + p)\ <^ €p. In addition, we set a constraint T <C eF for 
Eqs. (190)-(193). 
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at the energy spectrum of a particle excitation. That sin- 
gularity lies within the spectral continuum. Apart from 
integrable models, however, 1/Tp{k) 7^ at finite k and 
therefore the particle peak in A{k, e) is broadened within 
the energy range defined by 1/Tp{k). We notice here, 
that in a generic case this range scales to zero faster than 
[k — kpY with k — > fci?, cf. Eq. (188), while the deviations 
from the linear spectrum occur on the scale {k — kF)^/'m. 
This justifies the consideration of power-law singularities 
in the limit k ^ kp. 

A finite temperature trivially broadens the singulari- 
ties in A{k, e) even within the linear LL description (Gi- 
amarchi, 2004). It would also broaden the singularities 
in S{q, uj) even in the absence of relaxation mechanisms 
(with a possible exception of the finite-temperature be- 
havior of S{q, cj) at (7 — > which we briefly review in 
Sec. IV. B). This broadening comes from the smearing 
of the edge of the Fermi distribution. Relaxation would 
manifest itself in the time evolution of the distribution 
function of excitations (thermalization) and in a num- 
ber of transport phenomena. Here we concentrate just 
on the elementary processes of relaxation of particles (p) 
and holes (h). 

Turning to the case of small finite temperatures, we no- 
tice that the above consideration of the zero-temperature 
particle relaxation rate remains valid as long as the par- 
ticle energy ^(fc) s/cpT, where = kp/{2m) is 
the Fermi energy. At smaller £,{k), the phase space of 
left-moving excitations participating in the collision (p^, 
Ip^I) is not controlled any more by the small transferred 
momentum of Eq. (187), but rather by thermal smearing 
~ T/vp of the momentum distribution function. As a 
result, the factor oc {k — kp)"^ coming from integration 
over PL is replaced by a factor oc mT, yielding 

— ^ (X \A\^mTUk)/vFf , T « C(fc) « 

Tp[K, 1 ) 

instead of Eq. (188). 

The finite-temperature effect is more dramatic for 
holes, since it makes their relaxation possible in the first 
place, see Fig. 14b. Due to thermal smearing, a coun- 
terpropagating particle can give up an energy of order of 
T. Thus, a hole can relax its energy with a character- 
istic energy loss of Ae ^ epT/\^{k)\. It means that an 
energetic hole "floats" towards the Fermi level in many 
steps small compared to |^(fc)|, as long as the hole energy 
remains large compared to \J epT . Under this condition, 
application of Fermi's golden rule yields the rate 

a l^lWe^^^, |C(fc)|»v^, (190) 

for a single step of the relaxation process; this rate defines 
the lifetime of a state with given energy ^(fc). 

Within the perturbative treatment, the evaluation of 
the relaxation rates was generalized to the case of spin- 



1/2 fcrmions by Karzig et al. (2010). Targeting the ex- 
periment by Barak et al. (2010), the evaluation was per- 
formed for electrons in a quantum wire of a small width 
a ^ i/kp interacting via a Coulomb potential which is 
screened by a gate at some large distance compared to a 
and l/kp. The relaxation rates were found to be 

(191) 

for particle excitations. Here, k is the dielectric constant 
of the host material, and A(^) = In |l/2A;Fa| In |^/4eir|. 
At finite temperatures and excitation energies |^(A;)| <^ 
\/epT, particles and holes relax with the same rate, 

\m\ « (192) 

where the numerical constant is ci = 41n2 — 1. 

Note that the rate Eq. (192) with |^(fc)| ~ T also 
determines the thermalization time of the system to- 
wards a boosted equilibrium distribution function [cf. 
Eq. (227)]. For the relaxation of spinless fermions such 
thermalization process was addressed rigorously by solv- 
ing a linearized quantum Boltzmann equation exactly in 
Ref. (Micklitz and Levchenko, 2011). Interestingly, the 
equilibrium is not characterized by an cqiiipartitioning 
of the injected excitation energy. Most of the injected 
energy rather stays within the, say, right moving branch 
because the energy transfer between right and left movers 
is suppressed as can be seen from Eq. (187) and Fig. 14: 
if a right moving excitation loses an energy ^ |^(fc)|, only 
a fraction ^ \S^{k)/eF\ of that energy is transferred to the 
left moving branch. 

The perturbative treatment of scattering of a spin-1/2 
electron off an electron in the Fermi sea requires that the 
incoming electron has energy ^(fc) ^ mvpV{q — >■ 0)/h. 
This is the applicability condition for the Born approxi- 
mation. One may view this condition as the one allowing 
the electron to preserve its integrity without separating 
into spin and charge modes in the collision process. In- 
deed, in the weak-coupling limit, the difference between 
holon and spinon velocities is Vc — Vg ~ V{q — )■ 0)/h, 
so one may recast the condition for the applicability 
of the perturbation theory as Vc — Vs ^ Aw, where 
Av = ^{k)/{mvp) is the difference of the velocities of 
the colliding particles; in other words, holon and spinon 
have no time to separate in the course of the electron 
collision (Karzig et al., 2010). We note that the pertur- 
bative result for particles at the boundary of its appli- 
cability, ^(fc) ^ mvpV{q 0)/h, matches the estimate 
of the holon relaxation rate evaluated in the limit of low 
energies, see Sec. IV.A.2 (Schmidt et al., 2010b). 

The asymmetry in the relaxation rates of particles and 
holes is a direct consequence of the nonlinearity of the 
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FIG. 15 (Color online) (a) Schematic representation of the ex- 
perimental setup used by Barak et al. (2010). A bias applied 
to lead 1 injects the current Ji through the left tunnel junc- 
tion into the grounded wire; the current I2 through the right 
tunnel junction is collected in lead 2. Depending on the bias 
polarity, particles or holes are injected. The injection occurs 
in a window of momenta (marked by shaded regions in panels 
(b) and (c)) around a value k controlled by the magnetic field 
B. (b) A hole (0 < fc < kp) injected from lead 1 cannot relax, 
and will be collected in lead 2. (c) The relaxation of a parti- 
cle (fc > kp) injected from lead 1 results in the formation of 
additional particle-hole pairs. Since only the particles are ex- 
tracted into lead 2, the collected current exceeds the injected 
current; the difference, drawn from the ground corresponds 
to the hole current sinking into the ground, see the dashed 
line in panel (a). 



excitation spectrum. It naturally explains the results of 
the experiment (Barak et al., 2010) in which electrons 
were injected in and extracted from a quantum wire. In 
the experiment, two tunnel junctions designed to have 
a momentum-dependent tunneling rate, were attached 
to a grounded quantum wire, see Fig. 15a. Because of 
the device constraints, it was possible to inject particles 
or holes within some band of momenta, with the cen- 
ter of the band controlled by a magnetic field applied 
perpendicular to the wires comprising the device. When 
holes were injected through the left junction, the cur- 
rent collected by the right junction was equal to the in- 
jected current (blue and red dots follow the same curve in 
the left portion of Fig. 16). That is naturally explained 
by the absence of hole relaxation: a single hole injected 
through the left junction is extracted with the right one, 
see Fig. 15b. However, once the junctions (and the ap- 
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FIG. 16 (Color online) Injected current I\ (dark dots) and 
collected current I2 (light dots) as a function of the magnetic 
field, which controls the momentum of injected carriers. The 
two currents coincide in the case of hole injection. If particles 
are injected, on the other hand, the collected current exceeds 
the injected one due to relaxation. The experimental setup 
and physical explanation are shown in Fig. 15. Adapted from 
Barak et al. (2010). 



plied injection bias) are tuned to inject and collect parti- 
cles, the collector current exceeds the injected one. This 
striking behavior can be explained by the relaxation of 
an injected particle, which creates a number of particle- 
hole pairs. Particles of these pairs are "scooped" by the 
collector, while holes are allowed to sink into the ground, 
see Fig. 15c. A simple set of rate equations explained 
quantitatively the observations (Barak et al., 2010). 

The relaxation processes considered above involve only 
low-energy excitations and do not change the numbers 
N]^ and iV^ of left- and right-moving particles. Changing 
those numbers bears consequences for the conductance, 
the thermopower, and the thermal conductance (Karzig 
et al, 2010; Levchenko et al, 2010, 2011a,b; Matveev 
et al, 2010; Micklitz et al, 2010; Rech and Matveev, 
2008; Rech et al, 2009). At low temperatures, the re- 
laxation of the difference Nb. — involves states close 
to the bottom of the band, see Fig. 17 (Lunde et al, 
2007; Matveev and Andreev, 2011). We define the corre- 
sponding relaxation time rjv by relation d{Nfi—NL)/dt = 
— {Nr — Nl)/tn, assuming that the temperature is the 
same for the left- and right-movers, while their chemi- 
cal potentials are slightly different. Because a "deep" 
hole is involved in the relaxation, the rate is exponen- 
tially small at low temperature, l/r^v c>c exp(— cf/T"). 
The pre-exponential factor scales as a power of temper- 
ature, with an exponent depending on the type of inter- 
action potential and the presence of spin degeneracy. If 
one assumes a smooth (in real space) potential and sets 
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FIG. 17 (Color online) A small-momentum relaxation process 
leading to a change in the numbers of left- and right-movers 



= ^0(1 ~ Q^/Qo) small q, while Vj,>fep — 0, then 
(Lunde et at, 2007; Mickhtz et al, 2010)" 




The temperature dependence of the pre-exponential 
factor comes from the phase space constraints on the 
scattering event (yielding a factor cx T"^), and from the 
partial cancelation of the direct and exchange contribu- 
tions to the scattering amplitude, similar to the one oc- 
curring in Eq. (189), which provides an additional fac- 
tor oc T^. Note that the latter factor is not present 
in higher-order terms with respect to the inter-particle 
interaction potential. Therefore, in the generic case 
1/tn oc exp[£th(0)/T], where £th(0) < is the en- 
ergy of a hole at the bottom of the band, renormalized 
by interactions. 

We should emphasize that the above estimates of ta? 
refer to an "elementary act" of changing Nh—Nl. In that 
act, a hole in the fermion distribution near the bottom 
of the band changes the direction of its motion. The 
exponential factor in 1 /tat comes from the probability for 
the existence of such a hole, and the prefactor comes from 
the inverse lifetime l/rdh oc of the existing deep hole, 
see Eq. (190). The characteristic variation of the hole 
momentum in the scattering event depicted in Fig. 17 is 
Ap ^ T/vf, while its energy variation ~ T"^ /ep is small 
compared to the characteristic change of energy (^ T) 
in each of the involved particle-hole pairs near the Fermi 
levels. This is why the hole dynamics may be viewed as 
diffusion in momentum space with the diffusion constant 

B ^ iAp)yTd^ « (194) 

and is described by a Fokker-Planck equation (Cas- 
tro Neto and Fisher, 1996). The proportionality coeffi- 
cient missing in Eq. (194) was found for the case of weak 
and strong interactions in a system of spinless fermions, 
respectively, by Micklitz et al. (2010) and Matveev et al. 
(2010). For arbitrary interaction strength, the results are 
discussed in Sec. IV. C. 



To conclude the discussion of relaxation processes 
within the perturbative treatment of interactions, we 
mention here a peculiarity of the scattering processes 
for spinless fermions on a lattice (Pereira et al., 2009). 
The free-particle spectrum ^{k) oc — cosfc of a tight- 
binding model allows for two particles with momenta k 
and fc3 = TT — fc to scatter into two other states, fci and 
^2 = TT — /ci. If the chemical potential is shifted from 
the middle of the band, this process apparently yields a 
finite decay rate for a range of possible particle momenta 
k within first-order perturbation theory. However, it is 
not fully clear if the lowest-order perturbative treatment 
is applicable in that special case: for a given state fc, 
the mentioned states fci, fc2, fcs involved in the relaxation 
process are also involved in the formation of two-particle 
bound states (Pereira et al., 2009). 



2. Spinful fermions at arbitrary interaction strength: holon 
lifetimes 

In this section, we will go beyond the weakly interact- 
ing limit for spinful fermions, and consider the lifetimes 
of holons in the low-energy limit. 

As we already mentioned, the phase space argument 
applied to interacting spinless fermions in ID leads to 
the estimate of the lifetime in Eq. (188). This estimate 
is valid at small energies of excitations whose dispersion 
relation resides within the particle-hole continuum. A 
similar argument applied to a decay of a holon into two 
spinous in a ID spin-1/2 fermionic system would lead to 
a decay rate oc |fc — fc^^j, possibly contradicting the notion 
of a well-defined holon branch at small k — kp- 

A combination of the methods described in Sees. II. B 
and II. C should allow us to express the proportionality 
coefficient in Eq. (188) in terms of higher derivatives of 
the dispersion relation with respect to momentum and 
particle density; similarly, these methods should allow 
one to reliably evaluate the decay rate of a holon. Such a 
program was not implemented yet for spinless fermions, 
but an attempt was made to evaluate the broadening of 
the holon branch of excitations in a spin-1/2 ID fermion 
system. There is some disagreement in the conclusions 
of Schmidt et al. (2010b) and Pereira and Sela (2010). 
We find that the decay rate of a holon close to a Fermi 
point {+kF, for definiteness) scales to zero faster than 
\k — kp\^ , as we discuss next. 

In order to elucidate the possible decay processes for 
holons, it is convenient to start again from a descrip- 
tion in terms of refermionized quasiparticles. The band 
curvature of the physical fermions leads to interactions 
between the quasiparticles. Away from the Fermi points, 
it is advantageous to classify the interaction processes by 
their relevance in the RG sense and to consider all possi- 
ble interaction operators which are allowed by SU{2)- 
symmetry and Galilean invariance. Due to its built- 
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in S'?7(2)-symmctry, non-Abelian bosonization (Gogolin 
et al., 1998) is a convenient tool to achieve this. Ex- 
pressed using the left- and right-moving holon densities 

Jq(.t) and spinon densities Ja{x) {a = L, R), the Hamil- 
tonian of the linear LL reads Hq = He 



Hg , where 



He = 2ttVc j dx[J%{x) + Jl{x)], 
H,^^ I dx[jlix) + Jl{x)]. 



(195) 

The holon and spinon densities are related to the physical 
fermion operators by 

(7 

Jaix) = ^tpl^{x)Saa''ipaa'{x), (196) 

and Sarr' denotes the vector of spin matrices (half of 
the Pauli matrices for spin- 1/2). The operators Ja{x) 
are related to the physical charge density by Pc{x) = 
1\fK~c {Jl{x) + Jr{x)). This Hamiltonian emerges at the 
low-energy RG fixed point and is valid in the narrow- 
band limit. The leading correction for increased band- 
width is an interaction between left-moving and right- 
moving spin densities, (Gogolin et al., 1998) 

Hg = -2TTVsg / dx Jr{x) ■ Jl{x). 



(197) 

Note that when expressed in terms of the Abclian spinon 
fields 0s and 9^, the operator Hg generates the sine- 
Gordon term (Giamarchi, 2004). The band curvature 
of the physical fcrmions leads to interaction operators 
which are cubic in spin and charge densities, 

Hr, = ^Jdx [r?_(4 + J£) - v+{JpL + JIJr)] , 

k_{JrJI + JlJI) 



H. = 



47r2 



dx 



+ k+{JrJI + JlJ%) 



dx {Jl + Jr)Jr ■ Jl- 
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Note that these operators represent all cubic terms which 
are compatible with SU (2)-symmetry. In particular, this 
symmetry prohibits terms linear in the vector opera- 
tors J„(.t). Interaction operators containing quartic and 
higher-order terms in Ja{x) and Ja{x) do exist but their 
contribution is subleading for small bandwidths. 

The prefactors g, C, k± and r]± can be fixed phe- 
nomenologically by relating them to other observable 
quantities. The modification of the constants of the 
Hamiltonian Eq. (195) in response to a small density vari- 
ation yields the relations (Pereira and Sela, 2010; Pereira 



et al, 2006; Schmidt et al., 2010b) 



K_ -|- 



(199) 



The difference — k+ can be related to the mass m of the 
physical fermions by considering a charge current varia- 
tion of the Galilean-invariant system. One finds (Nayak 
et al, 2001; Pereira and Sela, 2010) 



(200) 



It is known that upon a bandwidth reduction g flows 
logarithmically to zero (Gogolin et al., 1998). Assuming 
the initial bandwidth to be of order kp, for a smaller 
bandwidth of order k — kp the effective coupling con- 
stant will flow to g{k) = 1/ ln[kp /{k — kp)]. As the 
chemical potential /i is proportional to fc^, the derivative 
dg/dfi (X —g^jkp, as noted by Schmidt et al. (2010b). 
The derivative dvs/dfx, on the other hand, remains finite 
for small bandwidths. Therefore, in leading logarithmic 
approximation, dg/djj, can be neglected and the coupling 
constants k± and ( can be related as 



-g{K- +K+). 



(201) 



Holons can relax via the creation of low-energy spinous. 
Let us investigate the decay of an initial state \i) = 
|fc)c|0)s which contains an additional holon with momen- 
tum above the Fermi edge and no spinon excitations. Re- 
laxation of the holon to a momentum k' < k can happen 
via the creation of two spinon density excitations with 
momenta pl < and pR > 0. This final state will be 
labeled |/) = \k')c\PL,PR)s- For momenta k close to the 
Fermi point, momentum and energy conversation for this 
process read 

k = k' +PR+ PL, Vck = Vck' + Vs{pr - Pl), 

and have nontrivial solutions (k ^ k') for Vc > Vs- 

The holon lifetimes associated with this decay channel 
can be calculated using Fermi's golden rule. Two com- 
binations of operators from the interaction terms (198) 
have a nonzero matrix element between the states \i) 
and I/). To first order in the interaction, {f\H(^\i) is 
the only such term. To second order, only {f\HgHK,\i) 
and {f\Hf^Hg\i) are nonzero. 

For the first-order matrix element 
finds (Schmidt et al., 2010b) 



{f\Hc\i), one 



2l' 



'k-k' 



PL -PR 



(202) 



The matrix elements {f\Hg\i) and vanish be- 

cause Hg and H^ do not couple spinous and holons. The 
remaining first-order matrix element {f\HK,\i) = be- 
cause it contains only terms of the form J^{x), which do 
not create spinous on opposite branches. 
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To the second order, cross-terms of the operators Hg 
and may couple the same initial and final states as 
above, yielding the amplitude 



k — k'—pL—PR 



PlPrI 



(203) 



Other second-order terms exist but they contain higher 
powers of PL and pn and are therefore subleading com- 
pared to Treg for holon momenta k near the Fermi points. 
According to Fermi's golden rule the rate is 



Tholo 



= 2TrJ2\Tc+T.g\^Sief-e.), 



(204) 



l/> 



where ej and e/ are the energies of th(; initial state |i) 
and the final state \f), respectivc^ly. The sum over all 
final states |/) translates to a summation over the mo- 
menta PL < 0, pn > and k' e [kp,k]. It can be 
seen from Eqs. (202) and (203) that each of the decay 
channels taken individually would lead to a decay rate 
1/Thoion oc {k — kp)^- However, Fermi's golden rule (204) 
contains the square of the sum of the probability ampli- 
tudes and T^g. The prefactors of both amplitudes are 
related according to Eq. (201) and one finds T^+T„g = 0. 
Therefore, the decay rate vanishes^ up to terms propor- 
tional to g^{k — kp)^, in the calculation of l/rhoion per- 
formed to the second order in g = 1/ lii[kp/{k — kp)]. 
Retaining in Eq. (199) the derivative dg/dfj. oc g^/cp ex- 
ceeds the accuracy of our calculation. It is not clear if 
the evaluation of l/rhoion to order g'^ would yield zero. 
Possibly, in that order the distinction between integrable 
and non-integrable systems emerges. 

In the limit of weak backscattering, V{2kp) <C V'(O) <C 
vp, the universal logarithmic dependence for g(k — 
kp) is reached only at very low energies, while its 
bare value g oc V{2kp)/vp is applicable as long as 
[V{2kp)/vp]ln[kp/{k - kp)] < 1. In that limit, which 
includes weak Coulomb repulsion, Eqs. (199) and (203) 
yield 



'Hiolon 



OC ep- 



Vp 



V{2kp) 



Vp 



k — kp 



ki 



(205) 



This estimate should be viewed as the result of perturba- 
tion theory in V{2kp) in the basis of well-defined holon 
and spinon modes with linear spectrum, which sets a 
limit on the holon momenta, k — kp < ml/(0) (we also 
used Vc — Vs ~ V'(O) in the derivation). Curiously, the lat- 
ter estimate for l/rhoion at the limit of its applicability, 
k — kp ^ mV{0), matches the estimate of the relaxation 
rate of a spinful fermion evaluated in the basis of free 
fermions perturbatively, see Eq. (191). 



3. Relaxation of excitations in a weakly interacting ID Bose gas 

Within the integrable Lieb-Liniger model, see 

Eq. (159), the excitations of a ID Bose gas do not relax. 
The DSF exhibits a power-law singularity at the Lieb-I 
mode and a power-law behavior converging to zero at the 
Lieb-II mode, see Sec. II. E. As discussed in Sec. III.B, 
in the limit of weak interactions the dispersion relation 
for the Lieb-I mode approaches the Bogoliubov quasi- 
particle spectrum, while the dispersion of the Lieb-II 
mode corresponds to the spectrum of "dark" solitons. 
A weak perturbation breaking the integrability leads to 
finite lifetimes of the excitations (Muryshev et at, 2002). 
At zero temperature, the finite relaxation rate of the Bo- 
goliubov quasiparticles smears the singularity in the re- 
sponse functions at the Lieb-I mode (Tan et ai, 2010). 
At T ^ 0, the relaxation rate of a dark soliton prepared 
in some high-energy state also becomes finite (Gangardt 
and Kamenev, 2010; Mazcts et al., 2008; Muryshev et ai, 
2002) . The relaxation rates of Bogoliubov quasiparticles 
and dark solitons in ID strongly depend on temperature. 
The theory of dark soliton relaxation was also extended 
to include the dissipative dynamics of the so-called de- 
pletons forming around an impurity imbedded in a ID 
Bose gas or a spin-fiipped particle in a spinor ID Bose 
gas (Gangardt and Kamenev, 2009; Schecter et ai, 2011). 

The leading corrections to the Lieb-Liniger model 
(159) have the form of a three-body interaction term. 



V = 



a 

9m 



dx : p^{x) : 



(206) 



These terms of the Hamiltonian H = ffLiLi + V can be 
derived explicitly by a projection onto the lowest sub- 
band of transverse quantization in a confining poten- 
tial with cylindrical symmetry. For a model in which 
the interaction in 3D is described by a pscudopotential 
VaDir) = 47r(a/m)(5(r), where a is the scattering length 
(Pitaevskii and Stringari, 2003), and with the amplitude 
of radial zero-point motion = (mur)^^^'^ ^ a one 
finds 4 (Mazets et ai, 2008; Muryshev et al., 2002; 01- 
shami, 1998; Tan et al, 2010) 

7 = 2mc/p = 2a/ [pal) , a = 181n(4/3)(a/ar)^ . (207) 

Here, p is the density of ID Bose gas. The limit of weak 
interaction means 7 <C 1. 

A finite three-particle scattering amplitude which leads 
to a damping of the Bogoliubov mode appears already in 
the first order in a ^ 1. The evaluation of the corre- 
sponding relaxation rate is especially simple for quasi- 
particles with energies ei{q) ^ 'ymp'^, so that Ei{q) i=s 



^ This conclusion of Schmidt et al. (2010b) differs from the one of 
Pereira and Sela (2010). 



* The value of a in (Mazcts et al, 2008) contains a spurious factor 
of 4 compared to Eq. (207). 
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I2m. In addition, wc assume £i{q) T. To the low- 
est order in a, the differential rate of inelastic scattering 
in which a quasiparticle with momentum q loses energy 
u is given by 

^,(w) = ^^ j dp5[uj-e,{q)+e,{q-p)]g{p,co) . (208) 

The Fourier transform g{p,co) = J dxdte'"^-'P^g{x,t) of 
the correlation function 



small energy transfer, < max{T. 7p^/to}. These pro- 
cesses are Bose-enhanced by the high occupation factors 
of the final states of the quasiparticles receiving the en- 
ergy uj. Finding the full dependence of l/r^ on temper- 
ature within this range would require a full knowledge 
of the correlation function Eq. (209), which is still not 
available. However, matching the results obtained at low 
and high values of T/Tg, one may see that l/r^ reaches 
a maximum (Tan et al, 2010) 



g{x,t) = {:p^{x,t) ::p'(0,0) :) 



(209) 



should be evaluated for the Lieb-Liniger model, 
Eq. (159). In terms of (Jq{oj), the inverse lifetime for 
a given momentum q is given by 



1 



dujaq{u)) . 



(210) 



The set of equations (208)-(210) allows one to evalu- 
ate the temperature dependence of l/r,. The charac- 
teristic temperature scale for the variation of 1 /r^ is the 
quasi-condensation temperature Tg = accessi- 
ble experimentally (Armijo et al., 2011). By the order-of- 
magnitudc, this is the temperature at which the chemical 
potential of the ID interacting bosons crosses zero. 

In the limits of low (T < T^) and high (T » T^) tem- 
peratures, one may use the proper asymptotes of the cor- 
relation function Eq. (209) to evaluate a-g(w) and l/r,. 
It turns out that in these limits the lifetime given by 
Eq. (210) is controlled by scattering processes with en- 
ergy transfer in a broad range, w < ei{q)/2 and is inde- 
pendent of q (Tan et al., 2010) 



3\/3 



m 



92{T) . 



(211) 



Here 52 = (: p^(0,0) is the two-particle correlation 

function (normalized by p^). For weak interactions 52 
decreases monotonically with T from (72 = 2 at T ^ T, 
to ff2 = 1 at T < (Kheruntsyan et al, 2003). The 
presence of 32 in Eq. (211) is due to the fact that the 
two particles receiving the energy u! in the three-particle 
collision must be at the same spot. 

A more detailed analysis actually indicates that 
Eq. (211) yields the dominant contribution to the relax- 
ation rate only outside the range of temperatures 



7' 



,3/8 



£1(9) 



1/4 



< ^ < ^-3/4 



1/2 



(212) 



That range is broad as long as the energy of incoming 
quasiparticle is not too high, ei{q)/Ts <C 8/7^/^, and 
includes some parts of temperature intervals where, re- 
spectively, T <SiTs, and T ^ Tg. Within the range given 
by Eq. (212), relaxation is dominated by processes with 



max 
Q 



m y 7~^/^£i(g) 



(213) 



at Tmax ~ 1.6Ts. Under the assumed condition on £i(g), 
the temperature dependence of l/r^ is not monotonic, 
the maximal rate (213) significantly exceeds the limits 
given in Eq. (211). The kinetic equation accounting for 
the small-energy transfers effective in the temperature 
interval (212) was considered by Mazets (2011). 

The relaxation of dark solitons in the interacting ID 
Bose gas is very similar to the relaxation of holes in the 
interacting Fermi gas which we considered in Sec. IV.A.l. 
Dark solitons correspond to the excitations over the 
ground state with the minimal energy £2 (g) at given mo- 
mentum q. The soliton velocity Vg {q) at any q is smaller 
than the soimd velocity in the Bose gas, v — (p/m)y^, 
see Eq. (171). Therefore, the relaxation of solitons is pos- 
sible only at finite temperatures in a Raman-like process 
in which two phonons arc created (each of the two partic- 
ipating phonons replaces a low-energy particle-hole pair 
in the case of the fermionic hole relaxation, see Figs. 14 
and 17.) 

The two-phonon processes lead to a typical momen- 
tum transfer Aq ~ T/v. The soliton's velocity is zero 
at q = irp, so the transferred energy in the elemen- 
tary act of relaxation of a dark soliton is on the order 
of T^/{\M*\v'^) ~ T'^/Tg and much smaller than T. Here 
M* = —Ap/v is the soliton's (negative) effective mass; 
following Gangardt and Kamenev (2010), we consider 
T <C Tg. The smallness of the energy allows one to use 
the Fokker-Planck equation (Landau and Lifshitz, 1980) 
to describe the time evolution of the momentum distri- 
bution function f(q) of a dark soliton. The problem is 
similar to the previously considered kinetics of a heavy 
particle in a linear LL (Castro Neto and Fisher, 1996) 
and to the diffusion of a deep hole in an interacting elec- 
tron gas, see e.g., Matveev and Andreev (2012b). Using 
the notations of the latter work, we write the Fokker- 
Planck equation in the form of a continuity condition for 
the distribution function in the momentum space, 



dtf 



-dqj , J 



B{q) 



1 dE 
T~dq 



f (214) 
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with J being the corresponding current 

The apphcabihty of Eq. (214) is confined to the vicinity 
of g = Trp in order to satisfy the requirement of the small- 
ncss of the energy transfer. Thus, one may use the ex- 
pansion £2((?) = £2(0) + (g — 7rp)^/(2M*) for the soUton's 
energy, and replace B{q) by a constant, B = B{q = wp). 
After that, the meaning of Eq. (214) becomes quite clear: 
it describes the motion of a "particle" (dark soliton) sub- 
ject to a Langevin random force (yielding the second term 
in the brackets) and a viscous force 



F = -XV, X = B/2T, v = d,e2{q). 



(215) 



The viscous force leads to a particle acceleration as M* < 
for the dark soliton. Once the soliton's velocity Vs{q) 
becomes of the order of v, Eq. (214) becomes invalid. 
However, if the initial velocity was on the order of the 
thermal one, Vs (q) a/|M*|T, then the Fokkcr-Planck 
equation describes the longest part of relaxation process, 
which takes a time t ^ \M*\/x- 

The viscosity coefficient was evaluated by Gangardt 
and Kamencv (2010), ^ 



X{T) 



1215 



T 



^ ,T«;^T,. (216) 



To obtain this result one needs to include not only lin- 
ear but also quadratic terms in dx^f and d^d in the in- 
teraction Hamiltonian of the quantum impurity with the 
LL describing the low-energy excitations. Such terms are 
necessary because they account for the three-particle col- 
lisions which are needed to capture two-phonon processes 
and the resulting soliton relaxation. At q = np the terms 
allowed by symmetries are 

H3 = Jdx [Veeid^ef + V^^id^^f] d\x)d{x). (217) 

An extension of the phcnomcnological approach of 
Sec. II. C leads to the couphng strengths 



V, 



1 / 1 d^e2 



96 



2 \m dq^ 



27r2 V 



(218) 
(219) 



The evaluation then proceeds by removing from the 
Hamiltonian the terms linear in d^'^ and d^O using the 
unitary transformation (71), and then treating the re- 
mainder within perturbation theory. As was shown ex- 
plicitly by Gangardt and Kamenev (2010), the coeffi- 
cients of the interaction Hamiltonian yield a vanishing 



^ The function B{q) at arbitrary interaction strength was recently 

evaluated (Matveev and Andreev, 2012a). 
® The numerical coefficient in Eq. (216) corrects an error in the 

original publication (Kamenev, 2011). 



relaxation rate for the integrable Lieb-Liniger model. In 
the weakly interacting regime, the lowest-order correc- 
tion appears in the order due to corrections to e2{q) 
coming from three-particle interactions (206). Later on, 
the outlined approach to the relaxation of dark solitons 
was generalized to the case of a depleton, the dressed im- 
purity state in a quantum liquid (Schecter et al., 2011). 



B. Conductivity and Drude weight for interacting particles 

So far in this section we concentrated on the relaxation 
of excitations. A related question which has attracted a 
lot of attention recently is the effect of such relaxation on 
transport, and specifically the relation between transport 
properties and integrability (Grossjohann and Brenig, 
2010; Heidrich-Meisner et al., 2003; Herbrych et al, 2011; 
Jung et al., 2006; Jung and Rosch, 2007; Karrasch et al., 
2011; Prosen, 2011; Rosch and Andrei, 2000; Sirker et al, 
2009, 2011; Znidaric, 2011; Wu and Berciu, 2011). 

The current response j{x,t) to a force f{x,t) — 
fq,ui exp{iojt—qx)e^^* applied to a linear hh is easily eval- 
uated by solving the corresponding classical equation for 
the density waves in the liquid. It yields for the Fourier 
components of the current 



i(9,w) = -i 



Kv /, 



TT u! — qv — iO 



(220) 



In a Galilean-invariant system, the factor Kv/tt equals 
p/m, where the average particle density p and the mass 
m are independent of interactions. The linear response 
of the current to a force field applied to a spatially- 
homogeneous system is characterized by the conduc- 
tivity, j{q,uj) = (T(q,uj)fq_^, which is a complex fmic- 
tion, a{q,Lj) = a'{q.uj) + ii7'\q,uj). Its real component, 
a'{q,oj), is the dissipative part of the response. The re- 
sponse of a Galilean-invariant system at (jr = is purely 
inertial and independent of interactions, since it is noth- 
ing but a center-of-mass motion caused by an applied 
force uniform in space. According to Eq. (220), the 
corresponding dissipative part of conductivity, cr'(w) = 
(T'{q,bj)\q=o, is 



a'{ixi) = 2-kD5{uj). 



(221) 



The magnitude of the "Drude peak" D in the conduc- 
tivity <j' {<jj) is given hy D = p/{2m). In the presence of 
a periodic lattice potential, a{q,co) is still well-defined 
for wave vectors much smaller than the size of the Bril- 
louin zone (Sirker et al, 2011). By continuity equa- 
tions, the conductivity is related to the susceptibility; 
the dissipative part of the latter is related to the DSF by 
the fluctuation-dissipation theorem (Doniach and Sond- 
heimer, 1998). Therefore, for small q 



Siq,uj) 



2q^ 



w(l - e-l^'^) 



(j'{q,Lo). 



(222) 
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As wc have seen in the previous parts of the review, the 
delta peak in S{q,oj) at finite wave vectors q becomes 
broadened in a nonlinear LL, and according to Eq. (222) 
leads to a finite width of the peak in the dissipative con- 
ductivity. At T = 0, the peak at lv = vq has a width 
which scales as a higher power of wave vector, oc q'^ in 
the absence of particle-hole symmetry or cx \q\'^ in the 
presence of the symmetry (Khodas et al., 2007a; Pereira 
et al., 2006). In either case, taking the limit q ^ one 
recovers the Drude peak in the conductivity in the limit 
0. 

At T 7^ 0, the universal nonlinear LL theory leading to 

Eq. (147) is insufficient for understanding the fate of the 
Drude peak. Indeed, the applicability of Eq. (147) in the 
limit g' ^ requires that T scales to zero not slower than 
q: one is not allowed to take the limit g — > at fixed T. 
However, in a Galilean invariant system, the existence of 
a Drude peak is protected even at finite temperatures, 
since the constant uniform external field causes the same 
center-of-mass motion irrespective of the temperature. 

The real part of the conductivity (at q = 0) can be 
written as (Sirker et al., 2011) 



a'{co) = 2'JtD6{uj) + f7reg(w). 



(223) 



Invoking the notion of a finite relaxation time r, one 
may expect cx (1/t)/[cj^ -|- (1/t)^]. In order to 

reproduce the correct zero-temperature limit, we have to 
set 1/t(T) = at T ~ 0. In a generic system, one expects 
1/t{T) finite at T 7^ which means zero Drude weight, 
D{T) = at any finite temperature. In the special case 
of integrable models, one may conjecture 1/t{T) = and 
consequently D{T) ^ even at finite temperature. 

One way of checking this conjecture relies on the rig- 
orous Mazur inequality (Mazur, 1969; Zotos et al., 1997) 



D > 



2LT 



E 



(224) 



Here L is the length of the system, I is the spatial in- 
tegral of the current density operator, and the operators 
Qfc form a set of commuting conserved quantities, or- 
thogonal to each other {{QiQk) oc 5ki). Moreover, if the 
set includes all conserved quantities Qf., then equality 
is reached in Eq. (224) (Suzuki, 1971). Finding at least 
one conserved quantity Qn with a nonzero overlap with 
I allows one to prove D{T) 7^ at finite temperature. 
This is the case, for example, for charge transport in the 
Hubbard model away from half-filling and spin transport 
in the S' = 1/2 XXZ model at finite magnetic field, where 
a proper local conserved quantity can be found. At zero 
field, all local conserved quantities for the XXZ model 
are even under the transformation 5| ^^j^ ^ '^f^ 
while I is odd, and the corresponding overlaps equal zero. 
It seems that a conserved quantity which has finite over- 
lap with X was recently constructed (Prosen, 2011), sup- 
porting the claim of finite D(T) for the XXZ model at 



zero field. Recent numerical results seem to support this 
claim (Karrasch et al., 2011). 

Alternatives to the investigation methods based on 
Mazur's inequality are reviewed in a recent excellent pa- 
per by Sirker et al. (2011). 



C. Conductance of interacting fermions in ID 

The conductance of ballistic quantum wires adiabati- 
cally connected to leads is quantized in units of e^/ {ttH) 
at low temperatures. This experimental observation (van 
Wees et al., 1988; Wharam et al., 1988) was first under- 
stood in terms of adiabatic transport of free fermions 
(Glazman et at, 1988). It was realized later that inter- 
actions, taken into account within the framework of the 
linear LL theory, do not affect the quantization of adia- 
batic transport (Maslov and Stone, 1995; Ponomarenko, 
1995; Safi and Schulz, 1995). Finite-temperature correc- 
tions to the quantized conductance, whether in the pic- 
ture of free or interacting fermions, are associated with 
the electron states near the bottom of the conduction 
band. In the case of free fermions the correction is easily 
evaluated and shows no dependence on the wire length L, 
but an activated temperature dependence 5G oc e~^^l'^ , 
see Eq. (226) below. Interactions do not alter the ac- 
tivated nature of the temperature dependence for rela- 
tively short wires (Lunde et al., 2007). In sufficiently 
long wires, however, equilibration facilitated by the scat- 
tering of holes near the bottom of the band, sec Fig. 17, 
ultimately leads to a much bigger correction, oc 
(Rech et al., 2009). We will mostly concentrate on the 
transport of spinless fermions below; one may think of 
fully spin-polarized electrons, the corresponding conduc- 
tance quantum is /{2-Kh). 

In the absence of interactions, the distribution func- 
tions of fermions in the wire, see Fig. 18, keep the mem- 
ory of the distribution in the lead they originated from, 



/W(fc) = 



-+ 



0{-k) 



exp(e^(fc)/T) + l exp(^«(fc)/T) + l' 



(225) 



Here ^^'^{k) = k"^ /{2m) — e^'^ are the energies of elec- 
trons coming from the left (L) or right (i?) leads, re- 
spectively. The difference between the chemical po- 
tentials is determined by the bias V applied to the 



wire, 



ep + eV . Evaluating the current / 



e / {dk/2TTh){k/m)f^°^{k) at smah bias {V -> 0), one eas- 
ily finds 



Go 



27rfte-«F/'^-|-l' 



(226) 



which is equal to the quantum e^/(27rfi), up to a cor- 
rection proportional to e~'^^l^ , which is exponentially 
small at low temperatures (e^ = = at F = 0). 

The presence of a current / 7^ means that there is 
some finite average velocity of electrons in the wire. An 
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while its momentum is 




FIG. 18 Schematic picture of the quantum wire of length L 

which is formed hy confining a 2D electron gas with gates 
(dark regions). Electrons in the left and right leads are 
described by Fermi distribution functions characterized by 
a temperature T and chemical potentials tp and ep, re- 
spectively. As particles propagate from the left lead to 
the right one, some get reflected due to the momentum- 
conserving electron-electron interaction, creating dNR/dt ^ 
0, see Eq. (233). Adapted from Micklitz et al. (2010). 



equilibrium distribution function in the rest frame, at 
given drift velocity, chemical potential, and temperature 
{u, T, and jj, respectively) would be 



m = 



exp[/c2/(2TO) - uk - il]/f + I 



(227) 



At r = T = 0, the distribution function (225) may be 
brought to the form of Eq. (227). Therefore, we do not 
expect equilibration (caused by electron-electron interac- 
tion) to bring any corrections to the ballistic conductance 
at T = 0. Besides, each of the two parts of Eq. (225) rep- 
resents a distribution describing equilibrium separately 
within the left- and right-movers with N]^ ^ N^. At 
T = 0, the relaxation rate I/tat = 0, and the distri- 
bution (225) is stable. However, at finite temperature 
the processes depicted in Fig. 17 cause a redistribution 
between N^. and N^. To the lowest order in I/tat, the 
correction 6G ^ L/{vpTjq) to the ballistic conductance 
(226) increases linearly with the wire length L (Lundc 
et al, 2007; Micklitz et al, 2010). This defines a new 
characteristic equilibration length, Zoq ~ vftn oz e'''^'^ . 
The resulting resistance reduces G by an L-indcpendent 
amount, SG ^ —{e'^ /h){T/ep)'^ which can be found es- 
sentially from particle number and energy conservation 
laws (Rech et al, 2009). The initial consideration of 
weak interactions (Rech et al, 2009) was generalized later 
to the case of strong (Matvccv et al, 2010) and arbi- 
trary (Matveev and Andreev, 2012b) interactions. We 
address next that latest development. 

The Hamiltonian (41) in a finite-size system written 
in terms of individual bosonic modes reads (Haldane, 
1981b) 



(228) 



(229) 



where vj ~ vp and wat = vp/K'^ as a consequence 
of Galilean invariance. An eigenstate of the system is 
described by boson occupation numbers, and the total 
numbers of the left- and right-movers with respect to 
the ground state, Np^^^ = {N ± J)/2. It is worth not- 
ing that an excitation of the boson modes does not con- 
tribute to the current: a creation of bosons corresponds 
to exciting particle-hole pairs within the branches of left- 
or right-movers. So, the electric current I is related 
by / = evj{J/L) to the difference Np — Np (Haldane, 
1981b). On the other hand, conservation of energy and 
momentum determine the form of the equilibrium dis- 
tribution, e-'-"-'^^^/'^ /Z with some parameters u and 
T (here Z is the proper partition function). Using here 
Eqs. (228) and (229), we find the average value of J to 
be irkpLu/vj. Using that together with the relation be- 
tween / and J, we see that u is nothing but the drift 
velocity: 



I 1 



e nkp 



(230) 



In equilibrium, the very same velocity "shifts" the boson 
distribution function: 



niq) = {h\hq) = \e^^\i\-^i)/T - l 



1 -1 



(231) 



That distribution does not carry any charge, but does 
create energy current, 



Je = {7r /'S){Tyv)u. 



(232) 



As we already discussed, at zero temperature the dis- 
tribution e~*^^~"^^/^/Z with a finite drift velocity can 
be viewed as two coimter-propagating fluxes of particles 
with different chemical potentials, resulting in the quan- 
tized conductance Go = e'^/{2Trh) (Maslov and Stone, 
1995; Pononiarcnko, 1995; Safi and Schulz, 1995). At 
finite temperatures, the accommodation of the distribu- 
tions in the leads to the drifting one in the wire may 
cause backscattering, see Fig. 18, which leads to relation 



GoV 



dNp 
' dt 



(233) 



The crucial observation is that dNp/dt is related to 
the energy redistribution between right- and left-movers 
(Matveev and Andreev, 2011; Rech et al, 2009). Indeed, 
backscattering of a right-mover corresponds to a change 
ANfi = — 1, which in the limit of w — )■ (linear response 
regime) does not affect the total energy, see Eq. (228). At 
the same time, by momentum conservation, bosons must 
acquire the additional momentum of 2kF, see Eq. (229). 
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This is only possible if momenta and energies given by 
^Pr,l = kp and AEr^l = ivkp are transferred to the 
left- and right-moving bosons. Therefore, 



dER 



-vkp 



dNR 



dt dt 
By energy conservation, one has 



JE = 



dE. 



dt 



(234) 



(235) 



Equations (235), (234), and (232) relate dN^/dt to the 
drift velocity u. Using then relations (230) and (233) , one 
finds the corrected conductance, (Matveev and Andreev, 
2011) 



G = Go 1 



3 v^kjr 



(236) 



A full equilibration of the bosons to the distribution 
(231) requires energy equilibration and the adjustment 
of the velocity to the correct drift value (230). As 
we discussed for the example of weak interactions, see 
Eqs. (188)-(192), the energy equilibration rate scales as 
some power of temperature, while the adjustment of the 
velocity requires a variation of Nr — Nl. The correspond- 
ing rate has an activated temperature dependence and 
happens on a much slower scale, see Eq. (193) and the 
discussion around it. Therefore, there is an exponentially 
wide interval of wire lengths L for which full equilibra- 
tion of the energy does occur, but u{L) does not reach 
the value (230). Considering u{L) as an adjustable pa- 
rameter replacing the distribution (231), and using the 
momentum conservation Eq. (229), one finds (Matveev 
and Andreev, 2011) 



1 dNu 
L dt 



u{L)-u 



3 v^kf 



(237) 



eq 



which generalizes Eq. (236) to finite values of L/l^q, 



G = Go 1 



TT^^ L_ 

T v'^k% L + k 



(238) 



The evaluation of the equilibration length Zeq = '^vtn 
is beyond the linear LL theory, and should account for 
the processes involving holes at the bottom of the band, 
see Fig. 17. Like in the case of dark solitons considered 
in Sec. IV. A. 3, this problem can be reduced to the one of 
the kinetics of a mobile impurity and entails a solution 
of the Fokker-Planck equation (214). It yields (Matveev 
and Andreev, 2012b) 



eq 



(3/2)fc|g 
n'^vV2TrM*T 



£th(0)/T 



(239) 



where the parameters eth(O) < and M* > are deter- 
mined by the energy spectrum of the impurity, £th(fc) ~ 



£th(0) + k'^ /(2M*), and depend on the equilibrium elec- 
tron density p in the wire. The extension of the phe- 
nomenological treatment of Sec. II. C aimed at including 
the Raman scattering processes, sec Eqs. (218), allows 
to express the coefficient B in terms of the functions 
deth{0)/dp and dv/dp: 



B = 



+ 



Aw^p'^T^ f d2£th(0) ^ 2 dv rfeth(O) 



(&th(0)/rfp)' 

M*u2 



dp^ 



V dp dp 



(240) 



The approach leading to Eq. (238) was also applied to 
spin-1/2 fermions (Matveev and Andreev, 2011). In the 
most interesting case of strong interactions {Vg <C Vc) the 
result is 



G 



7rh\ &vlklL + &^ 



eq 



but the corresponding equilibration length Zeq for spinous 
has not been evaluated yet. 



V. CONCLUSIONS 

The linear Luttinger liquid (LL) theory has been in 
use for decades by now, and provides an effective tool 
to describe the low-energy properties of ID quantum liq- 
uids in terms of quantized linear sound modes. Despite 
its spectacular successes, the linear LL theory has limita- 
tions constraining its applicability even in the low-energy 
physics of quantum ID systems. Replacing the generic 
spectrum of particles with a linear one does affect quali- 
tatively the momentum-resolved dynamic responses and 
introduces an artificial particle-hole symmetry. Further- 
more, being a free-field theory, the linear LL description 
is devoid of any intrinsic mechanisms of relaxation and 
equilibration. 

This review exhibits ways of studying ID quantum liq- 
uids outside these limitations and beyond low energies. 
For this purpose, the representation of the linear LL the- 
ory in terms of the fermionic quasiparticles introduced 
by Mattis and Lieb (1965) turns out to be more beneficial 
than the standard bosonization treatment. The quasi- 
particles share many features with their counterparts in 
the Fermi liquid theory. If the constituent particles of 
the liquid have a nonlinear dispersion relation, so do the 
fermionic quasiparticles. Similarly to the Fermi liquid 
theory, the interactions between the quasiparticles are 
weak. That, in principle, should allow one to build a full 
kinetic theory valid at low energies. The difference from 
the Fermi liquid theory comes in the relation between 
the measurable degrees of freedom and the quasiparti- 
cles: the corresponding transformation is rather nonlin- 
ear. While a nonperturbativc kinetic theory based on 
the fermionic quasiparticle representation has not been 
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developed yet, some elementary relaxation processes in a 
nonlinear LL are understood. 

The realization of links between the physics of a non- 
linear LL and the Fermi liquid theory makes available an 
arsenal of methods existing in the latter. One of them, 
the theory of the Fermi edge singularity, facilitated the 
development of new methods for the evaluation of the sin- 
gularities in the dynamic response functions of ID quan- 
tum liquids. The new paradigm which emerged is the 
description of the many-body dynamics in terms of ef- 
fective models of quantum impurities moving in linear 
LLs. 

We reviewed the existing tools for the investigation of 
nonlinear LLs and some results obtained with these tools. 
As we mentioned, building a kinetic theory of ID liquids 
remains an open question. Other questions closely re- 
lated to the review include the kinetic theory of weakly- 
nonintegrable systems relevant for cold atomic gases, and 
the dynamics of edge states (chiral and helical) relevant 
for electrons in solids. The field remains wide open be- 
yond these few problems, with a variety of practically 
and conceptually important questions to be answered. 
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